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Abstract 


An experiment in two-dimensional objective analysis, conducted at The Institute for Ad- 
vanced Study is described. The method consists of fitting a second degree polynomial by the 
method of least squares to the meteorological data in a limited area around each point in a 
rectangular grid, and then interpolating for the height of the pressure surface at the grid point. 
The results of machine computations give analyses which are comparable in quality to the 
subjective analyses ordinarily encountered. When used in a three-parameter atmospheric 
model for a numerical forecast, the objective analyses led to a more accurate forecast than care- 
fully prepared subjective analyses did. Also a satisfactory representation of the field of vertical 
stability (over North America) was obtained by the two-dimensional analyses, made for three 
different pressure surfaces. 

The greatest obstacle encountered in producing a completely objective analysis arises from 
errors of computation and transmission of raw meteorological data. This obstacle has not been 
overcome, and necessitates a preliminary subjective examination of the data for detection of 
errors. 


J. Introduction 


Ever since the meteorology group was set 
up at Princeton with numerical weather 
prediction as one of its major aims, the prob- 
lem of preparation of the initial data has 
been under consideration. It was realized 
that in the forecasting models the initial data 
would be required at a set of regularly spaced 
grid points whereas the observing stations 
have an almost random distribution. 

One method of getting appropriate grid 
values would be via the normal process of 
analysis, namely by subjectively drawing 
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isopleths and then by subjective interpolation 
obtain the values at grid points. This method 
has in fact been used to prepare the data 
presently being used. The preparation of this 
data, however, involved an experienced synop- 
tic meteorologist in several days of work 
preparing data for one observation time. In 
the preparation of the data considerable care 
was taken to ensure continuity both in space 
and time and to remove any erroneous obser- 
vations. This work would be very difficult to 
compress into a time compatible with any 
routine forecasting system. Two questions 
then arise. Is so much care necessary? Does 
the subjective analyst do anything that a 
computing machine could not also do in 
much less time? 
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The only way to find answers to these 
questions is to propose alternatives to the 
methods used by the subjective analyst and to 
use these alternatives in the preparation of 
forecasts. This paper proposes one method 
of avoiding the subjective analysis and dem- 
onstrates its effects, firstly, on the type of 
analysis produced, and secondly, on numeri- 
cal forecasts obtained from data so prepared. 

The method used has evolved during the 
last few years as a result of informal discussions 
held at the Institute for Advanced Study 
(I.A.S.)!. Initially an attempt was made to 
fit polynomials to the observed data over 
fairly large arcas. Such an attempt is described 
by PanorskY (1949). This method was not 
used in conjunction with a numerical forecast. 
Rather, as a result of certain discrepancies 
being revealed by the Panofsky investigation, 
the method described below was evolved. 


Here we are primarily interested in ob- 
taining an accurate representation of the data 
and not in smoothing. The only smoothing 
we require is that necessary to remove the 
random errors in the observations thus ensur- 
ing continuity in the horizontal derivatives 
of pressure. We seek a “most probable” 
representation of the data. If one attempts to 
fit a polynomial over a large area one auto- 
matically implies the use of a smoothing 
technique. If, however, we follow Charney’s 
suggestion and only consider data immedi- 
ately surrounding the grid point in question, 
we imply a minimum of smoothing. This 
procedure agrees with ordinary synoptic 
practice in that it does not pay any attention 
to data in remote parts of the region as would 
any method involving a Fourier analysis 
or polynomial fitting to the whole of the 
region. 

Dr. J. SMacorinsky (1954) of the U.S. 
Weather Bureau has performed independently 
a considerable number of hand computations 
using a slight variant of the method described 
here with comparable accuracy. Smagorin- 
sky’s work, however, was definitely limited 
by his lack of a high speed computing machine, 
a facility which fortunately was available to 
us during the present experiments. It was 


1 The other participants were Professor J. von Neu- 
mann, Mrs M. Smagorinsky, Drs J. Charney, G. Platz- 
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felt at the I. A.S. that an attempt such as 
described below should be made for the 
following reasons. Firstly, if the attempt 
succeeded, we would be relieved of a lot of 
work in future data preparation. Secondly, 
it was thought that the basic two-dimensional 
interpolation process should be tested in a 
series of machine computations in order to 
discover its capabilities and limitations. It was 
further felt that any further extensions of 
machine computations of analyses should wait 
until this basic process had been fully tested. 


II. Method 


The method is purely two-dimensional 
and is based on the fitting of a second order 
polynomial to the observations in an area 
surrounding the grid point at which the value 
of the pressure or height is required. In this 
paper we will be concerned with the height 
of given pressure surfaces expressed in terms 
of height deviation, D, from standard atmos- 
phere. The polynomial can be expressed in 
the form: 


D=a9+40% + 4o1 Y + A1 XY + 420 X? + 402 y” 
or in condensed form: 
i+j=2 


D= 2 a; x'y 
ai) 


(1 
Thus, if the origin of the coordinate system 
is taken to be the grid point under considera- 
tion, the D value we seek will be 450. In 
order to determine the six coefficients aj we 
clearly require the equivalent of six independ- 
ent observations of D at known positions 
(x, y). 

The observations give us wind values as 
well as D values. The wind values are incor- 
porated in the scheme by using the geostrophic 
assumption. Thus, one wind observadon 
gives us values of 0D/dx and 9D/dy at a 
point. Provided that there is at least one D 
value within the region being considered, 
each wind observation is equivalent to two 
more D values or pieces of information. 
The errors introduced by the geostrophic 
assumption are not random. However, due 
to our lack of knowledge of the geostrophic 
deviations, we have had to assume that 
the errors from this source are random. 

Since, in general, the observations from 
nearby stations are not in exact agreement, 
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we usually require more than the minimum 
number of pieces of information and then 
fit the second order polynomial to these 
observations by the method of Least Squarest. 
Thus we choose the a; such that we mini- 
mize 


I : I : OD 
De D, =D)» +— a 
a Pal pe ox 
7 oD, \2 . oD, 2 aD, Si 
Ox dy dy | 


where D,, 2D,/0x, 0D, [dy are from the obser- 
vations and the corresponding primed quan- 
tities are those computed from the polynomial 
(1). op is the standard deviation of the errors 
in the r observations of D and o,, is the standard 
deviation of the s wind observations in the area 
over which the polynomial is being fitted. 
The ratio op/o, is called the weighting factor. 
It can be shown that, if the errors have a 
Gaussian distribution, by choosing the a, 
in this way we get the most probable value 
of the required D. We thus have to solve the 
set of six linear algebraic equations: 


Ish 


= "C0 
da; à 


for the coefficients a; and in particular the 
coefficient ago. These equations can be written 
in the form: 


y’a=z 


where y and z are a (6x 6) matrix and a 
6 component vector respectively, the elements 
of both being known from the observations 
and a is a vector whose elements are the 
unknowns 4;;. 

By triangularizing the matrix y we can 
reduce the equations to the form: 


Y 


in which y* is a lower diagonal matrix. 

bts a0, 2%, are the leading ele- 
ments of y*, a and 2* respectively we obtain 
directly that 


x. x 


4a=? 


eg u Po ih 
D=499=20' /Yoo - 


1 A description of this method may be found in most 
books on statistics or numerical methods, e.g. Whittaker 
and Robinson’s “Calculus of Observations”. 
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Fig. 1. Location of 19 x 19 grid of points used in compu- 
tations. 


II. Weighting Factor 


It was not possible, as in some physical 
experiments, to look up the standard deviations 
of the errors in the observations in standard 
references, as the analysis must include data 
determined by a variety of different instru- 
ments and methods of observation. The error 
in measuring wind, for example, depends 
on the wind speed, and varies from one place 
to another over the map. Errors also arise 
from the introduction of the geostrophic 
assumption in the derivation of the expression 
for S. Consequently, the appropriate value 
of the weighting factor was determined by 
the trials described below. 

A series of objective analyses was computed 
for the soo mb surface for November 25, 
1950, at 0300 GCT. The grid used is shown in 
Figure 1. Each analysis was made with a dif- 
ferent weighting factor. The results of each 
analysis were then compared with the mean 
of three subjective analyses for the same map, 
made by three different analysts. The results 
considered were the D values for 222 points 
and Y?D, the finite difference Laplacian of 
these D values, for 152 points. The points used 
were selected prior to the experiment as being 
those where sufficient data were available to 
permit a fair comparison. The means of the 
standard deviations of D and y 2D of the ob- 
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Fig. 2. Standard deviations between objective and 
subjective analyses for various weighting factors. 


jective analyses vs. the subjective analyses are 
shown in Figure 2. It is evident from Figure 2 
that the objective analysis has the bestagreement 
with the subjective analyses, so far as height 
is concerned, when wind receives a relatively 
small weight. However, a heavier weighting 
on wind produces a better agreement when 
V 2D is considered. It was believed best to 
use a weighting factor giving good results 
for vorticity calculations, in view of the fact 
that the objective analyses are principally for 
use in numerical prediction. It may be noted 
that in numerical prediction one is usually 
more interested in expressions such as \ 2D 
and other derivatives of D rather than the 
absolute value of D. Consequently, a weighting 
factor corresponding to a ratio of so ft. height 
error to 10 knots wind error was selected for 
routine use in upper air analyses. 

The above type of test was not performed 
for other levels, but a smaller weighting factor, 
corresponding to a 30 ft. height error for a ro 
knot wind error, was selected for use in 1,000 
mb analyses. This decision was made because 
of the existence of greater errors in the geo- 
strophic assumption near the ground. At 
high levels wind observations are more 
accurate with regard to direction than to 
speed. This would suggest that in future 
analysis schemes we should use separate weight- 
ing factors for the two parts of a wind obser- 
vation. 


IV. Search Procedure 


An analyst, when drawing an isobar, looks 
at all the wind and pressure data in an area 
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surrounding the small region in question. In a 
similar manner, the objective analysis program 
fits the quadratic surface as closely as possible 
to all wind and pressure data reported within 
a prescribed area surrounding the grid point. 
The limits of the area are defined by the re- 
lation 
lx] + lyl= lal 

where x and y are the coordinates of each 
piece of data with respect to the grid point at 
which D is to be computed. This defines a 
square (standing on one corner) with the grid 
point at the center. The search then consists of 
locating and assembling in the memory of 
the computer all the data reported within 
this area. These data are then used in the 
fitting of the quadratic surface and the com- 
putation of the D value at the grid point. 
The search is made for each grid point using a 
given sized area. When sufficient data are 
assembled for the computation, the D value 
is computed. If the data are inadequate in 
quantity, the point is passed and a search is 
made around the next point. In this way, the 
computation proceeds until all points in the 
grid have been considered. This leaves some 
points computed and some not computed. 

The next step is to go back to the points 
which were skipped over because of lack of 
data, and to try to assemble more data for the 
computation. This is done by assembling all 
reported data plus all previously computed data 
which lie in the search area. The search area 
may be enlarged also. Thus, some of the 
points which were skipped the first time the 
grid was scanned can now be computed. 
(Points previously computed are now passed 
over and left unchanged.) Some uncomputed 
points are still left uncomputed, due to lack 
of data. After the machine has proceeded 
through the grid a second time, the search 
area is enlarged again and the above process is 
repeated. The search routine which is currently 
in use has evolved through many trials and is 
summarized in Table I. 

The sudden increase in the size of the 
search area from the 4th to the sth search is 
due to the fact that by the end of the 4th 
search the only areas remaining unanalyzed 
are those with practically no data at all. 
Searches 5 through 7 are designed to fill in 


these areas with values which are not un- 


reasonable. 
Tellus VI (1954), 4 
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Table I 


Size of search , 
area (length Data included 
of side of in search 


square = a\/2) 


Passage 
through grid 


Ist I,ooo km Reported data only 
2nd 1,000 km 

d 2 x 
oth bee En Reported data plus 
5th 1.800 km |{Previously computed 
6th 1,800 km |[4ata 
7th 1,800 km 


The square search area was selected because 
of the great speed with which it can be used 
in the computation. However, it has become 
evident that with a rectangular grid, one 
cannot expand a search area in too large 
increments from one passage through the 
grid to the next, or discontinuities in the anal- 
ysis develop. On the other hand, when re- 
ported data are scarce, very small increments 


of a square search area will prove unsatis- 


factory because no significant increase of data 
will be observed except when the increasing 
boundaries of the search area all at once 
pass a number of grid points containing 
previously computed data. Thus, it now 
appears that a circular search area would be 
meteorologically superior to a square one. 
The number of additional grid points enclosed 
with the expansions of the area would then 
increase more smoothly than with a square 
search area. 

Strictly speaking, one needs six pieces of 
meteorological information (each wind gives 
two pieces of information and each height 
value gives one) to fit a quadratic surface in 
the manner described above in section II. 
However, in view of the fact that the data 
contain errors, a surplus of data is highly 
desirable in order that some smoothing can 
be done. Another reason for requiring surplus 
information is that, with some combinations 
of data, a nearly singular matrix y is en- 
countered in the computation of the D value. 
A truly singular matrix has not yet been 
encountered in any of the analysis runs, but 
nearly-singular matrices have been seen. 
These give solutions for D which are com- 
puted in the non-significant digits of the data 
and are therefore meaningless. These can 
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never be entirely avoided, but the probability 
of their occurrence can be reduced to a negli- 
gible figure by requiring enough data. Ex- 
perience to date has indicated that this desired 
result can be obtained if a minimum of ten 
pieces of information is required as a prereq- 
uisite to the computation of the D value. 


V. Machine Program 


The machine program can be outlined in 
the following manner: 


1. Part A code is read into the computer at 
the start of the analysis. 

2. Part A reads ın data cards, one by one 
converts the data from decimal to binary, 
and stores it. The data cards contain the 
index number of the station (latitude and 
longitude for ships), the wind direction, 
and speed, and the height of the pressure 
surface to be analyzed for each reporting 
station. Part A then reads in part B code. 

3. Part B code reads in the station locator 
cards, which have the latitude and longitude 
for the corresponding index number of 
cach station which can possibly report. 
It matches this information with the data 
already in the memory and converts the 
index numbers of the data to latitude 
and longitude. All data are next converted 
to the coordinate system of the grid. Part B 
code then reads the code for parts C, D, 
and E into the machine. 

4. Part C code is the search routine. Whenever 
enough data are found to enable the com- 
putation of a point, part C code transfers 
the control to part D. 

s. Part D is the matrix code, i.e. the section 
which triangularizes the matrix y, etc. It 
computes the D value at the grid point, 
using the data obtained in the search. It 
then transfers the control back to the 
search code. 

6. Part E is the code which converts the re- 
sults of the analysis from binary to decimal 
and punches the results on cards, which 
can then be tabulated. 


It can be seen that the program is entirely 
automatic. The process of making the objec- 
tive analysis thus consists of the following steps: 

1) The reported meteorological data are 
punched on cards. 
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2) The cards are run through a printer to 
produce a listing of the data. This listing is 
then inspected by an analyst with the aid of 
a plotted map of the data for the detection 
of any serious errors in the data. This step is 
unfortunately necessary because of coding and 
communication errors. However, it need 
not delay the procedure unduly, since the 
maps can be plotted at the same time as the 
cards are being punched. It is not difficult to 
postulate an addition to the machine program 
which would check the data and reject those 
containing serious errors. This has not been 
put into practice, since it appears that such a 
procedure would lengthen the machine com- 
putations to an unreasonable time. 

3) After any serious errors detected in the 
data are corrected, the data cards are added 
to the code cards, the resulting deck is placed 
in the computer, and the computation is 
started. 

4) The machine computation then proceeds 
entirely automatically. At the present, slightly 
under 11/, hours are required for the compu- 
tation of the D values on one pressure 
surface for a 19 x 19 grid of points. It is known 
how to rewrite the code with one slight 
internal change which would reduce the 
above computation time to under one hour. 
Eventual reduction to a time of less than 35 
minutes is visualized for the IAS computer. 

5) The output cards from part E are usually 
printed to enable an inspection of the analysis. 
These cards are also used directly as input 
information for the forecast problem. The 
printed results at each grid point consist of: 

a) the D value obtained from the analysis, 
in tens of feet, 

b) a two digit number which tells the 
number of pieces of meteorological infor- 
mation used in the computation of the D 
value, and ; 

c) a digit indicating at which passage 
through the grid the computation was made. 


VI. Results 


The subjective analyses available for com- 
parison with the objective analysis of soo mb, 
25 November 1950, 0300 GCT, were prepared 
by three different individuals, and will be 
referred to as A, B, and C. Analyses A and B 
were prepared very carefully by experienced 
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analysts, who took all the time they needed 
for their completion. Analysis C was made 
hurriedly under operational conditions. These 
analyses were then compared with an objective 
analysis (O) for the same map, using the 
weighting factor corresponding to fifty feet 
height error and ten knots wind error. The 
standard deviations of D and V?D in feet 
(for the same points referred to in section III) 
of the difference between the various analyses 
are shown in Table II. 


Table II 


Analysis pairs A-B A-C B-C O-A O-B O-C 
Std. Dev..D:l.… 50:11:67 ı 60% 600472 WE 
Std. Dev. V* D113) Ifo 1062018008 


It can be seen from this table that the objec- 
tive analysis did not agree with the subjective 
analyses quite as well as the subjective analyses 
agreed with each other. However, the rela- 
tively close agreement between A and B is 
something which cannot ordinarily be ex- 
pected between subjective analyses. It should 
not be supposed from these results that the 
objective analysis was inferior to the subjec- 
tive analyses, since in this case the forecast is 
the ultimate test of the worth of the analysis. 

Previously a 24-hour numerical forecast had 
been made on the IAS computer for the in- 
tense cyclogenesis of November 24—25, 1950 
starting from 1,500 GCT, November 24. The 
input data consisted of the 1,000 and soo mb 
charts prepared by analyst B (above). The 
same forecast was then made from objective 
analyses. The two resulting forecasts, ex- 
pressed in the form of 24-hour forecast 
changes of height of the 1,000 mb surface, 
were then correlated with cach other yielding 
a correlation coefficient of 0.96. Both fore- 
casts predicted with fair accuracy the cyclo- 
genesis which developed, neither forecast 
being observably superior to the other. From 
this experiment it was concluded that where a 
good data coverage is found, e.g. over North 
America, two-dimensional objective analyses, 
prepared for two different levels, are satis- 
factory for use in numerical prediction. 

When three pressure surfaces in the vertical 
are to be analyzed for use in a forecast, a 
new problem is encountered. This is the 
problem of representing accurately the hori- 
zontal variation of vertical stability. This 
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variation is the additional feature treated in the 
three-parameter forecasting model, as com- 
pared with the two-parameter model. In a test 
of the representation of stability by the above- 
described objective analysis scheme, maps were 
prepared for 1,500 GCT, November 24, 1950, 
of a quantity directly related to vertical stability, 
namely 


A?D = Do — 2Dz00 + Diooo 


where the subscripts indicate the pressure 
surface at which the D value was read from 
the maps. This quantity was computed from the 
data at the three pressure surfaces at each grid 
point from the objective and also from the 
subjective analyses. From the grid point values, 
values of A?D were interpolated at the loca- 
tions of all reporting radiosonde stations. The 
values of A?D were also calculated from the 
actual radiosonde reports. These were used as 
verification of the values obtained from the 
objective and from the subjective analyses. 
The values of /\?D obtained by objective 
analysis when correlated with the values 
from the actual observations yielded a corre- 
lation coefliıcient of 0.90. The correlation 
coefficient between the values obtained from 
subjective analysis and the reports was 0.94. 
From this it was concluded that the two- 
dimensional scheme of objective analysis 
gives an adequate representation of the hori- 
zontal variation of vertical stability where 
the data coverage is good and when the 
pressure surfaces analyzed are at least 300 mb 
apart. An improvement over these results 
could probably be obtained with a more 
elaborate scheme of objective analysis. It is 
doubtful, however, whether or not the slight 
improvement which might be expected would 
justify the extra computation involved. 

The next experiment consisted of using the 
maps prepared for 1,500 GCT, 24 November 
1950, in two forecasts with a three-parameter 
model (3)—one to be made from the subjective 
* analyses and the other from the objective anal- 
yses. The results, at 400 and 1,000 mb, of 
these forecasts are presented in Figures 3 and 4. 
The 700 mb results are very similar to these. 
These figures show the surprising result that 
the forecast made from the objective analyses 
verified better than the forecast made from 
the subjective analyses. Further confirmation 
of this is given by the correlation coefficients 
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of forecast change vs. observed change at 
400 mb, where the difference between the two 


forecasts was greatest. These are given in 
Table II. 


Table III 
Correlation 
coefficient 
Forecast change vs. observed 
change—subjective analysis. 0.63 
Forecast change vs. observed 
change—objective analysis..... 0.79 
Forecast change—objective anal- 
ysis vs. forecast change—sub- 
jective analysis); Ges. ro 0.82 


The results described above can be attrib- 
uted to the fact that when the objective anal- 
yses are made, the geostrophic approximation 
is forced on the analyses, thus giving better 
results when a quasi-geostrophic model is used 
to make the forecast. This forcing is most 
marked in the vicinity of pronounced troughs 
or lows at middle and high levels in the 
troposphere, where the speed and curvature 
of the flow are strong. It takes the form of 
diminishing the analyzed geostrophic wind 
from the actual geostrophic wind in the areas 
where the observed winds are markedly sub- 
geostrophic. This is due to the fact that ob- 
served winds as well as all the heights of the 
pressure surface are used in the objective 
analysis. The smoothing process incorporated 
into the objective analysis also contributes to 
the same effect, particularly when the half- 
wavelength of the wind pattern is not much 
larger than the dimension of the search area. 
The quasi-geostrophic prediction model, using 
the geostrophic approximation for measure- 
ment and advection of vorticity, then has a 
more useful representation of these quantities 
than it otherwise would. 


VI. Future Extensions of Objective Analysis 


A logical step in the extension of the above 
analysis scheme would be to use a forecast 
made from an earlier map in the analysis 
over areas where reported data are sparse or 
non-existent. This can be compared to the 
activities of an analyst who consults the 
previous map when analyzing an area which 
has poor data coverage. The program which 
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Rigen: 
ment of the low center. 


has been drawn up (and partly coded for the 
computer) to accomplish this objective goes 
through the following steps: 

1) The analysis is started in the manner 
described earlier. However, after the search 
routine has been completed with the use of a 
certain area the analysis is stopped, leaving 
unanalyzed the areas where data are very 
sparse. 
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400 mb initial, forecast, and verification maps. The dashed line indicates the observed 24-hour move- 


2) The forecast D values for the pressure 
surface under consideration are read into the 
computer. 

3) At all points where the analysis was 
completed, differences (A D) between ana- 
lyzed and forecast D values are obtained by 
subtraction. At the points where the analysis 
was not completed, no differences can be 
formed. 
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Fig. 4 1,000 mb initial, forecast, and verification maps. The dashed line indicates the observed 24-hour move- 
ment of the low center. 


4) The field of AD is filled in at the un- 
completed locations by fitting a quadratic 
surface to the completed points around the 
uncompleted ones and interpolating for the 
uncompleted points. This completes the AD 
field. The values of AD are then added to the 
forecast D values, giving a “corrected” fore- 
cast D value at each point where no analysis 
could be made. 
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5) The corrected forecast D values are 
given a certain weight and treated as data. 
The original analysis code then goes over the 
reported data, previously computed data, 
and the weighted and corrected forecast 
values. It fits the quadratic surface to all these 
values in the vicinity of the unanalyzed point 
and interpolates for the height of the pressure 
surface at the point in question. 
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VIII. Conclusions 


The results obtained above lead to the 
conclusion that a simple two-dimensional 
scheme of objective analysis is feasible and 
satisfactory for the purpose of transforming 
raw data into input data for a two- or three- 
parameter numerical prediction model. In 
fact such a scheme appears to be better suited 
for use in a quasi-geostrophic prediction 
model than information obtained from sub- 
jective analysis. The principal difficulty en- 
countered so far is introduced by incomplete 
and erroneous meteorological observations. 
The detection and correction of the errors 
at present must be done subjectively, i.e., a 
skilled analyst must inspect the data very 
carefully for omissions and errors. One could 
argue that this process practically amounts 
to a subjective analysis. It should be pointed 
out that in the preparation of input data for a 
forecast from a subjective analysis, the anal- 
ysis itself is only the beginning of the process. 
The reading, recording, and checking of 
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interpolated information at the grid points is 
the most tedious and lengthy part of the 
process. All this is unnecessary when an objec- 


tive analysis method is used. It is clear, how= 
ever, that the reduction and eventual elimina- 


tion of errors in computation and transmission 


of meteorological data is a subject which must 


receive high priority in future research and | 


development before a truly objective analysis 
is possible. 
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volves a minimum of multiplications. 


1. Introduction 


One of the objections frequently raised 
against the recently developed methods of nu- 
merical weather prediction is that they require 
calculation of second and third order deriv- 
atives of quantities that are measured only 
approximately and at widely separated points. 
Virtually all of those methods are based on 
vorticity equations for ‘“quasi-geostrophic” 
flow—i.e., vorticity equations in which the true 
wind and vorticity (but not divergence) are 
approximated by the geostrophic wind and 
the corresponding “geostrophic vorticity’. 
Thus, the horizontal gradient of geostrophic 
vorticity appears in these equations as a third 
derivative of the height of an isobaric surface. 
Coupled with the observation that the per- 
centage errors of finite-differences generally 
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A Note on the Integration of the Vorticity Equation for 
Quasi-Geostrophic Flow 
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Abstract 


. By the transformation of integrals arising from the solution of the “quasi-geostrophic” 
equations for barotropic and simple types of baroclinic flow, it is shown that the isobaric height 
tendency can be expressed in terms of first derivatives of the heights of isobaric surfaces. This 
demonstrates that the procedure by which the equations are solved must tend to cancel or 
“average out” the large percentage errors in the geostrophic vorticity advection—i.e. errors in 
second and third order finite-differences. The height tendency may be computed by a simple 
graphical procedure, which is carried out directly on analyzed contour charts and which in- 


increase as the order of the difference increases, 
this fact tends to show that the horizontal 
advection of geostrophic vorticity at individual 
points cannot be computed accurately from 
widely separated measurements. Accordingly, 
there has been considerable doubt that the 
quantities which enter into the “quasi-geo- 
strophic” vorticity equations can be calculated 
exactly enough for purposes of numerical 
weather prediction. 


One should suspect, on the other hand, 
that the inaccuracy of computed finite-differ- 
ences at individual points is not as serious as 
it might appear from an error analysis of the 
unintegrated equations. This stems from the fact 
that the solution of the quasi-geostrophic vorti- 
city equation can be expressed as an area (or 
volume) integral of the advection of absolute 
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vorticity or, stated in slightly different terms, 
as a weighted average over a considerable 
area (or volume). Thus, if the errors of the 
computed geostrophic vorticity advection at 
individual points are essentially random, they 
will tend to be “averaged out” in the process 
of solving the vorticity equation. For this 
reason, and because most errors in reported 
observations are nonsystematic, the percent- 
age errors in the solution of the vorticity 
equation are generally much less than the 
percentage errors in the computed vorticity 
advection at individual points. 

The quasi-geostrophic vorticity equation 
is not, of course, usually solved by a method 
of direct integration. For reasons of economy 
and simplicity, the most widely used method 
of solution is the so-called “relaxation” 
method. It is evident, however, that the proce- 
dure of “relaxation” also tends to smooth out 
nonsystematic errors in the computed geo- 
strophic vorticity advection. It can be shown, 
in fact, that the Richardson method of relaxa- 
tion is essentially equivalent to FJORTOFT'S 
(1952) method of repeated averaging. This 
suggests that the large nonsystematic errors 
introduced by second and third order differ- 
ences are generally “averaged out” by what- 
ever procedure one uses to solve the vorticity 
equation and, accordingly, that the percentage 
error in the solution is comparable with the 
errors of first order differences. This con- 
jecture would be confirmed, if the solution 
of the quasi-geostrophic vorticity equation 
could be expressed in terms of first derivatives 
of the height of an isobaric surface, thus 
showing that the errors in computing higher 
order derivatives must be cancelled out in 
the process of solving the vorticity equation. 

The purpose of this note is to show that the 
solution of the vorticity equation for quasi- 
geostrophic flow is expressible in terms of 
first derivatives of the isobaric contour height. 
For simplicity and clarity, we shall first 
consider the quasi-geostrophic equation ‚for 
barotropic flow, later indicating how the 
results in that case can be extended to certain 
simple types of baroclinic flow. Finally, we 
shall outline a graphical method for com- 
puting height and thickness tendencies from 
the original height and thickness charts. 
Although this method is not adapted to 
machine computations, it will be found 
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useful for manual calculation of tendencies at a 
small number of strategically located points. 


2. The Equations for Quasi-Geostrophic, 
Barotropic Flow 


In accordance with the usual procedure, 
the quasi-geostrophic vorticity equation is. 
derived by first eliminating the horizontal , 
velocity divergence between the equation of | 
continuity and the vorticity equation, and by 
replacing the true wind and vorticity with the | 
geostrophic wind and geostrophic vorticity. | 
At the same time making use of the hydro- 
static equation, we may put the vorticity ! 
equation for barotropic flow in the form: 

7 7 d 

< viy+J(p vy) + Er =o (1) 
in which t is time; x and y are rectangular : 
coordinates directed locally toward the east: 
and north, respectively, and measured in ı 
planes that are locally tangent to the geopoten- : 
tial surfaces. 


y is gzd-}, the geostrophic stream function 
in isobaric surfaces; z is the height of anı 
isobaric surface; g is the gravitational accelera- 
tion; À is the Coriolis parameter; 


_ du dv 
Ox oy 


2 2 “3 
N ee 4 
dy 


1068 dy ax’ 


R is the gas constant; and T, is a representative : 
value of the absolute temperature at the ground | 
surface. All lines normal to the (x, y) plane: 
have been treated as verticals, and all differ- - 
entiations with respect to x, y, or fare carried ! 
out with pressure held fixed. Stated briefly, . 
the problem is to integrate Eq. (r), beginning ! 
with known initial values of y at every point ! 
in the (x, y) plane. 

The key to the usual method for solving Eq. | 
(1) lies in regarding it as a nonhomogeneous 


linear equation, in which a is the dependent 


variable, and in which the terms not involving! 
time derivatives are regarded as known. 
Viewed in this light, Eq. (1) takes the form| 
=) = F(x, y, t)= known function (2) 
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_ Now, since the function F involves space 
derivatives only, its initial value at every 
point can be computed from the given initial 


values of y. Thus, the initial values of = at 


C 
every point can be computed by solving Eq. 
(2). Having computed the time rate of change 
of y at each point, we next extrapolate linearly 
from the initial values of y, to obtain the value 
of y at every point in the (x, y) plane and at a 
time slightly later than the initial moment. 
Finally, the predicted values of y are regarded 
as a new set of initial data, from which F 
and the y-tendency can be calculated by 
repeating the procedure outlined above, and 
from which one can predict the values of y 
at a still later time, and so on. 

The principal variations of this procedure 
arise from the use of different methods for 
solving Eq. (2). The method that will be 
applied here is an exact analytic method 
which is discussed in almost any elementary 
text on theoretical physics. As it turns out, 
equations of the same type as Eq. (2) frequently 
appear in the theory of forced vibrations of 


an elastic membrane, and their solutions are 
d 

well known. If we let = correspond to the 
C 

displacement of a membrane, » to the wave 

number of the membrane’s free oscillations, 

and F to a periodically impressed force, the 

theory of elastic membranes can be applied 

to the solution of Eg. (2) without essential 


modification. 


d Mi 
ihe NES Vee tp)? | 
-F(E, n, t) dédy (3) 


where £ and 7 are dummy variables corre- 
sponding to x and y, respectively, and Ky 
is the zero order Bessel Function of the 
second kind with imaginary argument. The 
indicated area integration extends over the 
entire (£, 7) plane. 


3. Transformation of the Vorticity Advection 
Integral 
The difficulties discussed earlier arise in the 


calculation of F, a process that is necessarily 
carried out by the method of finite-differences. 
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5 ; 0 
In particular, the computation of = from 
0 
Eq. (3) requires evaluation of a weighted 
integral I of the geostrophic vorticity advection, 


namely, 


T= JG ()J(v, v°y) dA (4) 


where G (r) = _ K, (vr); r is the distance 


between the fixed point (x, y) and the variable 
point (£&, 7); and dA is an element of area. 
The domain of integration A extends over the 
entire plane. Our present concern is to trans- 
form the integral I in such a way as to reduce 
the order of the derivatives of y. By making 
use of standard vector identities, we may 
rewrite I in the form 


T= [v-(GKx£vy)dA- fvG-KxlvpdA 
A A 


where £=v?y and K is a unit vector directed 
upward. According to Gauss’ Theorem, the 
first of the integrals in the equation above can 
be transformed into two line integrals, taken 
around a circle of infinite radius with center 
at (x, y), and around a circle of vanishingly 
small radius with center at (x, y). Both of 
these line integrals vanish, however; the first. 
because G vanishes exponentially at infinity: 
and the second, because y has continuous 
derivatives in the neighborhood of G’s singu- 
larity at (x, y). Thus, from the preceding 


equation, 


oo 


(dG ow 
=f {5 G4 dOdr (5) 
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where © is the angle between the x-axis and a 
line connecting the points (x, y) and (£, n). 
Written in this form, I is to be interpreted 
as a weighted integral of the radial transport 
of vorticity. 

As given by Eq. (5), I still involves second 
order derivatives of y. To reduce the order 
further, we shall next express & (the Laplacian 
of y) in terms of derivatives with respect to 
the polar coordinates (r, ©). Bearing in mind 
that G depends only on r, we may then rewrite 


Eq. (5) as 
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oo Qn 


(dG [T1 oy 
I= | 7 flaca 


0 0 


10 { op\ | oy 
ror (r +) | 20) de 
on 


” 1 dG f 2 f2y\? 
= | 2r? ef als) m 


0 0 


CO Qn 
1dG (9 [ dyp\ dy 
Je 5 or ( 4 10 0 (6) 
0 0 


In the first integral on the right side of Eq. (6), 
the integration with respect to © can be 


: P : ow . : 
carried out directly. Since =, is continuous, 


the first integral vanishes. The second may be 
integrated by parts, as follows, 


oo Qn 
ıdG f 2 { yoy 
f lt or 5) ES 
0 0 
27 


AG yo (aw 
À | Fal or or (35 ee 
0 


0 
27 CO 
N 1 dG 9 [ dy dy 
CE [fe dr or (2 3) dde = 
0 0 


CO Qn 
‘1 dG [9 (ay\? 
a u 
0 0 


In the second integral on the right side of 
Eq. (7), the integration with respect to © can 


: 4 Ow . : 
be performed directly. Since 7 is continuous, 
r 
the second integral vanishes. The first integral 
on the right side of Eq. (7) is again integrated 
by parts, this time by integrating first with 
respect to r. 


co 


27 
3 d (dGdy dy 
5 if or ( dr or 2) Hae 
0 


0 
27 © 
"fd (1 dG\ dy dy 
- ff in; E Z) or 00 drdO (8) 
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In the first integral on the right side of Eg: (8); 
the integration with respect to r can be carried 
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: dG 
out directly. That integral vanishes, for 7 


vanishes exponentially at infinity, and the 
derivatives of y are continuous in the neigh- 


borhood of the singularity of z at r=o. 


Finally, the integral I may be put in the form 


oo 27% 
er ap oy 
Te - #0 [2 RR) 
0 0 | 


d (1dG 
where H (r) miler (: i 
that Eq. (9) expresses I as a weighted integral 
of the radial transport of tangential (angular) | 
momentum, illustrating the connection 
between the principles of conservation of 
angular momentum and vorticity. 
Substituting the result of Eq. (9) into Eq. 
(3), we may now write the solution of the 
vorticity equation for quasi-geostrophic, baro- 
tropic flow as 


oo 27% 

init Dati, ap op 

M(x, y) = - [ro [FE doar 
0 0 


. It is noteworthy 


oo 27% 


+ | HE (| Bd (10) 


The points to be emphasized about this for- 
mula are that G (r) and H (r) are known ana- 
lytic functions which can be computed with 


any desired degree of accuracy, and that 7 
is expressed in terms of first order derivatives 
of y. This demonstrates that the large percent- 
age errors incurred in computing higher 
order derivatives of y must be cancelled out 
in the process of solving the vorticity equation 
for the y-tendency. As conjectured earlier, 
the percentage errors in computing the y- 
tendency from Eq. (10) are much less than the 
errors in the computed geostrophic vorticity 
advection at individual points. 


4. The Quasi-Geostrophic Equations for 
Simple Types of Baroclinic Flow 


Our next aim is to show that the solutions 
of the quasi-geostrophic equations for certain 
simple types of baroclinic flow can also be 
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expressed in terms of first derivatives of the 
heights of isobaric surfaces. If this is possible, 
then whatever conclusions one may draw 
regarding the errors in the solution of the 
barotropic vorticity equation can be extended 
to baroclinic flow as well. 

Either by direct specification or through 
finite-difference approximations, the “two- 
parameter” baroclinic models developed by 
Eapy (1952), ELIASSEN (1952), THOMPSON 
(1953), SAWYER and BusxBy (1953), and CHar- 
NEY and PHirres (1953) all contain the implicit 
assumption that the direction of the vertical 
wind shear (or the direction of the horizontal 
temperature gradient) is independent of height, 
and are all governed by equations of the 
same general type. For the present purposes, 


it will be convenient to deal with a system of 


equations whose mathematical form is similar 
to that of the barotropic vorticity equation. 
In an earlier report (THOMPSON, 1953), the 
author has shown that the equations for quasi- 
geostrophic flow in which the isotherm 
patterns in different isobaric surfaces are geo- 
metrically similar can be put in the form 


D he un in ob 
“v2 72 ear el tape 
av et (y, voy) tasty, V2) + Po ol) 


2 925 + J (p, 729) +] (9, 729) + BI (@, v9) + 


ot 
(12) 


Ap aT ian @ Po 
+ Bo ee, 


where the “barred” quantities are vertical 
averages, integrated with respect to pressure 
through the entire depth of the atmosphere; 
=RTA 1, T is the absolute temperature; a 
and b are nondimensional constants that are 
determined from a representative vertical 
wind profile; and y? is a constant that depends 
on a representative measure of the static stabil- 
ity through a very deep layer. The quantity 
y is the geostrophic stream function for the 
vertically averaged wind. Similarly, p is the 
geostrophic stream function for the difference 
between the vertically averaged wind and 
the wind at the ground surface. That is, 
owing to the hydrostatic relationship between 


y and y, 
(13) 


g = Y- Vo 
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where y5=gzA"! and 2, is the height of an 
isobaric surface next to the ground. 

It will be noted that the only terms of Egs. 
(11) and (12) which involve derivatives of 
second and third order, which do not involve 
time derivatives, and which do not already 
have the mathematical form J (x, V2x) occur 
in the combination 


Je, ve) +] (pvp) 
With the aid of Eq. (13), however, the ex- 


pression above may be rewritten as 
T (yp, V2y — 72) +] (9, v?p+ V2y) 


=J(y, 7%) +] (9 V9) - (yo, V°40) (14) 


The fact to be stressed is that the terms of this 
expression do have the mathematical form 
J («, v?x). Since the transformations by which 
one passes from Eq. (4) to Eq. (9) are valid 
for any scalar variable y, area integrals of the 
terms of Eq. (14) can be transformed to inte- 
grals of the type shown in Eq. (9). Finally, 
because Eqs. (11) and (12) are of the same 
form as Eq. (2), it follows that 2 and an 
be expressed in terms of first derivatives of y 
and y. This argument shows that the solutions 
of the quasi-geostrophic equations for simple 
types of baroclinic flow are not subject to 
the large percentage errors of second and 
third order finite-differences. 


5. A Graphical Method for Computing the 
Vorticity Advection Integral 


The most laborious and time-consuming 
step in the computation of the y-tendency 
from Eq. (10) is the evaluation of the integral J, 
given in Eq. (9). We shall next outline a 
simple graphical procedure whereby the 
integral I can be evaluated directly from ana- 
lyzed isobaric contour charts, with a mini- 
mum of numerical calculation. The first 
step of this procedure is to draw in the con- 
tours of the isobaric surface (e.g., 500 millibar 
surface) for unit intervals of y. The second 
step is to superimpose on the contour chart a 
transparent overlay, on which are drawn a 
set of concentric “‘coordinate” circles spaced 
dr apart, as illustrated in Figure I. The common 
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center of these circles is placed on the point at 
which the y-tendency is desired. Another 
set of “reference” circles is drawn midway 
between the “coordinate” circles. 

Thus far, the area of integration has been 
divided into a set of nested rings of width dr. 
The next step is to subdivide each ring into 
smaller elements of area, bounded by radial 
segments passing through the intersections of 
the contours with the appropriate “reference” 
circle and subtending angular increments dO. 
A typical element of area is the shaded area S 
in Figure I. It remains to find the contribution 
from each element of a given ring, the net 
contribution from a given ring, and finally, 
the total contribution from all rings. 

A fact that simplifies the procedure consider- 
ably is that Eq. (9) can be written in the form 


It 


CO 27 
fie - [#0 [ dapdow 
0 0 


where d,y is the variation of y over a radial 
displacement dr and doy is the variation of 
y over an angular displacement dO. Fixing 
attention on the typical area element shown in 
Figure I, we see that the increments dO have 
been chosen in such a way that doy is always 
unity, since the contours are drawn for unit 
intervals of y. Thus, the element of the ©- 
integration has the numerical magnitude of 
dy, the difference between the values of y 
at the midpoints A and B of the circular seg- 
ments which bound S. The sign to be affixed 
to that difference is positive or negative, 
according as the angle y at which the contours 
cross the circles is acute or obtuse. The ele- 
ment of the r-integration is obtained simply 
by adding up all the differences d,y for each 
ring (summing with regard to sign) and by 
multiplying each sum by the appropriate 
value of the weighting function -H(r). Finally, 
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the integral I is evaluated by summing up the 
contributions from all rings. Since H(r) 
decreases rapidly with increasing r, the #= 
integration need not be extended to infinity 
in order to insure that the error is less than 
any preassigned limit. 

It is worth noting that this procedure in- 
volves only as many multiplications as there 
are rings and, accordingly, that it can be 
carried out very quickly with only an adding 
machine. This method is not, of course, 
suited to automatic machine computation, 
owing to the irregular size of the area ele- 
ments. It is, however, ideally adapted to 
manual computation of y-tendencies at a 
limited number of isolated points—e.g., points 
on trough- and ridge-lines and inflection 
points. 


Tellus VI (1954), 4 | 


CHARNEY, J. G., and Pris, N. A., 1953: Numerical 

Integration of the Quasi-Geostrophic Equations for 

Barotropic and Simple Baroclinic Flows. J. Meteor., 

10, pp. 71-99. 

Eapy, E. T., 1952: Note on Weather Computing and 

the So-Called 21/, Dimensional Model. Tellus, 4, pp. 

157 — 168. 

ELIASSEN, A., 1952: Simplified Models of the Atmosphere, 
Designed for the Purpose of Numerical Weather 
Prediction. Tellus, 4, pp. 145 - 157. 

FJORTOFT, R., 1952: On a Numerical Method of Inte- 
grating the Barotropic Vorticity Equation. Tellus, 
4 pp. 179-195. 


Tellus VI (1954), 4 
 2—410071 


IIERORATTONZOEZTITESVORTICETLVZEEOUMTION 


325 


REFERENCES 


Sawyer, J. S., and Busugy, F. H., 1953: A Baroclinic 
Model Atmosphere Suitable for Numerical Integra- 
tion. J. Meteor., 10, pp. 54-59. 


Tuompson, P. D., 1953: On the Theory of Large-Scale 
Disturbances in a Two-Dimensional Baroclinic 
Equivalent of the Atmosphere. Quart. Jour. R. M. S., 
79, pp. 51 — 69. 


THompson, P. D., 1953: Appendix I, Report of Progress 
for the Period 1 Feb. 53-1 Aug. 53, Joint GRD - 
AWS Numerical Prediction Project. AF Cambridge 
Research Center, Cambridge, Mass. 


Chemical Denudation in Sweden 


By ©. ARRHENIUS, Grödinge 


(Manuscript received June 9, 1954) 


Abstract 


Many analyses of groundwater have been made by the Swedish Health authorities. Of those 
about 10,000 have been used for an attempt to map the geographic distribution of various 
substances. The factors analyzed are: Total amount of inorganic substances, Ca, Fe, Mn, NH,, 
SO, CI, NO,, HCO;, PO,, pH and oxidation by permanganate. It is found that the water 
of wells with a depth less than 2 m are more diluted than that of the deeper ones. 

Scania and the plains of the later district have rich waters whereas the Highlands of South- 
ern Sweden and the northern part of the country (Norrland) have water with a low salt 
content. Drilled wells are richer than the shallow ground water wells. The chemical denuda- 
tion is much stronger than formerly calculated. One can estimate that the topmost metre of 
the soil has lost about 1/5 of its weight by denudation since the ice age. Great losses are also 
caused through the sewage waters which are carried away to the sea, by the burning of straw 
etc. 

Very much may be done in order to stop or to counteract the chemical denudation. A well 
planned crop rotation, planning of a forest culture which improves the soil conditions, taking 
care of all affalls are some of the roads leading to a decreased chemical denudation of the soil. 


The geochemical processes within the various 
parts of the cycle of the earth are so intimately 
related to each other that it is hardly possible 
to say whether an action takes place in the 
lithosphere or in, for instance, the hydrosphere. 
This applies above all to the part played by 
water in weathering. This action also greatly 
affects conditions within the biosphere, so 
greatly, in fact, that one is justified in saying 
that organic life is wholly dependent upon the 
direction weathering takes. An occurrence to 
which little attention is paid is that which 
takes place within that part of the biosphere 
which might be called the anthroposphere. 
Human beings interfere with and regulate 
nature, and their activities have become so 
potent that in certain cases the results may be 
compared with geological events. 

In the present work the author attempts to 
give a cross-section of the reaction known as 


chemical denudation, caused by percolating 
water, subsoil water and watercourses, and the 
contribution of man to these processes. 

Rain falls, water sinks into the ground, and 
all the substances in the soil are more or less 
dissolved. Soil moisture, the source from which 
all plants derive their nourishment, is formed. 
Salts are liberated and precipitated as the water 
sinks farther down. At greater or lesser depths 
the water collects as subsoil water which in 
places comes to the surface as springs, wells, 
watercourses and lakes. Subsoil water carries 
with it dissolved substances and deprives the 
top soil of part of its stock of nutrients. This 
is what is known as the chemical denudation 
of the land. 

This impoverishment of the land has up to 
now been measured by examining the amount 
of dissolved substances carried by rivers into 
the sea (9, 15, 21). This is in itself a very im- 
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portant phenomenon as it consists of a trans- 
port of substance from the lithosphere to the 
hydrosphere. But more important from a bio- 
logical point of view is the transport of dis- 
solved substances from the upper layers of 
soil to those below. 


To illustrate these problems we have an 
extensive material consisting of analyses of 
water collected at our chemical stations and 
by the State Institute of Public Health. The 
material at our disposal includes water from 
deep wells, lakes and watercourses, water- 
pipes, shallow wells and springs. It is only the 
two last named which in the present work are 
taken as being representative of subsoil water. 

Analysis instructions are to be found in 
Meddelande fran Kungl. Medicinalstyrelsen, 
1950. (31). 

This material suffers from one drawback, 
viz. the samples were in the main taken from 
populated areas. Thus these districts rather than 
the great uninhabited areas are represented. 
Some samples, however, were taken from 
unpopulated areas. In addition there is a very 
extensive material consisting of conductivity 
determinations from Vattenfallsstyrelsens’ main 
power-lines through Sweden. By means of 
these investigations a belt taxation of the soil 
conductivity within large parts of the unin- 
habited areas has been carried out, the results 
of which are in good agreement with the 
maps based on water analyses. 

The results are given below in mg/l. As far 
as Ca is concerned, this is calculated from the 
total hardness, but the result is given in 
mgCa/l. 


ZAHN. 
Depth 
The samples of water were taken from wells 


of different depths. 

The deepest wells are to be found in the 
clayey plains (Upsala, Kristianstad, Malmöhus, 
Skaraborgs län) and in Kopparbergs län and 
Västerbotten. The least average depths are 
found in Blekinge, Halland and Värmland. It 
is to be expected that the salt content of the 
water should vary with the depth. Table ı 
shows the variations in water from wells of 
different depths. The total salt content, as well 
as Cl and Ca, is doubled in the first two 
metres. 

The consumption of permanganate, on the 
contrary, is greatest in shallow wells, probably 
on account of the relatively greater amount of 
organic substance in both the water and the 
surface soil. 

It is of great interest to ascertain whether 
the relationship between the different ions and 
the total content of inorganic matter is constant 
or not. 

In Table 2 the total salt content has been 
fixed at 100. 

From this it can be seen that the relative 
composition of the salts is practically the same 
atrditterent depths as fan as SO, ClyCa, Fe, 
Mn and NO, are concerned. On the other 
hand, the water in wells two metres deep or 
shallower seems to be more diluted than that 
in deeper wells. 

One explanation of this may be that the 
water on its way down through the earth 
attains complete saturation only after a certain 
period of time. It may also be that the upper 


Table 1. Depths of shallow wells and the content of different substances in the water by mg/l 


h in Total salt 5 Permanganate 
et I content | a | oS | He | Ma NOs consumption 
I 5 118 13 Bae, 0.7 0.57 3.6 21 
2 4 722 13 4.1 1.0 0.08 5.4 20 
3 5 194 26 6.2 122 0.15 9 25 
4 7 199 2 6.6 1.6 0.21 13 27 
5 a 231 32 7.5 => 0.18 14 19 
6 8 237 29 8.2 1.7 ORT, 14 18 
7 9 224), 29 Gh ey 0.18 LUI 2 
8 10 245 32 7.5 12 0.23 9 20 
9 9 235 36 Feo 1.0 0.09 12 15 
10 9 211 23 6.5 52 0.30 6 18 
II 15 209 19 OF, 0.8 0.12 10 18 
12 8 226 21 8.5 0.7 GES 12 9 
Wee ES 
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Table 2. Figures in Tab. I referred to total salt 


content = IO00. 


Depth Total 
in SION salt C1|Ca| Fe | Mn | NO, 
metres content 
| 
I 4 100 MNT 00005) 
2 3 100 II] 34| 0.8] 0.07] 4 
3 3 100 7373210.0:012.0:071 5 
4 3.5 100 731023 80.8 Ce mc 
5 3 Ioo 14| 32| 0.5] 0.08] 6 
6 3.5 100 121034 2.0:512.0.07|0506 
7 4 100 13| 34| 0.5| 0.08] 5 
8 4 100 13| 31| O.5| 0.09] 4 
9 4 100 15| 30| 0.4| 0.04] 5 
Io 4 100 sus O5 
II 5 100. 9| 32| 0.4| 0.06| 5 
12 35, 7100 9| 37| 0-4] 0.06| 5 


layers are as a rule so exhausted that the soil 
moisture in them is weaker than that in the 


deeper layers. 
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The Factors Examined 


The factors usually investigated when water 
is analysed are sulphate, total dissolved in- 
organic substance, called below total salt 
content, chlorine, bicarbonate, Ca, iron, man- 
ganese, nitrate, nitrite, ammonia, phosphate, 
permanganate consumption (oxidizing capa- 
city) and pH. As the samples have, in some 


cases, to be transported rather long distances | 


one is justified to a certain extent in doubting 
the value of the pH determinations. 

Not all the samples, however, have been 
examined with regard to the above mentioned 


substances. In some cases sulphate determina- | 
tions have been omitted, in others Mn, and 


so on. 

In Table 3 all the figures have been grouped 
according to the total salt content and averages 
have been made for each group and factor. 
For each of the factors SO,, Cl, bicarbonate, 
Ca, Fe, Mn, NOs, NO,, NHs3, permanganate 


Table 3. Total salt content compared with the content of other substances. The figures give mg/l. 
| 

ar so, | cl |Hco,| Ca | Fe | Mn | No, | No, | NH, |*2™| po, | pH 
SERIE I 1.6 10 6 0.3 | 0.06 ©. Moxon fe) 25 0.06 5.9 
DAS g 2 2.5 15 9 SE || Orn 0.0 ox 0.1 20 0.01 5.9 
Seo Dans Wr 2 3.6 20 10 0.6 | 0.03 || NOT O.I 18 0.04 6.0 
ARE CO ONO'S 3 7 24 12 0.4 0.03 PES, 0.1 0.4 20 O.I 6.1 
N Re ES 5 6 30 15 0.9 | 0.08 2EQ) MET O 0.2 15 0.4 6.2 
eg 5 6 34 19 0.9 | O.II 4.011 207 O.I 16 0.07 6.1 
2 5 8 39 21 | TO Ol O7 0.2 20 0.3 6.1 
SO sistas adie 6 9 40 24 11.23 || OKO 5.5 0.9 0.4 16 O.I 6.3 
PR sates: 6 II 51 32 0.8 | 0.04 8 3 O.I 14 0.7 6.4 
ITR ER ee 6 13 63 34 1.5 0.15 8 4 0.9 24 0.6 6.6 
TA HINTERE 6 ne 81 42 1.5 0.14 9 4 0.5 BD, 0.5 6.7 
LEO wt custome 8 18 95 51 | sehr II 4 Te 24 0.4 6.5 
TOR Le 8 19 107 59 18310 || One 13 4 142 22 0.5 6.6 
DOCS F0 E voa 8 20 146 WP 17.2, | Poured II 6 0.6 24 O7 6.6 
DORA. 8 23 156 Ti 122 ORTS II 3 1.3 2D 0.6 6.5 
ZOOFSTERERER 8 23 182 OI 1.5 | 0.16 II 4 1.4 22 0.6 6.8 
DIT NE ae 9 25 211 104 | Coasts) nr Io 12 BE 0.3 6.9 
300 ee Io 25 223 III 0.8 0.15 13 6 163 22 0.8 6.9 
LÉ en as 14 34 22 120 3 0.32 16 7 32 22 0.9 6.9 
DHCP ee Io 29 266 126 CT 0.83 II 4 Beil 30 1.0 7.0 
SOL ee 9 38 279 134 A 0.23 13 Gr, 2.4 28 12 6.9 
BOO IE: 23 38 294 139 2x0 |. 0418 14 7 3.3 19 af it 7.0 
ADS ete see ices 25 41 258 146 1.4 | 0.20 2 2 1.8 26 Dey) 7.0 
VE rentes 10 60 293 151 ROMEO 19 18 Set 28 2.0 6.9 
AS An teen ati Io 82 266 141 O.9) ||| 0887 17 35 6.8 23 0.8 7.0 
SESES eke caves 23 64 312 154 2.0 170.20 12 te 4-4 2; 0.8 6.8 
BO ge 2 92 295 159 ZA 0.68 21 39 18 38 DEN 7.0 
7302028 eme 42 176 341 192 1e) 0.47 29 15 9 25 12 Yoke 
1,620 : 30 644 324 308 3.0 1.09 25 3 6 34 Dek Pod 
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Fig. 1. Sulphate and total salt content. 
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Fig. 2. Chlorine and total salt content. 
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Fig. 3. Bicarbonate and total salt content. 


consumption, PO,, and pH a curve has been 
plotted showing the relationship between the 
total salt content and the respective factor. 
Fig. 1 shows the relationship between the 
total salt content and sulphate. It is seen that 
the SO, content increases in proportion to the 
salt content to between 300 and 400 mg/l. 
. After that the curve becomes very irregular. 
Fig. 2 shows the relationship between the 
total salt content and the chlorine content. 
The figure shows that a linear relationship 
exists up to the salt content 4co—450, after 
which the curve changes its direction abruptly 
and the chlorine content increases more rapidly 
than before. This is probably due to the basic 
substance, calcium carbonate, increasing its 
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solubility very slowly at this concentration, 
and to the fact that the increase of dissolved 
salts is due to an increase of chlorides. 

Between the bicarbonate content and in- 
organic substance (Fig. 3) the relationship 
seems to be that the increase in bicarbonate 
content takes place somewhat more rapidly 
than that of inorganic substance. 

Ca (hardness), Fig. 4, shows a slight decrease 
with a high total salt content. This is in agree- 
ment with what has been said above con- 
cerning Cl-total salt. 

A comparison of the Ca/Cl relationship 
shows that the chloride content increases more 
rapidly than the calcium content. Gripenberg’s 
(23) values for the waters of the Gulf of 
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Fig. 4. Calcium and total salt content. 
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Fig. 5. Relation between calcium and chlorins in ground water and Bothnian and 
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Fig. 6. Fe in relation to total salt concentration. 
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Total salt 


Fig. 7. pH and total salt content (Units as in fig. 6). 


Bothnia and the Baltic are included for com- 
parison. From this we perceive that subsoil 
water contains relatively more Cl than sea- 
water at high concentrations (Fig. s) while 
conditions are contrary at low. 

The connection between the total salt con- 
tent on the one hand and Fe and Mn on the 
other is not particularly marked. This is 
illustrated in Fig. 6. Other conditions, such as 
pH and oxidation-reduction potential etc., are 
decisive here. 

No connection between the total salt con- 
tent and the permanganate consumption can 
be traced, but on the other hand there exists 
some relationship between Fe and Mn. 

There seems, however, to be some connec- 
tion between acidity (pH), so that the greater 
the acidity of the water, the lower is the salt 
content (Fig. 7). If this is due to the acid 
soils being poorer or to their exhaustion must 
be investigated later. 

The relationship with PO, seems to be 
good, but there are so few analyses that the 
results should be used with caution. 

The averages indicate that up to 300—400 
mg/l total salt content the composition of the 
water is the same, although diluted in varying 
degrees. 


Geographic Distribution 


Maps have been compiled of the different 
characteristics SO,, total salt content, Cl, Ca, 
Fe, Mn, NO, and permanganate consumption. 
In order to facilitate the work somewhat, 
parish averages have been calculated and these 
form the basis of the compilation of the maps. 
(The maps are found at the end of this issue.) 

The picture of the total salt content obtained 
shows that the country can be divided into 
four regions: Norrland, the Central Swedish 
agricultural area, Smäland and Skäne. (Map 1.) 

The Norrland region is bounded in the 
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south by a line drawn from the mouth of the 
River Dalälven to the southern part of Lake 
Vänern. The country to the west of this lake 
connects to the Smäland region. There are a 
number of small districts in the Norrland 
region with a rather high total salt content. 
The largest of these is the district round Lake 
Storsjön. 

The Central Swedish region shows great 
variation, and the plains proper have high 
contents. 

The Smäland region is bounded roughly in 
the north by a line running from east to west 
through the southern part of Lake Vättern 
and to the south by a line drawn from Hal- 
landsäs through Eslöv and Knislinge to Mör- 
rum. The coastal areas in the east and west 
have somewhat richer water. In the inland 
districts the western half is considerably poorer 
than the eastern. 

The Skäne region has exceptionally rich 
water, particularly in the coastal areas. 

The distribution of the calcium content is 
on the whole about the same as that of the 
total salt content. In the north the districts 
round Lake Storsjön and Lake Siljan are out- 
standing. Only in small areas is the calcium 
content in Halland high. (Map 2.) 

There is a striking difference between the 
provinces of Upland and Södermanland. It is 
of interest to compare this picture with that 
which Lunpauisr’s investigation give of the 
calcium content of the clays in the district of 
Stockholm (3). There, too, the southern part 
of the district with its deficiency of calcium 
is in sharp contrast to the northern part, which 
has a particularly high calcium content, due 
to the Gävle limestone. 

The geological explanation of the variations 
in calcium content can be arrived at only 
after a careful study of all the details. 

The water of Skäne is rich in nitrates (Map 6) 
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but, except for a small increase in the Öster- 
götland plain, their occurence diminishes to- 
wards the north. The production of nitrates 
may be assumed to be very dependent on the 
climate. Romell has collected figures from 
the State Forestry Research Board which show 
that the farther north one gets the lower is the 
nitrogen accumulation in woodland soil. The 
nitrate and nitrogen producers are clearly de- 
pendent upon the temperature of the soil and 
the length of the growing season. (Oral 
information.) 

The ammonia content of water seems to be 
in co-ordination with that of nitrate. It cannot 
be proved that nitrogen in water is in any 
way inherent to coastal arcas. This material, 
therefore, gives no support to the theory that 
nitrogen is chiefly a cyclic element produced 
by the sea (34, 41). 

Curiously enough, the consumption of per- 
manganate (Map 8) is greatest along the Baltic 
coast. These waters seem thus to be more 
deficient in oxygen than those in other parts 
of the country. An explanation of this may 
be that the coasts of the Baltic are flat 
(stagnating waters). 

The content of iron and manganese (Map 7 
and $) is also higher in these parts. The subsoil 
waters of Skäne are usually poor in manga- 
nese. 

The occurrence of sulphate (Map 3) is not 
following the distribution of calcium. It is true 
that Skäne is extremely rich in sulphates but a 
relatively high sulphate content is also charac- 
teristic of southern Smäland. The Central 
Swedish region has high values too, extremely 
so in certain alum-slate areas. It is of interest 
to note that the two most northerly provinces 
show a somewhat higher sulphate content than 
the two lying immediately to the south (cf. 
the alum soils frequently occurring there). 

Finally there remains the chlorine content 
of the subsoil waters (Map 4). As is the case with 
sulphate, Skäne (tertiary and other sediments 
in the moraine) is also rich in chlorine, and 
the area rich in Cl stretches into southern 
Smäland. The chlorine content is also rather 
high in the Central Swedish region. Particu- 
larly high contents are to be found round the 
southerly point of Lake Vänern. 

The waters of Södermanland are richer in 
chlorine than those of Upland. 

The coastal area of Halland has subsoil 
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water considerably richer in chlorine than the 
waters of the inland districts. The chlorine 
content diminishes very rapidly with the 
distance from the coast. 

Häruın has (31), in connection with corro- 
sion investigations, made analyses of rain- 
water on the west coast of Sweden. 

The results of this investigation are shown in 
Table 4. 


Tab. 4. The occurrence of Cl and SO, in rain- 
water and its conductivity (according to Härlin). 


Total 
Distance so Sp. inorganic 
4 
from sea g/l | mg/l conduc- substance 
in km mg} 2 tivity (calc. by 
the author) 
oO 22. 13 LOS TON 116 
2 5 9 71 49 
I 1.8 5 46 32 
135 0.7 6 29 20 


Similar results have been arrived at in 
Holland. 

This demonstrates how rapidly wind-borne 
salt decreases towards the interior of the 
country. 

Only ro per cent of the samples have lower 
or the same chloride content as two kilometres 
from the coast. If for the rest of the country 
we calculate with 0.7 mg Cl/l it will be found 
that this figure is a few per cent of the amount 
found in subsoil water. It seems hardly prob- 
able, therefore, that the chlorides in subsoil 
water are carried from the sea by wind. 

The following part will deal with deep 
waters, i.e. those obtained from the bedrock 
by boring. There we shall see that where the 
bedrock is rich in chloride, the subsoil water 
is also rich in chloride. It must therefore be 
assumed that there exists a rather good circula- 
tion between subsoil and deep water, and that 
this is the most important source of the chloride 
content of the soil. 

It might be thought that urine and manures 
would affect the chloride content of subsoil 
water. This, of course, often occurs purely 
locally. In very many cases the surface water 
runs into the water supply, but generally 
speaking the addition of chloride emanating 
from human beings is so small that it cannot 
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be of any significance. Every human being 
emits about 5 kg Cl per annum (46). Thus 
the whole population of Sweden should pro- 
duce about 35,000 tons of Cl. Our import of 
potassium salts is the equivalent of about 
60,000 tons of CI (39), making a total of about 
100,000 tons. This is equal to 2.5 kg per hectare 
a year. This would, assuming that about one- 
third of the precipitation became percolating 
water, be the equivalent of about ı mg/l. For 
the country asa whole the addition is negligible 
It might be considered that the high chloride 
content of Skäne water is due to the density 
of the population and the intense use of 
manures. It can be assumed that a farmer with 
one hectare of land, growing sugar beets and 
kitchen produce, would use 200 kg potash. 
His family consists of five persons. Together 
this gives 125 kg Cl/hectare and year or about 
8o mg/l in the subsoil water. It is only in 
such extreme circumstances that the amount of 
CI found locally in normal Scanian subsoil 
water can have been contributed by human 
agency. 


Water from drilled deep wells 


Water from drilled deep wells gives informa- 
tion of the composition of water in mainly 
unweathered parts of the crust of the earth. 
The influence of cyclic elements, organisms 
and so on is much reduced at these depths. 
The type of bedrock must be of primary 
importance. The material obtained from the 
State Institute of Public Health included anal- 
yses of water from about 2,000 deep wells. 
A compilation has been made in the form of 
maps (Maps 9—14). On account of the limited 
extent of the material the results have been 
collected by “härads” and the average cal- 
culated. 

This material shows that all the ions with 
the exception of NO,, NO,, NH, and PO, 
occur much more richly than in subsoil water. 
In order to give a better survey of the rela- 
tionships a table has been compiled giving the 
averages by provinces for both deep wells and 
shallow wells (Table 5). Water from deep 
wells is much richer in chloride. Nitrogen 
compounds, on the contrary, occur in smaller 
amounts. It is impossible on the basis of this 
material to say whether they originate from 
percolating surface water or from the bedrock. 
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People of today seem to be of the opinion 
that lakes and watercourses were created so 
that all kinds of filth can be deposited in them. 
This is what is known as modern hygiene. 
What is commonly called filtered sewage 
contains large amounts of nutrients that disturb 
the balance for the flora and fauna of water- 
courses and lakes. At the same time large 
amounts of nutrients are transferred from the 
lithosphere to the hydrosphere, and are thus 
withdrawn from that part of nature’s economy 
upon which human beings are dependent for 
their livelihood. 


When modern hygiene has become uni- 
versal we can assume that all refuse from man 
will go into the watercourses. 


Sewage from the food-producing industries, 
the waste of the timber industry, the destruc- 
tion of straw and humus by fire must cause 
the loss of valuable substances. 

An attempt has been made to draw up a 
survey of how people in Sweden waste 
nutrients. 


Survey of losses and gains of nutrients in tons 
per annum. 


N | 120), || INO 
Refuse from human be- 

INS ee ee 28,000] 8,000] 7,000 
Household water (asses- 

Sete) tel. one ee 7,000] 2,700| 2,300 
Sitarchetaceores tte) CE 250 140 480 
Dairies (rinsing water) (20, 

Be 80.0 0 do cbr rade 238 60 70 


Welten) (GAN a Glos a oo ae oe 20—200| 8—50| 4—80 
Silo fodder (press water) 

(2 BO Svat Sek 130 60 410 
Bark on floated timber? 

(EOS teten 1,350 380 790 
Export timber and cel- 

InlosesToa nee 20,000| 2,000] 25,000 
If all straw is burned (36)| 21,000 


1 Assuming that all Swedish homes are fitted 
with W.C. and drainage. 

2 Timber is to an increasing extent barked be- 
fore floating. 


The export of grain and foodstuffs causes 
loss of nitrogen, potash and phosphorus. 
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Table 5. Averages for provinces and characteristics (D — deep wells: S = shallow wells). 


a nm ni ke 


Total 
County SO, salt CHE IEICO, | Ca: Fe Mn: |" NO, | NO, PNET eee, 
content 

we robe Miam er are He, ise ele ee eee ee 
Stockholm D 42 377 88 213 89 1224 I 0.73 4-3 5.3 1.5 | 0.54 

S 5 248 34 82 94 1022 |00.80 6 6 I 1.64 
Up Saar Eee. De 326 80 252 121 1.9 | 0.10 4-9 4-9 2.0 | 0.61 

S Io 293 31 236 1H TNT 0.0.76 7 7 4 1.10 
Södermanland...... D 42 349 70 202 85 PAS || OEP 8.8 5.7 1.720022 

S 15 DG 32 277 92 214] MOTO INTZ 6 2 0.91 
Östergötland....... D 55 374 48 245 124 1.4 | 0.18 9.5 TE 0.6 || 030 

SE 263 34 165 93 72220-202076 8 I is 
iOMKO pase year D15 170 7, 109 60 1.5 | 0.59 3.7 18.3} 1.6 | -0.0008 

SEO 132 19 40 34 0.5 | 0.08 18 2 I 1.63 | 
Kronoberg. ner. D 34 ea) 20 107 57 Tan |) :0:52 2.8 0.7 1.0 | 0.7500 

S 18 139 23 36 32 Ox7 :0.082 27 5 1.0 | 0.3700 
AMAR Re ere Dr5 361 93 141 76 2.6 | 0.19 BiB 1.5 1.72 Onin 

Sis 213 36 98 57 11892 0.68 Tel 3 3 0.89 
LILO. Wn a AO D 42 652 249 129 148 | GE ZX 1-7 | 71.0088 | 

S 46 404 109 79 88 0\6N|MO 154163 2 0.4 | 6.454 | 
Kristianstad.. on... 1D) 57 378 33 310 162 wo) | Koyasmr 15.6 1.9 0.4 | 0.65 

S 46 334 52 130 104 0.5 0.09 | 62 4 it I 
Malmohusiae. „Man. 1D) 2X0 549 103 355 159 5.4 | 0.33 6.2 1.4 4.5 | 2.50 

S 60 527 66 280 178 0.72 | oe 258 6 5 
Hallandesgı Eur. 2x: D 42 225 138 105 53 RS I Ko} 7.0 33 3.2 | 0.50 

S 15 236 68 80 54 1:92 20.472 0.74 4 6 1.39 
Göteborg & Bohus..| D 22 304 91 119 42 ets), LOSTT fe) 0.1 1.5 | 0.00 

Se 190 29 igen 43 2.2 | 0.38 6 57 ee I.II 
AIveborer eee ce D 34 ZI 29 108 62 0.9 | 0.45 5.7 2.4 0.5 | 0.26 

SP 16) 171 26 93 47 0.8 | 0.38 12 2 2 1.58 
Skaraborö.n..n..c.- D 42 680 313 279 139 TIM 0:20 8.0 4-5 4-5 | 0.62 

S730 292 61 185 94 TOOL II 2 I 1.70 
Varmlande"" D 18 237 45 156 71 1.12 15,0,29 2.3 2.5 0.6 | 0.25 

Sr 3 103 I4 61 28 || ACTES 6 2 I 0.25 
Nacken rer |e TS 363 85 180 96 1.4 | 0.28 Pops B® I.4 | 0.20 

S 10 163 17, 103 60 Os, Mes 2 T. 0.56 
Västmanland......| D 30 318 88 187 88 0:64 12:0,21 6.0 4.6 0.5 | 0.65 

S 10 94 20 129 70 1.820,29 9 4 I 1.00 
Kopparberg........ D 26 160 16 132 62 1.6 | 0.45 3.9 3.5 2.6 | 0.40 

SES 130 15 80 41 7.297 02544 17 5 2 0.2 
CENA ORS PP EEE EEE D 34 274 58 531 52 1.9 | 0.16 3.0 1.6 1.2. || 03 

Seo 119 18 58 31 280100 720 8 6 2 0.23 
Västernorrland..... D 15 187 19 115 52 3448057 0.2 1.5 IT 1.81 

SC 110 12 61 34 177472015 5 3 2 0.40 
émane rer ESS D 26 219 12 212 105 0.70 90:22 5.3 4.1 4.1 | 0.08 

Se 150 Io 128 64 0.7 | 0.08 8 I I 0.03 
Vasterbotten war: D 26 190 35 100 45 3.0 | 0.28 255 BZ 3.0 | 0.86 

SN 81 13 2 22 0.9 | 0.14 5 5 I 0.12 
Norsbotten yen eer D 10 202 2 137 48 2.6 | 0.18 1.0 97, 3.0 | 2.90 

Sagt 129 50 53 41 23503 7 4 6 0.50 

Import 31 tons per ann: This survey shows that the loss of nutrients; 
| x in industry is of subordinate importance (even! 
_N |PR& | KO | if it may be extremely important locally). 
| The great losses are in sewage and in agri- 
ONCAKCS (SO) ARE 900 340 — 1 att, 
Artificial manures (39) ...| 50,000] 90,000] 30,000 culture and forestry. 
(15,000 from the We waste more nitrogen than we import,’ 
atmosphere) 


we import seven times as much phosphoric 
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Table 6. Concentration of soluble inorganic substances, sulphates, chloride and nitrate in some Swedish 
rivers, and in the water which flows into the different rivers. 


No’ information about the nitrate in the rivers is available, but judging from other (28, 13, 38) in- 
vestigations the occurrence of nitrate in river water seems to be very low and is given in this table = o. 
It is easy to see that it is not only the total inorganic substance that decreases on the way from sub- 
soil water to the river mouth, but also the SO, Cl and NO,. 
a Er 


Totalt salt 
SO; ort Cll NO; |Ca+Mgl 


RARE UT 


4 river water 0.9 212 3.T fe) 5 

SSIES EH Se ONE EC ats 5 iG 3 se 
river water 1.5 20.9 3.3 fo) 6 

MÉCANISME aM en An nee F2 7 
u river water — 20 32 fo) 4 
RN © OR, Sennen eC ER Se Sa, 6 a = a 6 
Mapermanalven....<....-..+e.0s+s- | Te neater 5 ae 2-5 4 8 
Se Fe ee subsoil water 4 107 II 9 32 

4 river water Pail 32.6 4.5 o 15 
HEGMSAIVEn PAR PET Se aes een | le. 5 en : : AR 
: river water 1.2 28.4 3.0 fe) 12 
TAD Agrar ee ee ee | näher 5 a ” 2 a 
2 river water 1.6 24.6 2.5 fo) 8 
BASANT Re talc. «eos Severs es dou or abel angen oe AD x 54 n ê a8 
> river water 2.0 253 2.8 o 7 
(DIEU Zoo. ars SER EN ee ORR | a ter 28 Be ar = Ro 
ne river water 9.6 43.9 4.9 fo) Io 
NÉE SEE IMMO EE too ays er es 25 5 = > 
re | river water 12.7 58 bo Oe oO 22 
EAC ODI SS Tt acest ee he en er es 258 = co of 
ae river water 27.9 138.4 9.6 (6) Aus 
OT EST 22, «yh Las er 35 AE 58 28 naa 
: river water eG) 49.6 Gen o 20 
SEAMPAM See ses cs so sha 20 aie sie core à | ee 5 | 233 ae 17 83 
F river water gen 40.1 6.8 fe) 14 
Ae a a Se TR eas ltr 5 153 ot _ at 
= river water 5.3 30.6 7-9 o 8 
MÉTRO D MO Te Other ay 157 28 18 us 
river water 7-4 51 ep, fe) 14 

Bean ed ra | een 35 ats 33 38 28 
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acid as we throw away and a certain balance 
exists with regard to potash. 

And we consider that we can afford to 
waste one or two hundred million kronor a 
year for the pleasure of enriching the sea with 
nutrients. 

At the same time we destroy our water- 
courses with our waste matter so that the fısh 
die out, our water supply is contaminated and 
our opportunities for swimming are limited. 
In order to remedy some of these drawbacks 
we add chemicals to the water, which then 
becomes sterile. The equilibrium of our waters 
is disturbed, and they lose their ability to 
purify themselves. 

We have seen above how much the subsoil 
water carries away. On an average for the 
whole country every litre of percolating water 
removes from the upper layers of the soil: 
Totalt salt content SO; Cl Ca Fe Mn NO; PO, 


122 or A or GO) a 


If we calculate that one-third of the pre- 
cipitation evaporates, one-third runs away and 
one-third soaks into the ground, and an average 
precipitation of 600 mm it makes per square 
metre (200 ]): 
Total salt content 

26.4 


SOACIRCAREE 
1.87 226 


Mn NO, PO, 


84.0:2400:.0344 1.2 Onl 07 


For a square kilometre the figures will be 
the same but will then indicate not gr but ton. 
If we calculate that denudation began 10,000 
years ago and that it was not greater then than 
now (in reality it was probably 2—4 times 
as great then, and has since fallen to its present 
values), we arrive at the following figures for 
leaching (in 1,000 tons per square kilometre): 


Total salt content SO, Cl Ca Fe Mn NO, PO, 
264 ue) PAO) A One TA I 


OS ARR ERE INTIS 


A square kilometre one metre deep con- 
tains about 1,500,000 tons of soil. Since the 
Ice Age, therefore, at least one-fifth of the 
total weight of the soil has been washed away, 
and only the least soluble materials remain. 

For the whole country the annual loss 
should be (in 1,000 tons): 


Fe Mn NO; PO, 
575 41 


Totalt salt content SO, Cl Ca 
10,800 740 1,070 3,280 99 14 


These figures represent incredible amounts. 
The fertility of the soil is continually dimin- 
ishing. 


Leaching 


Stone recently ground down by the ice 
naturally produced a much more concentrated 
water-solution than does the present denuded 
stone-powder, the soil. As a minimum of the 
salt content of the soil-moisture in those days 
we must take the concentration of water from 
deep wells. On account of the great surface 
of finely-ground powder exposed to the action 
of water it was probably even greater. Leaching 
has continued since then. An example of how 
leaching progresses is given below. Different 
kinds of rocks were ground in a colloid mill 
and an attempt was made to extract phophorus 
acid from them. Extraction was carried out once 
by means of HCl, 3—6 times by 2 % citric 
and 3—6 times by 0.01% NaCl solution. The 
sum is given for the first two extraction 
mediums. The results of the various extractions 
is given for the water extractions. Averages 
are given for the different types of rocks. (5). 


From this it is clear that water dissolves only 
a small part of what citric acid and hydro- 
chloric acid can do in the same time. But the 


Extraction from different types of rocks (average) by means of HCl and citric acid, and fractionated 
with water (0.01% NaCl solution). 


The figures give: P,O, in 1/1,000 per cent of dry weight 


HCl 2 % citric Water 

acid leu eres ee a VI 
Gneiss RED 138 97 0.91 0.22 0.19 0.20 0.06 0.07 
Granite RL 145 129 0.63 0.39 0.52 0.24 0.II 0.05 
Limestone. .....= 353 20 0.22 0.12 0.41 0.15 0.0 0.0 
Sandstone...... 49 17 0.25 0.13 0.38 0.10 0.10 0.0 
SIA RE: 165 115 m2 0.15 0.45. 0.35 0.15 0.20 
Hälleflinta...... 125 94 0.15 0.0 0.15 — — == 
Porphyry nt Re 73 65 0.12 0.20 0.60 0.0 0.0 0.0 
Dieritee ner: 286 ZE 120 0.73 0.40 0.6 0.06 0.06 
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gradual decrease in the yield of phosphoric 
acid is of the greatest significance. 

In the previous part we gave a Table illus- 
trating losses in the upper layers of soil. The 
rivers carry much smaller quantities to the sea 
(6). Table 6 gives a comparison of the con- 
centrations of river water near the mouth of 
the river, and the subsoil water concentration 
in the precipitation areas of the rivers (32). 

No information about the nitrate in the 
rivers is available, but judging from other 
(28, 13, 38) investigations the occurrence of 
nitrate in river water seems to be very low and 
is given in this table = o. It is easy to see that 
it is not only the total inorganic substance 
that decreases on the way from subsoil water 
to the river mouth, but also the SO,, Cl and 
NO,. 

The total amount of salt washed out is 3—4 
times greater in subsoil water than in river 
water (Table 7). On its way to the sea the 
water must, by precipitation, absorption, the 
action of organisms and so on loss large 


Table 7. Relation Subsoil water: River water. 
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amounts. It is of interest to note that this 
applies also to such easily dissolved ions as 
SOM Clend-NOe 

Generally speaking it can be said that lime 
is precipitated most effectively in the large 
rivers of Norrland. Nyköpingsän, Stängän, 
Mörrumsän, Olman and Gullspängsälven also 
show great precipitation. 

The three rivers of northern Halland, 
Mjölnaän, Tommarpsän and Kävlingeän, and 
Tidan and Lidan in Västergötland, on the other 
hand, show only little precipitation. 

The shorter a watercourse is, the smaller is 
the difference between subsoil water and river 
water. 

As far as sulphateion is concerned, there is 
no more to be said than that Tidan, Lidan, 
Tommarpsän and Kävlingeän show high con- 
tent and little loss. The same applies to the 
inorganic substance. 

With regard to chloride, no conformity can 
be found, but it is surprising to find that this, 
too, diminishes so greatly when subsoil water 
becomes river water. Nitrate and phosphate, 
two materials so important to organisms, de- 
crease particularly rapidly, which is, of course, 
quite natural, for plants take up with great 
avidity every trace of these substances in the 
water. 

If we take into consideration the disap- 
pearance of the various substances from the 
subsoil water we find that the sulphateion 
undergoes the least change, while chloride, 
Ca + Mg and the total salt content disappear 
in about the same proportions. 

We see therefore that denudation is par- 
ticularly great. If only the conditions during 
the most recent decades are taken into con- 
sideration we find that about ro times as much 
lime and nitrogen disappear as is added to the 
soil. Equilibrium seems to exist with regard to 
phosphoric acid. 

As regards the removal of calcium, a calcula- 
tion was made on the basis of two extensive 
investigations, one made in 1930 and the other 
in 1940 (7). It shows that the sugar-beet area 
of Skäne, which on an average is given 120 kg 
CaO per hectare and year, has become 0.13 

H units more acid. If it is calculated that the 
soil of Skäne requires on an average 3,000 kg 
CaO to change the reaction one pH unit, it 
would mean that 60 kg would disappear per 
hectare and year. The plants take away 90 kg 
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and denudation about 100 kg. This would 
mean a loss of 190 kg and an addition of about 
120 kg. A deficiency of 70 kg as against a 
calculated 60 kg. 

If we consider the figures for leaching during 
10,000 years we find that 2,600 tons of easily 
dissolved nutrients have disappeared per hectare 
or, as mentioned above, at least one-fifth of 
the total weight of the top metre of the soil. 

The yield of a cultivated plant is usually 
considered to be an indication of its well-being 
and within small areas, where the climate can 
be ignored, it is the available plant food that 
is decisive (Map 15). A study of the map of the 
yield per hectare of winter wheat in Sweden (8) 
shows that it is identical with the one giving 
the distribution of calcium in the water. There 
is every reason for assuming a relationship 
between soil rich in plant food and greater 
yields (Fig. 8). 

If the yields of parishes, for which water 
analyses are available are compared with the 
calcium content of the water, it is found that 
there is excellent agreement (Fig. 8). This 
does not mean that the calcium content is the 
deciding factor with regard to the fertility of 
the soil, but rather a sign of how much the 
soil is impoverished. 

A similar distribution: occurring sparsely in 
the Norrland and Smäland regions and very 
frequently in Skäne and the Central Swedish 
region is found in most grove plants, our 
finest deciduous trees and most of the in- 
habitants of herbaceous slopes (örtbackar). 

The localities occurring sporadically in the 
Norrland region are to be found as a rule on 
what Gunnar Andersson called “south moun- 


Yield 
da 
30 
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tains’. HALDEN (24, 25), GRAPENGIESSER (22) 
and WisTRAND (48) have all shown that the 
name “south mountain” is misleading, for 
many of the localities where “relicts” are found 
slope towards other points of the compass, 
even towards the north. Their main charac- 
teristics are extreme weathering and water rich 
in nutrients. 

Gunnar Andersson has taken the retreat of 
the oak and the hazel towards the south as the 
basis of a theory on the deterioration of the 
climate. 

The present writer has shown that the lower 
parts of peat layers are richer in mineral 
nutrients than the upper parts of the same kind 
of peat (4). This indicates that formerly the 
peat was growing and deposited in richer 
water than is the case at present. 

Sandegren has been able to demonstrate that 
in certain deposits in Bohuslän only the oldest 
are rich in calcium. Sandegren has provided, 
orally, the following information: “During 
S.G.U.s peat survey, samples from Bohuslän 
often contained freshwater calcium mud in the 
deepest layers of the bog. All such occurrences 
were at a high altitude, and according to 
pollen analysis they belong to the oldest of the 
post-glacial periods, the Boreal period.” 

HESSLAND (27) has found similar cases. 

The retreat of certain plants towards the 
south may just as well be explained as being 
caused by a change in the factors of the 
habitat as by a change in the climate. 

In view of the fact that agricultural research 
has produced certain races of plants suitable 
for different climates, particularly winter 
wheat, one can say with regard to this plant 
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Fig. 8. The calcium content of the ground water and the yield of winter wheat in the 
lake district of Sweden. 
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that it is the soil factors that are responsible 
for its varying degrees of thriving. 

Thus the retreat towards the south and the 
present “climatic limits of plants” can no longer 
be advanced as a proof that the climate has de- 
teriorated. The proof must be sought elsewhere. 


What can be done to counteract denudation? 


Denudation affects most seriously countries 
with humid climates, that it to say the world’s 
best arable areas, and is therefore a great, if 
invisible danger, probably greater than the 
visible mechanical erosion (10). 

The question rises involuntarily: Can any- 
thing be done to counteract this destruction 
of the soil? If nothing is done, the soil of the 
country will gradually become waste. Instances 
of such soil devastation may be found in many 
places. Rome’s campagna, which, during the 
early years of the history of Rome was its 
granary, was neglected more and more, the 
soil was denuded and the pH value sank to 5. 
It was not until the present century that the 
land could be cultivated anew. In Scotland 
and England the forests were laid waste, arable 
land was turned into pasture for sheep with 
the result that in some parts the cultivation of 
the land again seems impossible, while in other 
parts it is costing enormous sums. During the 
sixteenth century in western Java flourished 
the Sultanate Bantam. The Deutch described 
the country as very fertile, but now the soil 
is exhausted and to a large extent covered 
with brushwood and the country continues 
to decline. 

This is what will happen in Sweden, although 
more slowly, if nothing is done to prevent 
impoverishment of the soil. 

One practicable way is to increase evapora- 
tion so that less water percolates through the 
soil. We have done this unconsciously. Our 
yields, both of agricultural and forests, have 
increased. These consume more water. BERG- 
STEN’s (11) investigations have shown that the 
discharge of our rivers has diminished con- 
siderably during recent years. In many places 
farmers are of the opinion that wells are run- 
ning dry and that the subsoil water is subsiding. 

This is a very dangerous development, for 
although we may prevent leaching, we may 
at the same time cause a general aridity. 

Another, perhaps more difficult but better 
and more certain solution is to seek help from 
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plants to a greater extent than hitherto. As is 
well known, plants have an incredible ability 
to collect the food necessary for them from 
very diluted solutions (26). Thus Burp (14) 
has shown that plants in both good and bad 
soil reduce the concentration of soil moisture 
to the same low level and keep it there. If 
the land is then allowed to lie without vegeta- 
tion the amount of dissolved substances in- 
creases rapidly. If plants are allowed to keep 
the land free from nutrients, the percolating 
water cannot take much away. 

During summer hardly any water soaks 
through the soil (no water runs through under- 
ground drains then), and when the ground is 
frozen no percolation can take place. Spring 
and autumn are the seasons when leaching 
takes place. 

Such crops as winter cereals, winter rape 
and pasture are some of the plants that might 
be used. Bare land should clearly be avoided 
during the winter. Investigations must be made 
to find out which ways are most practicable. 

A shorter period of rotation in the forests 
and more frequent thinning out would prob- 
ably be an effective way of preventing un- 
wanted leaching, perhaps it might even be 
replaced by an accumulation of plant foods and 
better nutrients conditions and an increase in 
fertility. It appears as if the number of needles 
produced by a spruce attains a maximum when 
the tree is so—60 years old, after which it 
decreases (19). If a tree is allowed to live 
60—70 years it need never pass through the 
later period, when its needles become fewer 
and it requires more nourishment, without 
being able to produce enough needles to 
enrich the soil. 

If a wood is allowed to remain undisturbed 
only dead needles will fall to the ground and 
become humus. When trees are felled living 
needles are mixed with litter and humus. 

EBERMAYER (19) has shown the relationship 
between the different kinds of needles. 


Living 


Dead 
Current year |Previous year 
% % % 
Total ashes . 1.55 1.89 1.52 
Ioacousoaot 0.62 0.41 0.14 
Ca ere 0.19 0.49 0.44 
ee Roc 0.30 0.24 0.06 
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The living needles contain about five times 
as much of the important plant foods potash 
and phosphorus as the dead ones. 

If we can grow trees with large crowns, 
the falling needles will be reduced to a mini- 
mum, and when the trees are felled the living 
needles will be incorporated into the ground. 
Birches and other deciduous trees contribute 
more to the soil than do pine trees. An ad- 
mixture of leaves is necessary to keep up the 
productivity of the soil. 

Trees derive their nourishment from the 
lower layers of soil, and are thus able to 
redeem nutrients that have sunk down below 
the reach of the roots of herbaceous plants. 

Finally it must be stressed that the ground 
must be kept occupied with trees in order to 
prevent leaching. 

It may be said, therefore, that forest, prop- 
erly cared for can do much to eliminate the 
harmful effects of chemical denudation. 

We must also cease our present waste of 
refuse. Modern techniques are moving in the 
right direction by trying to utilize all the by- 
products of manufacture. 

But town waste matter, too, must be made 
use of, for it represents great values, and at 
the same time the threatened contamination 
of our water can be avoided. 

This may be done by spraying meadow- 
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land with sewage, as is done at Lund. In this 
way the plant food is effectively saved, and 
the water continues on its way completely 
purified (44, 45). 

At Munich (35, 47) the Hofer system has 
been applied since 1912. After the sewage has 
been roughly filtered and diluted it is led into 
dams containing plants and fish. Bottom vege- 
tation develops rapidly and serves as food for 
the fish. The water runs outcompletely purified. 

For about seventy years so-called rieselfelder 
(47) have been in use in Germany. The sewage 
is led on to these fields where it sinks through 
the soil. After a certain time the fields are sown 
and the plants absorb the plant food. 

These are some ways to solve the problem 
of how to extract the greater part of the 
nutrients from sewage and at the same time 
purify the water so that it no longer destroys 
the watercourses. 

The burning of forest land and straw is a 
great waste of nitrogen. 

The waste water from our pulp mills should 
also be saved. 

Bark should not be allowed in the rivers, but 
should remain in the forests. 

Much effort must be expended at once to 
solve all these and allied problems, so that our 
most valuable heritage, our soil, shall not be 
laid waste. 


BIBLIOGRAPHY 


I. ANDERSSON, G., 1902: Hasseln i Sverige fordom och 
nu. Sveriges Geolog. Unders Ser. Ca 3. 

2. ARRHENIUS, O., 1953: Vissa ämnens fördelning i 
marken i Kopparbergs lan. Sveriges Geolog. Unders. 
Ser. C. Arsbok 44. 


3. — 1949: Artfördelningen i Stockholmstrakten, 
Svensk Botanisk Tidskrift. 
4. — 1940: Fosfathalten hos svenska torvslag. Sveriges 


Geol. Unders. Ser C. 
5. — 1934: Fosfathalten i skänska jordar. Sveriges 
Geolog. Unders. Ser. C. 


6. — 1952: The chemical denudation of the soil. 
Tellus 1952. 
7. — 1950: Markreaktionen hos sydsvenska jordar. 


Sveriges Geolog. Unders. Ser. C. Arsbok. 

8. ATLAS OVER SVERIGE, 1953: Utgiven av Svenska 
Sallskapet for Antropologi och Geografi. Red. Fil. 
dr M. Lundqvist. 

9. BAULIG, HENRI, 1910: Ecoulement 
denudation. Annales de Geographie XIX. 

10. BENNET, H. H., 1939: Soil Conservation, Mc Graw 
Hill Book Co, New York and London. 


fleuvial et 


11. BERGSTEN, F., 1950: Contribution to study of evap- 
oration in Sweden. Medd. Ser. D 3. Sveriges Met. 
Hydrol. Inst. 

12. Bortskaffelse og anvendelse af ensilagesaft. Statens For- 
sogsvirksomhed i Plantekultur, 486. Medd. 21 aug. 
1952. Stencilerat av Fiskeristyrelsens Sötvattens- 
laboratorium. 

13. BRAADLIE, O., 1930: Om elvevannets sammensetning 
i Tröndelag. K. Norsk Videnskabsselskabs Skrifter. 
14. Burp, JOHN S., 1918: Water extraction of soils as 
criteria of their crop producing power. Journ. Agric. 

Res. XII. 

15. CLARKE, F. W., 1924: The data of geochemistry. 
United States Geological Survey Bull, 770. 

16. CONVAY, E. J., 1942—43: The chemical evolution 
of the ocean. Proc. Roy. Inst. Acad. Dublin. 

17. Cooper, L. H. N., 1938: The nitrogen cycle in the 
sea. Journal of the Marine Biolog. Ass. Plymouth. 

18. CZAPEK, F., 1922: Biochemie der Pflanzen. Gustav 
Fischer, Jena. 

19. EBERMAYER, E., 1876: Die gesamte Lehre der Waldstreu. 
Berlin. 


Tellus VI (1954), 4 


| 
| 
| 


20. 


21. 


22. 


DE 


24. 


iS. 


26. 


AE 


28. 


29. 


30. 


32. 


32. 


33- 


CHEMICAL DENUDATION IN SWEDEN 


Erers, H., Saat, R. N.J., and Waarpen, M. v. D., 
1947: Chemical and physical investigations on dairy 
products. Elsevier. Amsterdam. 

ERIKSSON, JOEL, 1929: Den kemiska denudationen i 
Sverige. Medd. fr. Statens Meteorologisk-Hydrografiska 
Anstalt, 5. 

GRAPENGIESSER, S., 1934: Norrländska Vegetations- 
bilder. Svensk Botanisk Tidskrift. 

GRIPENBERG, S., 1937: The calcium content of baltic 
water. Journal du conseil, Charlottenlund. Vol. XI. 
HALDEN, B., 1931: Den lokala sedimentgränsen. 
Skogen 1931. 

—, 1928: Asken vid sin svenska nordgrins. Svenska 
Skogsvärdsföreningens Tidskr. 

Harvey, H. W., 1926: Nitrates in the sea. Marine 
Biol. Ass. U. K. Plymouth. 

HESSLAND, J., 1949: Calcareous fresh water sediments 
from Northern Bohuslin. Arkiv f. mineralogi o. 
geologi, I. 

HOFMAN-BANG, O., 1902—03. Studien über schwe- 
dische Fluss- und, Quellwässer. Bull of the Geolog. 
Inst. of the Univ. of Uppsala. VI. 

HOLMSTRÖM, S., och JOHANSSON, W., 1953: Jord- 
brukets ensilagesilor. Medd. fr. Jordbrukets Utrednings- 
institut 16. 

HULTÉN, E., 1950: Atlas över växternas utbredning i 
Norden. Esselte. 

HÂRDIN, AXEL, — : Samband mellan klimat, färg- 
kombination och medelförsämring. Korrosionsnämn- 
dens medd. to. 

Meddelanden frän Statens Meteorologisk-Hydrografiska 
Anstalt, 1914—51. I—14. De svenska vattendragens 
arealförhällanden. 

Medicinalstyrelsens föreskrifter angdende sättet för ut- 
fôrandet av fysikalisk-kemisk och bakteriologisk undersök- 


Tellus VI (1954), 4 
3—410071 


34 


35 


47. 


48 


341 


ning av vattenledningsvatten etc. Medd. fr. K. Medicinal- 
styrelsen 87—1950. 

. RANKAMA, K., and SAHAMA, Tu. G., 1950: Geo- 
chemistry. University of Chicago Press. 

. SCHILLINGER, S., 1927: Klär- und Reinigungsanlagen- 
für Münchener Abwässer in Zusammenhang mit 
dem Anbau des Wasserkraftunternehmens mittlere 
Isar A. G. Wasser. Verlag Chemie. 

. SCHNEIDEWIND, W., 1922: Die Ernährung der land- 
wirtschaftlichen Kulturpflanzen. Parey, Berlin. 

. SCHRÖDER, 1878: Fortschritte und pflanzenphysio- 
logische Untersuchungen. Tharandter Fortschr. Jahr- 
buch Suppl. 1. 

. SONDÉN, Kras, 1914: Anteckningar rörande svenska 
vattendrag I. II. Stockholm. 

. Statistik Ärsbok 1953. 

. Statistiska Centralbyräns Sockenstatistik 1951. 

. SVERDRUP, H. U., JOHNSON, MARTIN W., FLEMING, 
R. H., 1946: The Oceans. Prentice Hall. New York. 

. VALLIN, S., 1941: Betänkande angäende vattenför- 
oreningar II. Statens Offentliga Utredningar. 

. —, 1952: Förorening genom pressaften frän siloan- 
läggningar. Svensk Fiskeritidskr. 

. WENNSTRÖM, M. 1946: Konstgjord bevattning. 
Nordiska ingeniörsmötets Förh. Sthlm. 

. —, M., 1949: Biological purification of settled sew- 
age in shallow ponds. United Nations Scientific Con- 
ference, Lake Success. 

. VOGEL, J. H., 1896: Die Verwertung des städtischen 

Abfallstoffe. Arbeiten der deutschen landw. Gesellsch., 

ISL, IT 

— 1928: Abwasser. Ullmanns Encyclopädie der techn. 

Chemie. 

. WISTRAND, GUNNAR, 1934: Bidrag till kinnedomen 
om floran i Pite Lappmarks barrskogsrayon. Svensk 
Bot. Tidskr. 


The Carbon Dioxide Exchange Between the North Atlantic 
Ocean and the Atmosphere 


By A. NELSON DINGLE, The Ohio State University, Department of Physics and Astronomy* 


(Manuscript received June 28, 1954) 


Abstract 


An approximate distribution of carbon dioxide exchange rates between the atmosphere and 
the North Atlantic Ocean is obtained for each month of the year on the basis of oceano- 
graphic data from Böhnecke’s Atlas and the assumption that atmospheric carbon dioxide main- 
tains uniformly the volume percentage of .03. The results indicate the direction in which 
such an even volume percentage must be changed to approach equilibrium. 

Aside from other influences, the ocean-atmosphere exchange of carbon dioxide and the 
trajectories of tropospheric air masses combine to produce high carbon dioxide content in 
mT and cP air and low carbon dioxide content in mP air. Since water vapor is at a minimum 
in the continental polar anticyclones during the polar night, the carbon dioxide probably 
serves as an effective inhibitor of the outgoing infrared radiation. This inhibiting effect may well 
serve to limit the intensity of cooling during an ice age; it probably has much less significance 


during inter-glacial periods. 


Introduction 


Carbon dioxide as a quantitatively minor 
constituent of the atmosphere is seldom treated 
in detail in the textbooks. Meteorologists 
generally learn that carbon dioxide uniformly 
makes up about 0.03 per cent by volume of 
the atmosphere. Brief consideration of biolog- 
ical and industrial processes, however, leads 
one to observe that carbon dioxide must be 
quite nonuniformly distributed near the 
ground. The rates of influx or efux to or 
from the atmosphere vary with time in 
response to the diurnal cycle, the seasonal 
cycle, the work week, labor strikes, etc. These 


* Present affiliation: University of Michigan, Ann 
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effects produce such complexity in the carbon 
dioxide distribution near the earth’s surface 
that it can be accurately drawn only upon 
the basis of numerous observations at well- 
distributed stations over extended periods of 
time. Such observations are not available. 
There remains, however, the consideration 
of the carbon dioxide distribution in the air 
over the oceans. Since about seventy percent 
of the earth’s surface is occupied by oceans, 
conceivably the oceans may yield significant 
information about the broad scale distribution 
of carbon dioxide in the atmosphere. It is 
clear from simple gas solubility considerations 
that warm water will not retain as great 
concentrations of carbon dioxide as cold 
water will. It is to be expected a priori, there- 
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fore, that the warm parts of the oceans will 
tend to release carbon dioxide to the air, 
and the cold parts to absorb it. This suggests 
that on the basis of exchange rates one should 
find certain air mass types enriched in carbon 
dioxide relative to certain others, and that the 
synoptic distribution of carbon dioxide in the 
atmosphere may prove interesting. 

This synoptic distribution, in addition to its 
purely meteorological significance, is also 
pertinent to certain questions raised by at 
least three different writers during the last 
fifteen years. The first of these is G. S. Car- 
LENDER (1940) whose findings have been 
reiterated by E. GLUECKAUF (1951), i.e., that 
the artificial burning of fossil fuels since 1900 
could have produced an increase in atmospheric 
carbon dioxide equal to that which appears 
to have been observed up to 1935. The second 
is K. Buch (1939, 1942, 1948) who has attempt- 
ed to typify air masses of Northern Europe 
according to their carbon dioxide content. 
The third is M. Doze (1949) who has hypoth- 
esized that the observed excess of the mass 
18 isotope in atmospheric oxygen over that in 
the plants and waters of the earth with which 
the air should be in equilibrium, is maintained 
by the flow of carbon dioxide from the oceans 
to the air and subsequent release of O!8 from 
the carbon dioxide to atmospheric oxygen in 
exchange for O!6 under ultra-violet irradiation. 
In addition to these relatively recent works, the 
present study is certainly pertinent to the 
questions raised by T.C. Chamberlin and 


Svante Arrhenius almost sixty years ago, and 


perenially revived since, about the role of 


carbon dioxide in the radiational equilibrium 
of the earth. 


The Carbon Dioxide Exchange Rate 
The exchange of carbon dioxide between an 
aqueous solution and the atmosphere above it 
involves the following equilibria: 
_ (1) CO, (air) = CO, (aq) equilibrium 
L constant: & 
(2) H0 +CO, (aq) = H,CO, equilibrium 5 
k’ constant: K,= [m 
(3) HCO;=HCO;+H+ equilibrium 
Al constant: Ky 
CO, +H+ equilibrium 
constant: Ka 
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Clearly the process is quite complex even in 
its simplest form, i.e., using pure water as sol- 
vent. 

Generally speaking the reactions (1) and (3) 
are fast compared with (2). If one starts with 
this assumption, then at the water surface, 


[CO; (aq)] = x [CO; (air)], at equilibrium, 


and in the body of water [CO, (aq)] is deter- 
mined by reaction (2). If finally one assumes 
that H,CO; is uniformly concentrated through- 
out the volume of water communicating 
directly with the surface, then [CO, (aq)] (=) 
varies with x (depth measured from the sur- 
face) in a manner determined by the diffusion 
and hydration rates. | 


0) 
Ÿ a 
AX b 
A 
x Ÿÿ 
Figure 1. 


If N is the number of moles of CO, (aq), 
passing the level x within an area À, then 


ANa dN, _ ; 
ag Aka HCO) AEX 
d (dN ; 
oe (4) eco) 

Also IND D 2 assuming that the transfer 

dt dx 
is entirely due to diffusion. Thus 

2 
D at ee ECO, 


the solution of which is 


ik 


=X Zn 
C= (Co — Je D + € (1), 
where « is the concentration of carbon di- 
oxide in the water mass, below the region 
affected by the surface, i.e., at x, = co. Using 
equation (I), the curves of concentration, c, 
against depth, x, were computed for three 
typical situations (Figure 2). From these it is 
clear that the gradient of c with x has virtu- 


1°80 90 s00 NO 120 130 


Concentration,c,of COz in moles /dm? 


Figure 2. Carbon dioxide concentration, c, in ocean 

water as a function of depth, x, for three selected 

points, March. 
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ally vanished at a depth of 1 m in the case 

of pure diffusion. The effect of mixing is 

clearly to depress this level to greater depth. 
At the water-air interface, x = o: 


dN. = 

m VRD (ca — 6) (II) 
It is this equation which forms the basis for 
the calculation of the exchange rates. 

By definition, cy is the concentration of 
carbon dioxide (aqueous) at the surface in 
equilibrium with the overlying air, and hence 
it may be computed from the relation 


co = [CO; (aq)| =, [CO;(air)]. 


The atmospheric concentration of carbon 
dioxide over the oceans is not known in detail. 
It must vary widely with location from the 
tropical to the polar oceans. As a starting point, 
an assumed atmosphere may be used, and 
subsequent developments should reveal how 
such an atmosphere must be modified to 
approach an equilibrium condition. The ini- 
tial assumption is made, therefore, that carbon 
dioxide is uniformly present in the atmosphere 
in the amount of 0.03 % by volume. The 
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concentration, [CO, (air)], is thus directly 
proportional to the barometric pressure at the 
surface of the water. The solution coefficient 
a, is a function of the temperature and the 
salinity of the surface waters. The values 
derived by Harvey (1945) from the experi- 
mental data of Bour (1899) have been used 
in this work. 

The concentration c, of carbon dioxide in the 
ocean “mass” (that part of the ocean water in 
communication with the surface by diffusion 
and mixing) has been computed by means 
of the oceanographic technique as described 
by Harvey (1945). The result of Harvey’s 
derivation is: 


_ carbonate alkalinity -[H*]- «, 
; Ke 
K, S (: ar A) * 4,0% 


In oceanographic practice, the various charac- 
teristics of sea water such as “carbonate alka- 
linity”, “chlorinity” and “total boron” are all 
related empirically to the “salinity” which is 
observed. The above equation thus may be 
written: 


C 


a = 

ely (12 09 ee) [ a 
Te 

In this equation: 

CES = 220 (& 15 to 20) 


Ky is the apparent dissociation constant of 
boric acid in sea water (= 107? to 2.5: 107 4 


[H*] is the hydrogen ion concentration 
(~ 0.5: 1078 to 0.9: 10” 8) 
Kj is the first apparent dissociation constant of 
H,CO; (* 0.6: 10-7 to 1.1+ 1077) 
Ky is the second apparent dissociation constant of 
H,CO, (& 0.5-10-® to 1.3-10-9) 


ao is the activity of water in sea water 
(~ 0.98 to 1.00 hence taken to equal 1) 


%/&, is the ratio. of the solubility of carbon 
dioxide in sea water to that in pure water 
(& 0.8 to 0.9) 
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Figure 3. 


Each of these quantities has been related ex- 
perimentally to temperature and salinity 
(Buch, 1939). For computational purposes one 
therefore constructs curves and reads the 
appropriate values therefrom. 

Representative monthly mean temperature 
and salinity distributions in the surface waters 
of the North Atlantic have been tabulated by 
months and published in BOuNecke’s Atlas 
(1936). Monthly mean distributions of sea 
level pressure were procured from the United 
States Weather Bureau. 

After computing & and c as outlined above, 
the rates of exchange may be computed 
pointwise by means of equation (I). 
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Distributions of annul average carbon dioxide exchange rates, 1o—* mol m—? sec-1. 


— = VkD (cree) (II) 


The diffusion coefficient D is evaluated by 
means of a curve derived from data given in 
the INTERNATIONAL CRITICAL TABLES (1928). 
The hydration rate k depends upon tempera- 
ture and pH, but its variation with pH is 
relatively small. The pH of the surface waters 
varies between the limits 8.0 to 8.3 (Harvey, 
1945) hence for this purpose it is assumed 
constant at 8.15. k is then evaluated for each 
point by means of the thermodynamic relation 


a 
log, k=log0 ko — m te - +) 
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Figure 4. Distribution of monthly mean carbon dioxide exchange rates, 10-4 mol m-? sec-1. 


in which m was found from FAURHOLT’S (1924) 
determinations to be equal to 4220. 


The Results 


The present work is necessarily somewhat 
crude, and hence to be regarded as a prelimi- 
nary study. The network of points for which 
the computation has been made is limited to a 
minimum for the purpose, i.e., 48 points 
were chosen well distributed over the area. 
Data for some of these were missing for some 


months, so that the map for January is based 
upon only 38 points, and those for February 
and December upon 40 points each. 

The results are presented graphically in 
Figure 4. The map used is the equal-area 
projection of Mollweide, so that areal integra- 
tion might be accomplished by use of a planim- 
eter. The exchange rates are given in units 
of 1074 mole m”? sec-1. Application of an 
areal summation to the annual mean exchange 
chart (Figure 3) shows that a .03 % carbon 
dioxide atmosphere would have to gain a net 
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Figure 4. Distribution of monthly mean carbon dioxide exchange rates, 1074 mol m? secu’. 


amount of carbon dioxide in order to arrive 
at a dynamic equilibrium with the North 
Atlantic Ocean. A qualitative inspection of 
temperature and salinity relations in the other 
oceans of the world suggests that the world- 
wide ocean-atmosphere carbon dioxide equilib- 
rium requires still greater increase of atmos- 
pheric carbon dioxide than that over the North 
Atlantic, thus it is evident that the actual aver- 
age carbon dioxide content of the atmosphere 
must exceed 0.03 % somewhat. Itis noteworthy 
that no great adjustment is needed, however. 
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The charts of Figures 3 and 4 show numer- 
ous interesting features. Parallelism with the 
isotherm pattern is as anticipated. The reduc- 
ed rate of emission of carbon dioxide in the 
area affected by fresh water drainage from 
the Amazon is quite marked, as is the similar 
reduction near the Gold Coast of Africa. This 
latter effect is produced mainly by cold water 
from the Benguela Current and from up- 
welling. Other points may be made, however 
it is not intended that detailed interpretations 
should be constructed upon the basis of these 
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October 


charts which are not themselves based upon 
detailed data, but that a general picture may 
be presented. 


General Conclusions 


Reference was made earlier to the work of 
BucH (1939, 1942 and 1948) relating atmos- 
pheric carbon dioxide content to recent air 
mass history. Buch reports a minimum of 
carbon dioxide in air found near Spitzbergen 
which appeared to have a cool maritime 
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Figure 4. Distribution of monthly mean carbon dioxide exchange rates, 10-4 mol m? sec =. 


history, and relatively large amounts in con- 
tinental air from the East and in maritime 
tropical air. It is clear from Figures 3 and 4 
that maritime polar air in order to approach 
equilibrium with the surface waters must 
become relatively depleted in carbon dioxide 
as compared to average conditions, and that 
the reverse is true for maritime tropical an 
If one envisions a meridional interchange of air 
following the convectional model of tbe 
general circulation (HADLEY, 1735—6; RossBy, 
1941) it is reasonable to hypothesize that the 
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air subsiding in the polar anticyclones must be 
relatively heavily laden with carbon dioxide. 
With respect to atmospheric carbon dioxide 
therefore, the tropical oceans are the principal 
source regions and they serve to maintain 
higher than average levels of carbon dioxide 
concentration in both continental polar and 
maritime tropical air. It is obvious that all 
gradations of carbon dioxide content must be 
present depending upon the respective tra- 
jectories, path lengths, speeds of motion, and 
low level stability of particular air masses. 
Thus in order to classify air masses by means 
of their carbon dioxide content one must 
trace their histories rather carefully. For this 
reason the attempt of CALLENDAR (1940) to 
perform such a classification is only moder- 
tely successful. 

Implicit in these observations are the follow- 
ing tests for the general hypothesis set forward, 
1. In relatively stable conditions a vertical 

gradient of carbon dioxide content (volume 
percentage) should be observed directed 
upward over the tropical oceans and 
downward over the polar oceans. 

2. A nearly uniform content should be found 
upon sounding the subsiding air of the 
polar anticyclones. 

3. Precipitation formed in maritime tropical 
air should contain higher concentrations of 
dissolved carbon dioxide than that formed 
in maritime polar and transitional air 
masses. 

A carbon dioxide sounding device has been 

developed by Wicciams et al. (1953). Prelim- 

inary results obtained over New Mexico in 
summer show a rather sharp drop in carbon 
dioxide content upon ascending through the 
tropopause. Here evidently (over New Mexico 
in summer) the tropospheric stirring is sufhi- 
ciently strong that no appreciable inhomogene- 
ity occurs below the tropopause. That a definite 
decrease occurs upon passing through the 
tropopause suggests that mixing between the 
troposphere and the stratosphere is not as 
rapid as most meteorologists have believed. 

Unfortunately, the Williams instrument is not 

capable of giving quantitative information at 

this stage of development. 

CALLENDAR’S (1940) hypothesis that an 
increase of atmospheric carbon dioxide has 
caused a rise of temperature over much of 
the earth, is peculiarly challenging. In the 
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first place, cause and effect in meteorology are 
seldom that clearly discernible; and secondly 
there are numerous other possibilities, reason- 
ably plausible, which he appears not to have 
considered. It is clear, for example, that a 
warming of the ocean, from whatever cause, 
would require a shift toward higher carbon 
dioxide content of the air in order to preserve 
the dynamic equilibrium of the ocean-atmos- 
phere exchange. 

The possibility that long term shifts of the 
general circulation such as those implicit in 
secular trends of temperature and other weather 
elements might have altered the mean trajec- 
tories of air (and sea water) needs also to be 
considered. If the frequency of maritime trop- 
ical air invasions of Europe were to increase 
somewhat at the same time as the frequency 
of maritime polar air invasions decreased, one 
would expect an increase of carbon dioxide 
in the air of this region. If the trajectories 
were so modified that the maritime tropical 
air maintained a longer contact with warm, 
saline ocean water, and the maritime polar air 
was exposed more briefly to the cool and less 
saline waters, again an increase of the observed 
carbon dioxide should appear. 

In light of the knowledge that the North 
Atlantic Ocean has become warmer during the 
period to which Callendar refers, it appears 
that the carbon dioxide concentrations in the 
air of Western Europe would be directly 
affected because of its trajectory. It is not 
known why the North Atlantic has become 
warmer. A sustained circulation shift might 
well account for this and contribute additional 
effects such as those discussed above. Thus 
the view that the observed changes are account- 
able to secular changes of the general cir- 
culation as the basic phenomena in this chain 
of events seems most attractive. 

As regards the role of carbon dioxide in the 
radiation balance of the earth, two points are 
worthy of mention. The first is that prelimi- 
nary results (Howarp, 1954) of high disper- 
sion infrared absorption studies of carbon 
dioxide and water vapor reveal that the ab- 
sorption of these two gases is essentially 
additive. This would tend to refute HUMPHREYS’ 
(1940) well-known conclusion that carbon 
dioxide isn’t important because water vapor 
does a better job of absorbing the same wave 
lengths. The second point is that, according 
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to the model envisioned above, the polar 
anticyclones should be laden .with carbon 
dioxide. This would place a relatively large 
amount of carbon dioxide into a region 
where water vapor is seriously depleted, and 
it would seem that this is an especially oppor- 
tune distribution for controlling the flux of 
infrared radiation from the earth during the 
polar night. It appears conceivable there- 
fore, that atmospheric carbon dioxide may 
serve the very important role of limiting the 
outward flux of heat during periods of maxi- 
mum glaciation, since in such periods very 
little water vapor is available over the ice caps. 
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On the Structure of Some Cumulonimbus Clouds which 
Penetrated the High Tropical Troposphere' 


By JOANNE STARR MALKUS and CLAUDE RONNE 
Woods Hole Oceanographic Institution 


(Manuscript received 25 May 1954) 


Abstract 


Some extremely large oceanic trade-wind cumulonimbus clouds extending upwards of 40,000 
ft into a region of strong winds and intense vertical shear have been studied by means of time- 
lapse photography. A simultaneous still picture of the clouds taken a known distance and di- 
rection away from the motion pictures permits triangulation upon the clouds and calculation 
of the vertical and horizontal motions of several of the individual towers. By means of a nearby 
radiosonde observation, it is established that the air forming the strongest of these towers must 
have risen from near cloud base with little or no dilution of its buoyancy by mixing with the 
clear-air surroundings. 

The model of Malkus and Scorer (1955) concerning the erosion of cumulus towers is reviewed 
and tested upon these towers with satisfactory results. A minimum horizontal cloud dimension 
is apparently required for the production of undiluted towers. The horizontal motion of the 
clouds relative to the air is also estimated from the model and tested by the observations and the 
net upward transport of latent heat in water vapor is calculated approximately. 


I. Introduction fairly weak wind shear. Large, precipitating 
cumulonimbus have proved somewhat for- 
midable to investigate by means of aircraft, 
especially when accurate vertical draft measure- 
ments are desired. Fortunately, the new 
theoretical model of cloud towers is well 
suited to testing by measurements from time- 
lapse motion pictures, with a minimum of 
additional information, including a nearby 
clear-air temperature sounding. The present 
paper is a study, based on time-lapse photog- 
raphy, of some oceanic trade-wind cumu- 

1 Contribution No. 716 from the Woods Hole lonimbus clouds whose towers penetrated 
Oceanographic Institution. heights exceeding 40,000 ft and wind shears 
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Recently it has been found possible (MArkus 
and SCORER, 1955) to formulate quantitative 
laws for the rise and decay of individual 
cumulus towers and to describe their time 
“ dependence in terms of relatively few parame- 
ters. So far, however, most of the observe- 
tions used to check the theories on cumulus 
have been made upon rather small, non- 
precipitating trade-wind clouds, which have 
been ascending through air masses exhibiting 
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Fig. 1. Time cross section for San Juan, Puerto Rico, from 1200 GCT, March 31, 1953 to 0600 GCT April 2, 1953. 

The ordinate is height in feet. Winds are plotted in knots, a short barb indicating 5 knots, a long barb 10 knots, 

and a solid triangle so knots. The solid lines are isotachs of wind speed, drawn every Io knots, with J indicating 

maxima. The dashed line indicates the base of the westerlies, and the dotted line shows the top of the moist layer. 

The heavy solid line marked PT indicates the polar trough which passed San Juan moving eastward on 

March 31, and the solid line marked EW denotes a weak easterly wave which passed San Juan moving 
westward early on April 1. 


as high as 30 mps/km. In addition to testing 
the theory in this rather extreme case, some 
important differences between large and small 
trade cumulus are brought out which may be 
interpreted by means of the model. 

These data were taken on April 1, 1953, 
during a joint observational program in the 
Caribbean carried out by the meteorological 
group from the Woods Hole Oceanographic 
Institution and the Department of Meteorol- 
ogy ofthe Imperial College, London. Through- 
out this program, two main observation 
points were used, the first being the small, 
flat West Indian island of Anegada (18°so’ N, 
64°20° W; dimensions two by ten miles) 
and the second a PBY-6A aircraft, equipped 
as a meteorological tool, which was flown 
in the near vicinity. Supplementary infor- 
mation, such as radiosonde and radiowind 
data, was obtained from the U.S. Weather 
Bureau at San Juan, Puerto Rico (120 miles 


west-southwest from Anegada). The cumulo- 
nimbus clouds on this occasion were photo- 
graphed by a time-lapse camera mounted 
on the beach at Anegada, and sufficient addi- 
tional information was available to permit 
numerous calculations from the films. After 
a brief description of the prevailing meteoro- 
logical conditions and the equipment and 
techniques used in the measurements, the ma- 
jor portion of this paper will be concerned 
with these calculations and their interpretation. 

The synoptic situation on April ı was 
dominated by the passage on the preceding 
day of a strong polar trough. In the Caribbean 
area, the low-level easterlies were overlaid 
by westerlies beginning at 15,000 ft elevation. 
These westerlies were concentrated in a 
shallow, high-tropospheric layer with maxi- 
mum speeds of roo mph at an elevation of 
40,000 ft. The sequence of events is illustrated 
in tbe time cross-section for San Juan shown 
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Fig. 2. 


Still photograph from the PBY showing the clouds studied. The arrows indicate the points upon 


which triangulation was performed. The aircraft was being flown at 6,000 ft toward 293° True and its location 
is given in Figure 5. 


in Figure 1. At the time of the present observa- 
tions, 1220—1345 LST (1620—1745 GCT) 
the region was dominated by the northwesterly 
flow to the rear of the trough. By this time 
the overcast associated with the disturbance 
had broken sufficiently for several huge cumu- 
lonimbus clouds to be observed and photo- 
graphed from Anegada. Fortunately the same 
clouds were simultaneously photographed 
with a still camera from the PBY, which was a 
known distance and direction from Anegada. 
Triangulation therefore permitted calculations 
of the distances, heights, and sizes of the 
clouds and their alteration with time from the 
motion picture film. The resulting vertical 
and horizontal motions could then be com- 
pared with the high-level wind and radio- 
sonde observations (which were taken at 
1500 GCT by the Weather Bureau at San 
Juan) and with the relations prescribed by 
the model. 


II. Methods of Observation and Calculation 


The time-lapse motion picture camera on 
Anegada was a Bolex 16 mm, equipped with 
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an extreme wide-angle (58.3° horizontal, 
44.3° vertical) lens. Exposures -were made 
automatically at the rate of one frame every 
3.6 seconds or 1,000 frames per hour. The 
camera was mounted on a tripod and aimed 
toward 22° from true North during the entire 
run. 

The still camera in the PBY was a conven- 
tional Speed Graphic (50.0° horizontal, 42.0° 
vertical) which was mounted on gimbals so 
that it always remained level and was pointed 
at right angles to the fuselage. At the time 
the still pictures of the clouds were taken 
(1325 LST) the aircraft was flying at 6,000 ft 
elevation toward 293° True and was making 
a speed of approximately 130 mph over 
the ground. Since it was observed over Ane- 
gada on this heading at 1314 or 11 minutes 
earlier, a base leg for triangulation purposes 
of 24 miles is obtained. This distance was 
checked by additional knowledge of the 
time when the aircraft reached a location due 
north of San Juan. The base leg derived thus 
by interpolation came out 23 miles or within 
5 %. Correction of the base leg for wind 
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Selected prints from the time-lapse motion picture run used in the calculations. The individual towers 


studied are denoted by the arrows with symbols above. The bubble followed on Tower Y (one of the strongest) 
is shown by the horizontal line; Tower W (finishing its cycle at the beginning of the run and therefore not 
included in the six towers studied) is indicated by the square; Tower C (weak) is shown by the arrow alone; 
Tower Q (weak) is denoted by the inverted triangle; the streamer on Cloud II is marked by the erect triangle (used 
only as a check in the calculations); and Tower X (strong) is indicated by an X. The camera was mounted on the 
beach at Anegada and pointed toward 22° True throughout the run. The print number is given beneath each frame 
and consecutive prints were I.2 minutes apart (every 20 frames of the motion picture film, which was exposed at 
the rate of one frame every 3.6 seconds). 


drift of the aircraft was made, but was negli- 
gibly small. 

Figure 2 shows the Speed Graphic (PBY) 
photograph used in the triangulation. The 
time-lapse pictures were printed at intervals 
of 1.2 minutes, and the prints numbered from 
the beginning to the end of the run (1—s8). 
Selected prints from this series are reproduced 
in Figures 3 and 4. Print 44 was picked as the 
frame coincident with Figure 2. Triangula- 
tion was first performed for the three points 
marked with inverted arrows (on Figure 2 
and Print 44) on the far cloud, hereafter to 
be called Cloud I. The angle of each point 
from the optical axis of the camera was 
found by using the proportional distance of 


the picked point from the center to the edge 
of the negative. The half-width of the nega- 
tive and focal length of the lens were accu- 
rately known and these gave the tangent 
of the total half-angle covered by the camera. 

Results of this and a similar triangulation 
on one point on the near cloud (called Cloud 
II) are shown in Figure 5. By laying out two 
angles for each point at the ends of the base 
leg on a scaled map, it was determined that 
the three points in Cloud I lay not only in a 
straight line but almost exactly parallel to 
the image plane of the motion picture camera 
(i.e. in the plane of the photographs) and a 
perpendicular distance from it of 78 miles. 
The point on Cloud II was a perpendicular 


Tellus VI (1954), 4 


§ 
| 


STRUCTURE, OF CUMULONIMBUS CLOUDS 355 


Fig. 4. Additional selected prints from the motion picture camera on Anegada. On Print 40, Tower P (weakest 

tower studied) is shown by the arrow alone; Tower X is indicated by an X, and the streamer on Cloud II by an 

erect triangle. Print 44 was the one used in triangulation (compare to Fig. 2). The arrows alone indicate the points 

triangulated upon (see results in Fig. 5). The arrow alone on Prints 52 and 58 (the last one in the run) shows 
Tower R. The dissipation of the entire cloud mass may be noted on the last print. 


distance of 66 miles from the image plane. 
After discussing the possible errors in tri- 
angulation, we may carry out calculations 
from this point forward using the motion 
picture frames alone. 

In the case of Cloud I, for which all the 
significant calculations were performed, the 
image location never exceeded an angle of 
+ 10° from the optical axis, hence errors in 
distances and separations due to distortion by 
the lens may be totally neglected compared to 
errors in triangulation. An error of 2.2 miles 
in the base line (9 %) would give an 8 % 
error in the distance to Cloud I, and it is 
doubtful that the error in the base line was 
this great. The largest source of error lies in 
the determination of the angles. A 2° error 
in any one of the angles would produce a 
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7 mile or 9 % error in distance to Cloud I. 
All the angles save one may be shown reliable 
within this accuracy. The greatest uncertainty 
lies in the angle between the optical axis of 
the Speed Graphic and the heading of the 
PBY, which was adjusted to 90° by sighting 
along the trailing edge of the wing (fortu- 
nately constructed at exactly 90° to the air- 
craft axis). That the errors in angle were very 
small, however, is suggested by several inde- 
pendent checks upon the triangulation. The 
first is that all three points on Cloud I lie in a 
straight line along the plane 290°—110°, as 
an analysis using the wind profile shows must 
be the case. The second is that triangulation 
gives horizontal separations on the map 
between these points departing by much 
less than 10 % from the separations measured 


to 
[m 
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CAMERA TO CLOUDI 
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15 MI 30% 
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Fig. 5. Scaled map showing results of the triangulation 
using Figure 2 and Print 44, Figure 4. The camera on 
Anegada was located at À, pointing in the direction of 
the dashed arrow (toward 22°). The PBY was located at B, 
heading in the direction indicated by the solid arrow. 
The three points indicated by arrows on Cloud I (Fig. 2 
and Fig. 4) lie along the line marked “Cloud I’ which 
lies very closely in the plane 290°—110° and 78 miles 
perpendicular distance from A. The single point tri- 
angulated upon in Cloud II lies 66 miles perpendicularly 
from the plane 290°—110° through A. 


on the photographs, using the derived distance 
to the camera of 78 miles. The final check is 
that triangulation on Cloud II gives its heights 
and velocities as extremely close to those of 
comparable portions of Cloud I. 

The vertical angular distance (on Print 44) 
between the horizon and the farthest left 
point on Cloud I is determinable from the 
ratio of this distance to the vertical half- 
width of the negative and knowledge of the 
focal length of the lens. Using the horizontal 
distance of 78 miles and the tangent of this 
angle, the vertical elevation of this point above 
the horizon was found to be 27,000 ft or 8.24 
km. This length could then be used as a basic 
distance scale for the calculations to follow. 
In order to obtain the elevation of this point 
above sea level, however, a correction had to 
be introduced due to the curvature of the 
earth. For a distance of 78 miles, this proved 
to be 4,000 ft or 1.2 km. In presenting actual 
altitudes of all points on Cloud I, this correction 
has been added. 

In order to apply the distance scale obtained 
from Print 44 to other parts of the cloud and 
to motions of portions of the cloud towers 
throughout the entire 70-minute run, it was 
first necessary to ascertain whether the cloud 
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had any component of motion into or out 
of the plane of the photograph and whether 
it might develop any slope normal to that 
plane. Inspection of the San Juan radio-wind 
observation (see Figures 1 and 9) shows that 
the wind directions throughout the height 
range of the calculations (30,000—47,000 ft) 
departed rarely, and never more than by 10°, 
from the plane of the photographs. Although 
the wind departed significantly from this 
direction at lower levels, it may be shown by 
consideration of even the minimum drag 
forces ordinarily acting on cumulus that the 
component of cloud motion normal to the 
photographs should have become vanishingly 
small by 9.5 km. Thus the cloud, which was 
shown by triangulation to have no significant 
slope perpendicular to the plane of the photo- 
graphs on Print 44, should not develop any 
such slope with time, nor should it move any 
appreciable distance toward or away from the 
camera during the 70 minutes considered. 
The distance scale for Print 44 was therefore 
used throughout without correction, with 
extremely consistent results. 


II. Basic Features of the Major Cloud 


The main dimensions of Cloud I are now 
readily obtainable from Print 44. At 31,000 ft 
(9.5 km), the width of the left-hand tower 
(containing the left-most point triangulated) 
is 20,500 ft or 6.3 km. The greatest elevation 
above sea level attained by any portion of the 
cloud is reached by the streamer extending 
out to the right, which achieves a maximum 
height of about 15.3 km or 46,600 ft and a 
horizontal extent of 45 miles. The striking | 
feature of this cloud’s changes with time on 
Prints 1—58 (70 minutes) is its production | 
of six successive towers, each of which is: 
seen to rise vertically, spread laterally into a. 
long streamer and then dissipate gradually. . 
The towers were found to vary widely in: 
their dimensions, maximum ascent speeds, | 
and heights penetrated. The two most vigorous ; 
towers are marked by a horizontal line and! 
an X on the prints shown in Figures 3 and 4.. 
A typical one of the weaker towers is denoted | 
by the inverted triangle. It is important to 
note that each of these towers underwent| 
an initial period of rapid rise which petered | 
out into a period of much slower rise some- 
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Fig. 6. Tephigram showing San Juan radiosonde obser- 
vation, 1500 GCT, April 1, 1953, reduced from the 
original data! to obtain more closely spaced points than 
that of the teletype transmission. The sounding is given 
by the solid line. The dashed line is the moist adiabat 
intersecting the sounding at cloud base (950 mb) used to 
obtain the buoyancy for the strongest bubbles, and 
the dotted line is the moist adiabat used in obtaining the 
buoyancy of some of the weaker bubbles. 


where between about 11.0 km and 12.5 km 
(or about 3—5 minutes after the tower first 
emerged from the main cloud mass). The 
period of slow rise averaged about 10 minutes, 
while another 10—1$ minutes elapsed, in 
general, during the dissipation of the streamer. 
In the case of the two most vigorous towers, 
close inspection of the films reveals that the 
transition from rapid to very gradual ascent 
took place near 12.6 and 12.2 km, respectively, 
where the appearance of the towers also 
changed markedly and sharp outlines became 
lost. This transition could be in no way related 
to the transition between supercooled water 
and ice crystals which must have occurred at 
no higher than 10.9 km, where even a maxi- 
mally buoyant cloud element reached — 40.0°C. 

The San Juan radiosonde observation for 
1500 GCT is reproduced in Figure 6. The 
sounding has a most striking feature. It shows 
that any cloud “element which is finally to 
reach the vicinity of 12.5 km with vanishing 
(and not negative) buoyancy, must have 
risen moist adiabatically and undiluted all the 


1 The writers are grateful to Dr. Herbert Riehl and 
Mr. José Colon of the University of Chicago for ob- 
taining and reducing these data and also those for the 
radiowind observation in Figure 9. 
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way from cloud base at 950 mb. Any mixing 
or dilution by the clear-air surroundings on 
the way up, even a small fraction of that 
encountered by STOMMEL (1947), Markus 
(1952, 1954) and others in smaller trade 
cumulus would have resulted in prohibitively 
large negative buoyancies at this level. The 
next section is intended to shed light on this 
point. 


IV. The Erosion and Life Cycle of the Cloud 


Towers 


a. The model to be used. — Markus and 
SCORER (1955) showed by theoretical and 
observational argument that the tops of 
actively rising cumulus towers should obey 
the following differential equation: 


a + Kw?=gB (1) 


where w is the ascent velocity of the tower; 
Kis a drag coefficient to be specified presently; 
and gB is the buoyancy acceleration which is 
; ES be ‘ : 
given by g (257) minus a correction for 
v 

suspended hydrometeors, where T, is the 
virtual absolute temperature within the tower 
and T, is the virtual absolute temperature of 
the surroundings. The Kw? term represents 
the drag exerted by the surrounding air upon 
the tower and has been called form or profile 
drag by aerodynamicists. When it is small 
or zero, equation (1) returns to the familiar 
equation for the non-interacting buoyant 
parcel. When the vertical acceleration dw/dt, 
hereafter called w, becomes zero, the buoyancy 
and drag are balanced and define a limiting 
velocity, wr, which is given by 


w =(gB/K)! (2) 


The drag coeflicient, K, may be shown by 
dimensional argument to be inversely pro- 
portional to a linear dimension of the tower 
element; thus, the time dependence of some 
linear dimension must be specified in order to 
proceed further with the differential equation. 

In obtaining an expression for K, Markus 
and Scorer (loc. cit.) have attempted to 
construct a model of cloud towers which 
permits physical interpretation of some of the 
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Fig. 7. Schematic diagram of a single buoyant element, 
showing the bubble cap with rounded upper surface 
followed by a turbulent wake. The radius of curva- 
ture of the cap is denoted by R and the cap diameter 
by D. The angular half-aperture is @m which is by defi- 
nition equal to sin-+ D/2R. For the cloudy bubbles 
studied by Markus and Scorer (1955), D and R were no 
the average nearly equal, as shown here. This also proves 
to be the case for the large elements in the present study. 
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processes at work. They have envisaged a 
cumulus tower to be composed of one or 
more of the bubble-like elements shown 
schematically in Figure 7. Each element is 
hypothesized to consist of an undilute cap 
(which rises moist adiabatically) and a trailing 
turbulent wake. The upper surface of the cap 
is maintained quasi-spherical and its lower 
surface is flattened and irregular. The turbulent 
wake consists of moist air shed from the cap 
and drier air incorporated by entrainment 
from the surroundings. Using this model, 
equation (1) may be derived fairly rigorously 
and it may further be shown that 


K = (3) 


where R is the radius of curvature of the round- 
ed cap. The only crucial assumptions are 
that the actual ascent velocity, w, of the cap 
departs only little from the limiting velocity, 
wr, defined by equation (2) and that the cap 
retains its shape (or ratio of diameter D to 
radius R) during its lifetime. Both of these 
assumptions are justifiable theoretically if the 
rate of decrease of R with t (supposed due to 
shedding of the outside layer of the undilute 
cap into the wake) is sufficiently slow and 
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observationally they have been supported by 
measurements upon towers of small and me- 
dium-sized cumulus clouds. Alternatively, 
the theory so far can be derived by analogy 
with work on air bubbles in liquid performed 
by Davis and TAYLOR (1950). 

In the case of cumulus towers it was found 
that the relation 


R = EgB(-t) (4) 


could be made plausible physically and that 
it fit the available observations well. Malkus 
and Scorer called equation (4) an erosion law, 
where E is a constant of proportionality 
called the erosion constant, to be discussed 
presently, and f is chosen equal to zero when 
R becomes equal to zero and is negative 
during the elements active life. Equation (4) 
was only derived accurately for the case in 
which gB, the buoyancy of the element, was 
constant: that is for a wet-adiabatic environ- 
ment. Combining (4), (3), and (r), and solving 
for a constant proportional to E, we have 


SE CS (- 1)? gB(1 — w/gB)* (s) 


w 


If the hypothesized laws are valid, 2 E4 or 


G should be a constant for all cloud tower 
elements which are sufficiently isolated to be 
treated as individual bubbles. Equation (5) is 
ideally suited to testing by calculations from 
time-lapse motion pictures, if the buoyancy 
gB is known or can be estimated with the 
aid of a clear-air sounding and if it does not 
vary too rapidly with time. It was shown 
that (5) could be used with varying values 
of buoyancy provided that gB did not change 
by more than about 1/3 in 100 sec. An in- 
crease of buoyancy at this rate might overes- 
timate G by as much as 20 %. It should be 
emphasized that from the level of their 
formation upward, the cloud bubbles are 
hypothesized to be rising along a moist 
adiabat. The mixing with the environment is 
achieved by erosion and the «core of the cap 
remains undiluted until its final disappearance. 
The cumulus bubbles already studied were 
found to originate at different levels, some 
over land forming at or near the ground, 
others at cloud base, and still others up to 
the highest level within the cloud where 
the external lapse rate exceeded moist adia- 
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Fig. 8. Height of top of bubble in Tower X in km as a function of time in seconds. The time, —f, is chosen as 
zero when the bubble levels off and begins to appear in fuzzy outline on the film. Several points commonly 
occur on the same horizontal because readings could be made only to the nearest 0.05 cm on the microfilm 
viewer. Readings were made both in a forward and backward direction on the film and were found to coin- 
cide to this degree of accuracy. Velocities were obtained by reading first differences of the heights and were 
checked by average values read from the prints. Accelerations were obtained by taking second differences in 
heights (see inset graph) and were checked by taking first differences in velocities. The calculations made from 
these curves are presented in Table I. 


batic. No oceanic trade cumulus bubbles have 
been found which originated below cloud base. 
The ones formed at higher levels were gener- 
ally composed of smaller bubbles plus wake 
material, and thus were more “dilute”” or 
had less buoyancy than those forming lower 
down. 

The tests of the model performed by Malkus 
and Scorer upon small and medium-sized 
cumulus bubbles showed that G was constant 
to an excellent degree of accuracy, both during 
the life of a single cloud tower and from tower 
to tower, regardless of location or external 
conditions. The mean value of G came out 
slightly greater than 0.2 sec}, giving an 
E of about so+1o seconds. The highest pene- 
trating bubble studied was a trade cumulus 
tower reaching 16,000 ft which did not, 
‚however, become glaciated. 

b. The life cycles of the individual towers 
observed—The cloud towers in the present 
study not only reached great heights, but 
penetrated tremendously high wind shears 
and definitely contained ice crystals throughout 
a significant portion of their life histories. It 
should be of considerable importance to 
determine whether equation (5) and the 
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constant value of E obtained for the bubbles 
studied previously give here also a con- 
sistent relation between velocity, acceleration, 
buoyancy, and time. If so, it will be of interest 
to see how the present observations may be 
interpreted and clarified by means of this 
model. Most of the present towers have the 
advantage of being well isolated from inter- 
ference with one another due to great wind 
shear and thus form an especially valuable 
test case. 

To make this test, the original motion 
picture film was projected on a microfilm 
viewer, from which the heights of the bubble 
tops were read off at every frame (3.6 sec). 
Although the vertical resolution was still 
smaller than desirable, there being commonly 
only about 1 cm height difference between the 
emergence of a bubble and its death as an 
entity, the results of drawing a smooth curve 
through the points were extremely satisfactory. 
The measurements and method of calculation 
of velocities and accelerations for one of the 
strongest towers (called X on Figures 3 and 4) 
are shown in Figure 8. A dual criterion was 
used for picking the time of complete erosion, 
when (- 1) was to be set equal to zero. This 
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time was chosen when the levelling off of 
the ascent was indicated on the graph, coin- 
cident with the loss of the rounded upper 
surface showing itself on the film. This 
transition is readily identified since the hard 
upper edges rapidly become fuzzy and diffuse 
when the spherical shape is lost. In applying 
and testing equation (s) on Tower X, there- 
fore, the quantities w, w, and -£ were taken 
from the curves shown in Figure 8. It was 
then desired to see whether reasonable and 
observationally supportable values of gB, the 
bubble’s buoyancy, could give rise to values 
of G which both remained nearly constant 
throughout the bubble’s life and averaged 
about 0.2 sec”? as established for the smaller 
cloud towers. When the values of the buoyancy 
are chosen with the aid of the sounding in 
Figure 6, under assumptions to be described 
presently, the resulting calculations of G are 
obtained which are shown in Table I. Under 


these conditions the average G, or G, is 0.195 
sec? and the percentage spread in G, called 
A, is only 25 %. 

The values of buoyancy in the column 
marked gB in Table I may be obtained for 
this bubble from the sounding in Figure 6 
under several assumptions, the most crucial 
being that the San Juan radiosonde observation 


Table I. Tower X, April 1, 1953. Height range 
studied 10.42—12.24 km 


oe -t | (-?#)2 w wv eB G 
num! sec | sect mps ora) a sec? 
ber x sec? sec? 
I 288: | 17:02, 11700 2:21. 14703 MO 
2 2702| 70.5.0 1220 2 221017230 2027 
3 252 || 15:9 | 10.572 2207207 No 
4 225 Aa en 9.8 | — 2.2] 11.5 | 0.20 
5 Paley | a7) 2 || 20 ont || au 
6 198 14.I See || 2 9.0 0.18 
7 TSO) 013,5 7.2 | "2" 01 8.8 | 0.19 
8 162 12.8 6.4 | —3.0] 8.6 | 0.20 
9 144 | 12.0 5.8 | —3.0| 8.4 | 0.20 
10 126 112 GG) |) 2 S:2H| NO 10 
II 108 10.4 5.3 | — 2.2] 8.0 | 0.18 
12 90 | 95 4:7 | ==3 0078) | 0.19 
13 72 8.4 AVA ——- GeO RO OT 
14 54 74 3.90 0 3:2 87.62 On eS 
15 36 6.0 3.0 | — 4.8 7.6 0.19 
16 18 4.2 2.0 | — 6.2 Wao || 0722 
* Bubble diameter measured Ga Ontos 
as 1.33 km. A = 25% 
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is applicable to the clear-air environment of 
the clouds. This is readily established by pro- 
jecting the trajectory of the air surrounding 
the clouds at the appropriate heights (9.4—14 
km) backwards in time to the time of the 
radiosonde. It is thereby shown that the air 
in the cloud locality at the middle of the 
photographic run was 100 miles due north 
of San Juan at the time of the radiosonde. 
If we consider, for example, the weak stable 
layer occurring between 10 and 11 km (which 
coincides with the layer of high wind shear) 
we see that it would have to have a north- 
south slope of 1/150, or that common for a 
middle-latitude front, in order to occur at a 
10 % different elevation 100 miles to the 
north. In the trade-wind region, a north- 
south slope this great is out of the question. 
It is thus concluded that the observations on 
the sounding curve apply to the air in the 
cloud region with less than a +10 % shift 
in height. Local time changes within the air 
mass are certainly smaller still. The second 
assumption in obtaining gB is that the ascent 
curve of the bubble forming Tower X is very 
closely given by the dashed (right-hand) 
moist adiabat on Figure 6, which originates 
from the environment sounding at cloud- 
base level. 


That this curve delineates the 
maximum possible buoyancy for this tower is 
supported by other work on oceanic trade 
cumuli (Markus, 1954) showing that the 
cloud buoyancies at cloud-base level rarely 
exceed that contributed by water vapor alone. 
It will be noted that this adiabat recrosses 
the sounding curve just above 12.5 km, 
leaving the bubble some buoyancy at its 
demise at 12.24 km. Thus this adiabat, ac- 
cording to the bubble model, also represents 
near minimum buoyancy for the tower since 
isolated bubbles cannot continue to rise 
without positive buoyancy. In proceeding 
from the virtual temperature differences 
between the dashed adiabat and the sounding 
curve to the values of gB presented in Table I, 
a negative correction to the buoyancy due to 
the presence of suspended hydrometeors has 
been introduced. This correction allows for 
about 4.5 gm/m? of water or ice particles at 
10.42 km, diminishing to 2 gm/m at the 
bubble’s demise at 12.24 km. This correction 
amounts to nearly 30 % of the total buoyancy 
and is both uncertain and arbitrary: it is simply 
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delimited by that reduction in buoyancy 
necessary to proceed from the dashed adiabat 
to the values of gB in the table which produce 
values of G consistent with those of the 
previous study. That these figures are not 
unreasonable, however, may be seen as follows: 
if all hydrometeors were being retained in the 
bubbles from cloud base upward, the water or 
ice content at 10.5 km would be about 18 
gm/m?. These clouds were, however, observed 
to be precipitating heavily from their bases. 
Thus a hydrometeor content from 10—25 % 
of this figure is plausible, especially when it 
falls within the range of observed water 
contents in large cumuli. The latent heat of 
freezing has been ignored in the present 
computation, the assumption being that most 
of the freezing took place below 10.42 km 
where the bubbles began to be tracked. 


The other vigorous tower produced by 
this cloud (called Tower Y and denoted by 
the horizontal line in Figure 3) was also stu- 
died. It should be pointed out that Tower Y 
was composed of several bubbles, which 
accounts for its great width. Each of these 
was about the size of the single bubble com- 
posing Tower X and in Figure 3a four are 
plainly visible. The bubble studied is the 
right-hand one which ascended first. Its maxi- 
mum rate of rise was the same as that of 
Tower X, namely 11.7 mps. The buoyancy 
for Y was also obtained from the dashed 
adiabat in Figure 6, with the same magnitude 
correction for the weight of suspended 


hydrometeors. G came out 0.18 sec? with a 
percentage spread, A, of 45 %. It is believed 
that the presence of the neighboring bubbles 
disturbed the constancy of G, especially 
during the later portions of its life cycle. 
This was due to overly large ascent speeds 
acquired by this bubble’s being pushed up by 


successors, resulting in too low values of G. 
Next, one of the weaker towers was studied 


in the same manner and the results are pre- 


sented in Table II. This is the bubble called Q, 
which was denoted by the inverted triangles 
in Figure 3. The significant feature of Table 
II is that a considerably lower buoyancy had 
to be chosen for this bubble to give the same 


value of G as the others. The buoyancy for 
this bubble was taken from the moist adiabat 


shown by the dotted line (left-hand adiabat) 
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Table II. Tower Q, April 1, 1953. Height range 
studied 9.72—10.94 km 


Obser- ; 

vation| =? | (-#)# eS = gB ¢ 

bee sec | sect | mps or = sec”? 
re EOE | 4 8.9 o 13.0 | 0.21 
2 198 14.0 8.6 | —0.9| 12. 0.21 
3 180 | 13.4 8.4 | —ı.2| 12.0 | 0.20 
4 7524 072.3 8:1 | — 1-6) 11.1, | 0.18 
5 144 | 12.0 7.2 | —ı.9| 10.5 | 0.19 
6 126 baz 6.4 | —2.2 9.8 0.19 
7 108 | 10.4 5.9 | — 2.8 A, | 0.116) 
8 90 9.5 5.5 | —3.1| 8.8 | 0.18 
9 2 8.5 Bale |b 3.7 0058:44 6.0217, 
Io 54 TEA 4.1 | —4.3| 8.0 | 0.18 
II 36 6.0 3.0m) —— 5 7100| 70.179 
12 18 4.3 1.8 | —6.2| 6.0 | 0.20 
* Diameter of bubble measured G = 0.191 
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on Figure 6, with a correction for weight of 
hydrometeors of 3 gm/m? at 9.72 km dimin- 
ishing to 2 gm/m? at its maximum height of 
10.94 km. The resulting values of gB given 
in Table II average 0.55 times the values 
used for X and Y in the height range 10.42 — 
10.94 km. This latter bubble (Q) possessed 
lower ascent rates and perished sooner than 
the preceding two described. It can now be 
seen that it must have been a considerably 
more dilute bubble, its level of origin being 
perhaps about 5.9 km or 19,300 ft rather 
than cloud base. It presumably formed from 
smaller bubbles and dilute wake material near 
this level, while the other two formed near 
cloud base. Its concomitantly smaller required 
temperature excess (about 1.3° C compared 
to about 2.6°C for X and Y) over the sur- 
roundings at 5.9 km is consistent with this 
point of view. 

It will be seen by inspection of Figures 3 
and 4 that Tower Q rose on the extreme left 
side of Cloud I (upwind and upshear). In 
fact at the time it formed, this part of the 
cloud was almost entirely separate from the 
larger cloud mass to the right. The average 
width of the smaller cloud mass was only 
about 4 km, while Tower Y rose from a 
cloud mass at least 9 km in width. By the 
time Tower X rose from the left-hand por- 
tion of cloud, this section had grown until 
it, too, had attained an average width of 
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about 9 km. These data might thus suggest 
that, under given conditions, a minimum 
horizontal cloud dimension is required to 
permit emergent buoyant elements to rise 
all the way from cloud base without dilution. 
Since the undiluted elements will, in general, 
compose the strongest and most penetrating 
towers, it may well be that for cumulus to 
attain a sufficient vertical thickness to give 
rise to precipitation, a comparable horizontal 
size must first be attained. 

To investigate these points further, the 
other weak towers put up by Cloud I were 
studied. There were three such towers de- 
tectable on the film. One (P) was too feeble 
and short-lived to be treated quantitatively. 
One of the two remaining, called C, appeared 
on Print 12 just as the bubble studied in 
Tower Y was closing its life cycle. Although 
it appeared from the same large cloud mass as 
did Y, its buoyancy was the same as that of 
the diluted tower, Q. It had an active lifetime 
of 185 sec, a maximum updraft speed of 
7.2 mps, and a height range of 10.4—11.2 km. 

The final tower was the very last one 
appearing on the film sequence, called R. 
In order to have the same value of G as the 
other towers, it had to be slightly dilute and 
to follow a moist adiabat about half-way 
between the dashed and the dotted curve on 
Figure 6. It covered a height range from 
10.6—11.85 km, exhibited a maximum up- 
draft speed of 7.8 mps, and a life period of 
275 sec, just short of those of the strongest 
towers. At the time of its emergence, the entire 
cloud mass of Cloud I was apparently beginning 
to decay, and although the bubble came from 
the main portion of the cloud, the horizontal 
dimension of that body was then only a 
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little over 7 km. Thus it is seen that although 
the strength of the emerging towers and the 
dimension of the producing cloud mass ia 
general appear to be related, another factor or 
factors must also be operative, since the largest 
cloud mass produced weak as well as strong 
towers. A suggestion concerning one of these 
factors is made in later paragraphs. : 

The major features of the six towers put 
up by Cloud I are summarized in Table II. 
The towers are arranged in order of their 
appearance, the times shown being minutes 
elapsed since the beginning of the run. In 
the table, the measured diameter of each 
bubble (measured at the starred times in 
Tables I and II) is compared to the radius of 
curvature calculated from equation (4) using 
E=50 sec. The earlier study by Malkus and 
Scorer showed, as Figure 7 suggests, that D as 
measured and R came out very nearly equal, 
indicating aperture angles, 20,,, for cloud 
bubbles of about 60°. The good agreement in 
the present cases gives further support to the 
applicability of the model. The basic assump- 
tion that the actual ascent velocity w only 
slightly exceeds the limiting velocity wz is 
also supported. When the calculated R is 
used to obtain w, from the combined equations 
(2) and (3), it is found that on the average w 
exceeded wz, by 6 % at the starred times. 

The column marked Rx gives the approxi- 
mate dimension of each bubble as it first 
emerged from the cloud top, calculated by 
extrapolation from equation (4). It will be 
noted that the strong, undilute towers, X 
and Y, had maximum radii in excess of 2 km, 
while the weaker ones just exceeded 1 km. 
There is a corresponding variation in the 
maximum ascent velocities at the time of 
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emergence, given in the column marked wnax. 

Two of the four weaker towers emerged 
from smaller cloud masses while the other 
two emerged from the largest (— 9 km) 
ones. In the case of the latter, namely C and 
P, one noteworthy feature is held in common. 
In each case, the weak tower followed its 
predecessor’s emergence after an unusually 
short interval, 7.2 minutes in the case of C, 
and 6 minutes in the case of P. The strongest 
towers, Y and X, were separated by more 
than a half hour, while X and R were separated 
by nearly 20 minutes. If, as suggested by 
writers on the bubble model, the larger 
cloudy bubbles are formed by aggregation of 
smaller ones at low levels within the cloud, 
this fact becomes somewhat clarified. A certain 
time interval is apparently required for a 
sufficient number of small bubbles to accumu- 
late in the low cloud levels in order to form a 
large vigorous one which can rise to emerge 
undiluted from cloud top. 

It may be seen that each bubble once it 
emerged from protection by the main cloud 
body penetrated a vertical distance very 
nearly equal to its maximum radius at emer- 
gence. Furthermore, those bubbles which 
did not follow their predecessors too closely 
rose through a height within cloud (from 
their origin level to the top of the cloud body) 
very comparable to the width of the pro- 
tecting cloud body. Similar conclusions were 
also deduced by Malkus and Scorer for the 
smaller cumulus towers. It appears that the 
presence of ice crystals has no noticeable 
effects upon the bubbles during their active 
lives, although it may be responsible for 
slower dissipation of their wakes and remnants. 

The last column on the right shows that 
the very vigorous undilute bubbles had active 
lives of about 300 sec or five minutes (X 
could have been followed this long on the 
film) and successively weaker bubbles had 
successively shorter lives. Bubble P probably 
existed independently for 100 sec or less. 
As noted earlier, the remnants of the towers 
went on rising at 1/,—1 mps for about 10 
minutes after the round tops were lost, 
achieving heights greater than 13 km. Then 
dissipation set in and occupied another 10—1$ 
minutes without much further ascent taking 
place. 


Since the present discussion has in part 
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concerned undiluted cloud elements, some 
rising all the way from cloud base moist 
adiabatically, it must be emphasized that any 
resemblance to the hypotheses and results of 
the parcel method is strictly superficial. 
The parcel method considers cloud elements 
which in no way interact with the surrounding 
air and which would continue to accelerate 
upward wherever the parcel density is less 
than that of the environment. In the present 
case, a rough calculation shows that parcels 
becoming saturated at cloud base would 
emerge from cloud top at to km with upward 
velocities well in excess of 60 mps and still 
increasing! Since the bubble, on the other 
hand, is hypothesized to be travelling at very 
nearly its limiting velocity at all times (at 
least after emergence) its interaction with the 
surroundings is held to be absolutely vital. 
In fact, the buoyancy force is just over- 
balanced by form drag (produced by flow 
around the bubble cap) and is not producing a 
resultant upward acceleration. On the con- 
trary, the bubble’s upward velocity is in 
general decreasing with time due to erosion, 
despite the fact that the buoyancy is positive. 
To summarize: the parcel method prescribes 
no relations between the velocity and radius 
of a convective element and an erroneous 
relation between velocity, buoyancy, and 
time. The present model suggests quantitative 
relations between these parameters, which 
are tested here. 

c. The horizontal motions of the towers and 
the effects of wind shear upon their life cycles. — 
Figure 9 shows the San Juan radiowind obser- 
vation for 1500 GCT. The small triangles, 
squares, and X’s denote the horizontal speeds 
of several of the cloud towers at the heights 
plotted. These represent speeds in the plane 
of the photographs which, it will be remem- 
bered, fortunately coincided with the plane 
of the ambient wind. Each such symbol rep- 
resents an average horizontal speed from one 
print to the next (1.2-minute interval). It 
proved too difficult to get 3.6-second hori- 
zontal speeds from the microfilm viewer as 
was done in the case of the vertical speeds. 
Nevertheless several deductions may be drawn 
from Figure 9. It is first of all encouraging to 
note that the cloud towers began to shear 
horizontally at the level at which the intense 
shear in the ambient wind set in (~ 9.7 km). 
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Fig. 9. The wind profile plotted from the San Juan 
radiowind observation for 1500 GCT, April I, 1953. 
These points, marked by circles, were obtained upon 
reduction of the original data for this study. The curve 
on the left is wind speed in meters per second; that on 
the right is direction in degrees from true Morth. The 
layer of strong shear has been found on other soundingsin 
the area, and on previous and subsequent soundings at 
San Juan. The other symbols, such as triangles, bars, 
squares, and X’s are horizontal speeds of the various 
towers, as denoted on the figure. 


This fact supports the validity of the height 
calculations from the films and the applica- 
bility of the sounding to the cloud region. 
Second, it is noted that above heights of 
about 12.5 km, where even the most vigorous 
bubbles ended their active life cycles (made a 
transition to much weaker ascent rates), the 
horizontal speeds of the towers draw close to 
(but never exceed) the outside windspeed. 
More detailed consideration of the hori- 
zontal motions of the towers of Cloud I are 
possible, however. It was shown by Malkus 
and Scorer that in the case of a tower ascending 


through a shearing wind field equation (s) 


should be replaced by 


(2) gB(r — w/gB)? 
== : = 
2 w (wo? + u?)} 


(6) 


where u is the horizontal velocity of the bubble 
relative to the air. This equation is readily 
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derived by setting up a differential equation 
analogous to (1) for the horizontal motion of 
the tower (where drag forces only are oper- 
ative) and writing an expression for the 
total vector drag force which depends now 
on (u2+w?). When the vertical component of 
this total drag force is substituted into (1), 
K is eliminated and (6) is obtained. It has 
been shown that in most shearing wind fields 
the difference between the results of (5) and 
(6) would be only about 2%, but when u is 
of the same size as w, use of (5) gives about a 
20% overestimation of G. Figure 9 indicates 
that in the present case u is indeed comparable 
to w, since in the region of strong shear the 
horizontal tower velocities fall short of the 
wind velocity by 9—15 mps. However, it 
may also be shown that a rate of buoyancy 
decrease of the magnitude experienced by 
the present bubbles at these heights should 
lead to about a 15—20 % underestimation 
of G. Since these errors are compensatory, 
and since the present observations are probably 
not reliable to better than 15—20 % anyhow, 
no attempt will be made here to introduce 
such refinements. 

The data in Figure 9 also provide a means 
for a direct test of the horizontal equation of 
motion for cloudy bubbles. Analogously to 
(1), this equation may be shown to be 


er (u? + w?) 


(7) 


Di =. 
Vu? + uw? 


where 1 is the component of the bubble’s 
horizontal velocity along the axis determined 
by the external wind, ug; u is the relative 
velocity between bubble and surroundings, 
namely u=u,- ug; and the dot indicates time 
differentiation. The constant K was shown to 
be the same as K in equation (3), namely 
9/4R. Subtracting u; from both sides of (7), 
we have 


$ u 
v= — WU = 


— K (uv? + w?) (8) 
Vu? + w? 

since %=wU’=wdUldze and = 11 — Up, 
where U’ or dU/dz is the external wind shear. 

If the bubble has been rising through 
constant shear for- several hundred meters 
and if it is still fairly near the beginning of its 
life cycle so that K is changing only slowly, 
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u is very small, as is w, and (8) may be solved V. Conclusion 


approximately to give 


u „RU 


ee 


w OVu? + w? 


(9) 


after substitution of (3) for K. 


The region of strong shear on Figure 9 
corresponds to the earlier phases of the cycles 
of the towers studied. Hence (9) may be 
tested. The figure shows that Tower X 
possessed a # of about —15 mps in this region 
and that for Tower Q, u averaged about 
—12 mps. Assuming that these bubbles had 
been rising through the shearing layer long 
enough for ü to be very small, and taking the 
observed U’ of 3 x 10-2 sec, the radius of 
curvature R for Tower X comes out 1.8 km 
and R for Tower Q comes out 1.4 km from 
equation (9), compared to the calculated values 
of 2.1 and 1.3 km respectively which appeared 
in Table III. These independent results are 
consistent with the predictions of the model 
and agree within observational error, although 
it must be conceded that the measurements 
and the approximations made are such that 
this high degree of agreement is somewhat 
surprising. 

It is not clear from the observations whether, 
as the bubble continues its life cycle, the 
horizontal drag coefficient increases in the 
same manner as its vertical counterpart, 
which the model declares that it should, 
although this is at least roughly suggested by 
the fact that u=u, — ug seems to be decreasing 
in magnitude before evaporation sets in and 
makes it impossible to measure further. 
Changing wind shear coupled with the diffi- 
culty of tracing a single point in its horizontal 
travel precludes an accurate quantitative test 
of (8) solved for G (after substitution of K= 
9/4R and R=EgB(-1) such as was done for 
the vertical equation (1). 


In a shearing wind field a rising bubble 
should have a horizontal as well as a vertical 
component to its trailing turbulent wake. 
Earlier observations by Malkus (1949) suggest 
the presence of such a wake region in the 
horizontal, since several hundred feet of “dying 
cloud” were commonly found on the down- 
shear edge of cumuli. 
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The clouds studied in this paper are of in- 
terest for two main reasons: first, they are 
the largest clouds in conditions of strongest 
shear to which the techniques described here 
have been applied. Previous clouds investi- 
gated by the Woods Hole group have been 
small or middle-sized trade cumulus and fair 
weather cumulus of middle latitudes, none of 
whose tops reached the freezing level. Some 
similarities and some contrasts between these 
clouds and the smaller ones studied earlier 
are striking. The outstanding similarity is the 
applicability of the same erosion constant, E, 
to the towers. The predominant contrast is 
the relatively less important role played by 
mixing with the clear-air environment in 
the bigger clouds. This is not to say that 
interaction between cloud elements and their 
surroundings is any less important, since it 
was demonstrated here that the strongest 
bubbles studied emerged from cloud top with 
velocities less than 1/6 that of a non-interacting 
parcel. In fact, it has been indicated that a 
minimum cloud dimension is necessary to 
shield the innermost core from dilution. In 
the case of these clouds, a mass roughly 9 km 
across (in the plane of the wind and the wind 
shear) could produce at intervals of about 
20 minutes, undiluted bubbles 2 km in di- 
ameter which emerged from cloud top rising 
at 11 mps. The relation between radius, 
buoyancy, velocity, and lifetime of these 
elements agreed very satisfactorily with the 
relations prescribed by the model of eroding 
cloud towers. 


The greatest importance meteorologically 
of the clouds studied is, however, their role in 
the over-all picture. It must be supposed that 
clouds of this type and scale, although rarely 
observed from the ground due to prevailing 
overcasts under disturbed conditions, are the 
rule and not the exception when an extended 
polar trough reaches down to trade-wind 
latitudes. 

It is well established (RiEHt, 1950) that the 
strong and high-penetrating convective activity 
which occurs in these troughs is the means 
by which the upper troposphere receives the 
major fraction of its latent energy supply. 
A sample calculation will illustrate the magni- 
tudes involved. Assume that at the 25,000 ft 
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level, a single saturated cloud tower is active 
at a given moment in an area 30 km on a side 
(roughly correct for April 1, 1953). If the 
updraft speed is 10 mps and the draft is 2 km 
on a side, continuity will be met by sub- 
sidence at 4 cm/sec over the remaining area. 
Since even the more dilute towers at this 
elevation have about 3.0 gm/kg water vapor 
content and the ambient air from the sounding 
less than 0.5 gm/kg, the net upward transport 
of water vapor is probably not less than 
6x 107 gm/sec or 3.6x 101° cal/sec in latent 
heat. Distributed over the area involved, this 
is 400 cal/cm? per day. This figure should 
be compared to a several months’ aver- 
age evaporation rate of 250 cal/cm? per 
day calculated for a similar latitude in the 
Pacific tradewind region by RıEHL, YEH, 
Markus, and La Seur (1951). Clouds such 
as these thus play a major role in main- 
tenance of the energy budget of the 
trades and on that basis alone merit further 


study. 
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Charge generation in Thunderstorms by collision of ice crystals 


with graupel, falling through a vertical electric field 
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Abstract 


When an ice-crystal impinges on the lower surface of a particle of hail or soft hail part of 
the charge induced by the electric field is transmitted to the ice-crystal. The relaxation time for 
the exchange of charge varies from a few milliseconds to about 10-f seconds and is sufficiently 
short to enable this gain of charge. In a positive field the soft hail will therefore get a negative 
and the ice-crystal a positive charge. In an electric field E the final charge of a sphere with a 
radius Rise, = —2.12 ER?. At temperatures between — 18° and — 26° the number of ice crystals 
is sufficient to allow a time constant of about $ minutes for the charging process. The charge 
per droplet and per g of water derived from this process are in good agreement with measure- 
ments of Ross Gunn. The growth of precipitate and the formation of strong electric fields can 


be brought into correspondence with the timing and development of thunderstorms. 


Introduction 


During the years 1937—1941 G. C. SIMPSON 
and collaborators [1] published the results of 
experimental investigations on the distribution 
of electricity in thunderclouds and the asso- 
ciated electric fields. This initiated a new era 
in our understanding of the origin of elec- 


tricity in thunderclouds. The observations 


made it clear that electrical charge coincides 
with certain regions of temperature. Negative 
charge is centered in regions of temperatures 
below freezing point. This region of negative 
charge is marked by the 0° to- 20° isotherms. 
The positive charge is about 2—3 kms higher. 
It is roughly centered at the - 10° to -45° 
isotherms. A third and smaller region of 
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prevailing positive charge frequently, but not 
always, appears during thunderstorms between 
the +10° and o° isotherms. The free charge 
in this region is considerably smaller averaging 
about 1/5 of those in the other two regions. 
The measurements of Simpson have subse- 
quently been extended by others. According 
to E. I. WorKMANN and E. S. ReyNozps the 
centre of negative charge is between the 0° 
and -25° isotherms, but mostly at -—15° 
Pays ee 

The generation of electric space charge and 
electric fields is always associated with preci- 
pitation phenomena about which the following 
is known: It appears that a first precipitate can 
be formed in either the liquid or the solid state. 
Larger drops originate inside the cloud in a 
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temperature environment above the 0° iso- 
therm. According to Byers [2, p. 63] the exam- 
ination of 66 radar-echos shows that larger 
drops appear—in the average—660 m above 
the 0° isotherm (1e. at an altitude of about 
5,600 m), but in many cases only 300 m above 
the 0° isotherm, i.e. at temperatures of —4.0° 
to -2.5°. According to WORKMANN and 
ReynoLps [2, p. 147] the first radar-echo is 
associated with the -21° to - 24° isotherms, 
where it appears to be reflected by ice crystals. 
When the radar cloud reaches the — 28° to 
— 29° isotherm a noticeable downward move- 
ment of the top of the cloud begins simul- 
taneously with precipitation resulting in the 
separation of electric charge [2, p. 139]. The 
radar-cloud of a specific thunderstorm-cell 
will sink with approximately the same veloc- 
ity with which it grew (2.5—20 m/s). This 
marks the beginning of electrical discharges 
lasting 25—40 minutes. According to Byers, 
however, the first lightning-discharge is ob- 
served when the radar cloud passes the — 20° 
or an even lower isotherm, 10—1$ minutes 
after the appearance of the first radar-echo. 
The investigation of 35 cases has shown that 
radar-clouds were growing at a mean rate of 
5.4 m/s, 90 cm/s less than that of the visible 
clouds. At a height of 6,000 m a down-drift 
velocity of 90cm/s is reached by drops of 
about 100 u radius or by soft rime of about 
600 u radius which are responsible for radar- 
echos from the top of the cloud. 16 minutes 
after the first radar-echo has occurred and a 
mean cloud height of 10,500 m (-36° iso- 
therm) has been reached the average down- 
drift-velocity of the top of the cloud is 
3.6 m/s. The downward current of air is 
propagated from the lower part of the cloud 
to the top in connection with the formation of 
precipitation. It therefore commences a short 
time after the first radar-echo has appeared. 
10—1$ minutes later, when the first lightning 
discharge occurs the -20°— -28° isotherm 
is reached. The most frequently found width 
of that area in which the air flows downward 
is 1,500 m (with a maximum of 7,000 m). 
At an altitude of 7,800 m the mean down- 
drift-velocity is 6.7 m/s with a maximum of 
24 m/s (mean of 60 measurements) [3 p. 42]. 
Simpsons measurements by means of 
soundings have given evidence on the field 
strength within the clouds: In wide regions 
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it seldom exceeds some 100 V/cm. Measure- 
ments made by using aeroplanes [3 p. 80] 
yielded values up to a little above s00 V/cm. 
Ross Gunn obtained a maximum gradient 
of 1,600 V/cm at an altitude of 3,900 m, 
shortly before the plane was hit by lightning. 
This value refers to the unperturbed field [2 p. 
199]. (The field was of course disturbed by 
the aeroplane, the actual measured value being 
3,400 V/cm [2 p. 199].) 

The positive and the negative charge are 
by no means completely separated. According 
to Ross Gunn [2 p. 200] in the region from 
about +14.7° to -ı0° both positive and 
negative drops are mixed. But there are also 
cells extending over 2 km containing uniform 
negative charge. In these mixed regions the 
effective density of space charge is consider- 
ably smaller than the space charge density 
resulting from the charge on single drops. 
Positive space charge prevails in the region 
between 6,000 m and 8,500 m where precipi- 
tation with positive charge has been exclusively 
identified. 

The electric processes connected with the 
formation of ice crystals have been investigated 
by W. FINDEISEN, E. WALL, and R. ROSSMANN 
[4], H. Lueper [s], E. J. WorkMANN and 
S. E. REYNOLDS [2]. They have been critically 
examined by B. M. Mason [6a]. During the 
formation of a rime layer a negative charge is 
acquired by the ice. This process is capable of 
generating the required quantity of charge 
(B. I. Mason [6 b]). It must be supposed 
that several electrical processes will occur si- 
multaneously in a thundercloud. We shall 
now investigate the process taking place when 
rime and ice crystals come into contact and 
are again separated in an electric field. In the 
event of such a contact between rime and an 
ice crystal the crystal will slip off the ice grain 
after touching it for a short moment. At the 
same time the ice crystal removes part of the 
charge induced on the “bottom” of the ice 
grain by the electric field. In a positive field 
the small crystal will become positively 
charged, consequently leaving the bigger ice 
grain negatively charged. ELSTER and GEITEL 
discussed this mechanism for drops in 1913 [7]. 
In this manner a rain drop was to convey 
part of its induced charge to a cloud droplet. 
An air cushion between the drop and the 
droplet was to prevent coalescence of the two 


Tellus VI (1954), 4 


CHARGE GENERATION IN THUNDERSTORMS 


particles, and permit the positively charged 
droplet to recoil from the drop. The experi- 
mental investigation of T. E. W. SCHUMANN 
[8] proves, however, that the process of colli- 
sion between two liquid drops is different: 
Small droplets having a radius of less than 
180 y will coalesce with the raindrop. Larger 
droplets will travel through the raindrop and 
then disintegrate taking with them the charge 
induced on the top. The process then ends 
with a charge of opposite sign to that assumed 
in Geitel’s and Elster’s theory. The raindrop 
will receive a positive charge. Thus the rain- 
drop weakens the existing positive field. 
Consequently this mechanism does not in- 
tensify the field and will itself slowly die out. 
We assume, therefore, that the crystals 
slip off the soft hail in the form of small 
plates of varying thickness, as dendrites or 
as small frozen droplets. This will not—or not 
always—be the case with soft rime, at least 
not while the rime is still developing a conical 
habit through little supercooled droplets loose- 
ly depositing on a star-shaped crystal (spec. 
weight about o.r). Such forms of soft rime 
are found as well as soft hail, and later on we 
shall see that this is particularly significant. 


Relaxation Time of the Charge-Exchange 
Process of Ice 

If charge is to flow from an ice grain to a 
small crystal of ice the latter must have a 
certain conductivity. The relaxation time of 
the exchange process must be shorter than the 
time of contact. This time is given by: 


€ 
A9" 10 1-# 


(1) 


tr= 


Now neither the dielectric constant ¢ nor the 
conductivity x both of which vary with 
temperature, are defined values. Water always 
contains some impurities. This holds even for 
normally distilled water of a conductivity of 
‚about 1 x 1076 ( cm)! because of absorption 
of ammonia and carbondioxyde from the 
air. The dissolved substances are not “tightly” 
frozen. The ice is permeated by a network of 
very thin channels. In respect to its dielectric 
constant and its conductivity it behaves as an 
inhomogeneous dielectric. Both properties 
depend strongly on the frequency. Table I 


shows some data taken from an investigation 
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of C. P. Smrrm and C. P. Hırcacock [9]. 
The figures refer to ice obtained from water 
having a conductivity of 2 x 10% (Orne 
at the applied frequencies given in the left- 
hand column. 
The relaxation time obtained by this method 
ranges from 10°? to 10° s within the ranges of 
temperature and frequency given above. 
The relaxation time may also be obtained 
from other experimental data. It is assumed 
that at the instant of contact between ice 
grain and ice crystal the electric phenomena 
in the dielectric are similar to those at high 
frequencies. The dielectric constant will be 
determined by the relation &=3.40 +0.04 
(Pro) (Tt aineedes-s Cent.) "This empirical 
relation was found by H. Brommets [10] to 
hold for temperatures down to —18° (°C), 
and the corresponding lowest value of 3 for 
the dielectric constant will be used later on. 
The D.C. conductivity of ice like that of other 
ionic conductors, e.g. quartz, calcite, rock salt 
etc. is given by 


log x= 4-7 (2) 
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Fig. 1. Relation between conductivity of ice and tem- 

perature. 


Fig. 2. Relation between relaxation time of exchange of 
charge and temperature. 


with the two constants A and B and the abso- 
lute temperature T,. The constants A and B 
which determine the magnitude of the con- 
ductivity depend on the number of foreign 
ions absorbed by the water. Carefully distilled 
water has a conductivity of 1+2x10*% 
(2 cm)! while the conductivity of rain water 
is greater by about one order of magnitude. 
For temperatures ranging from 0° to - 19° C 
the values of the conductivity of ice are taken 
from JOHNSTONE [11], and for — 78° that value 
obtained by Banper at Lorg’s laboratory, 
namely 1074 (Q cm)-! (fig. 1). These data 
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refer to ice prepared from distilled water and 
are minimum values. Then the conductivity 
is given by 
| 3:78.x.10° 
Og %= 5.35 T, 


(22) 


This leads to a lower limiting value, because 
there are always very small quantities of 
sodium, magnesium, and chlorine dissolved 
in the ice crystal, particularly if they originated 
directly from supercooled droplets. Although 
no measurements are known, H. KÖHLEr’s 
investigation [13] carried on for a number of 
years suggest the following: melted soft rime 
shows a chlorine content of (0.067 --56.3) x 


10 6g per cm? of water with a certain frequency: 


distribution having a maximum at 3.5 x 10® 
g per cm?. This soft rime originates from a 
rather loose freezing (specific gravity 0.01 to 
0.1) of small and very small droplets with a 
radius of 8.8 u as the most frequent value. 
This value leads to 1.85 x 10714 g of salt per 
droplet, this being the most frequent conden- 
sation nucleus, and is responsible for a greater 
conductivity as that given by (2 a). The values 
published by H. Köhler may be assumed to 
hold also at 7,000 m altitude, because when air 
is driven up vertically from the ground it will 
carry with it condensation nuclei from ground 
level. The measurements of H. NORINDER and 
R. SIksNA [14] give further support to the 
assumption that the conductivity of ice 
crystals—i.e. soft rime, snow, soft hail, and 
hailstone—is greater than that given by (2a) 
(provided the crystals do not originate from 
sublimation on water-insoluble nuclei). These 
measurements, obtained with compact snow, 
are shown in fig. 1. 


Table I 


frequency 
kcycles 


dielectric 
constant € 


conductivity x 
ro IC) 


relaxation time fR 


(5) 


temperature °C 


BO 
. TOP 
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An approximate value for the conductivity 
of rime and hail will henceforth be assumed 
having an upper limit one order of magnitude 


higher than that given by (2a), i.e. 


3.78x 10° 


= (2b) 


log x=6.35- 


In this way the relaxation-time assumes the 
values shown in fig. 2 (shaded area): At 
temperatures between —10°C and -20°C 
the relaxation-time varies between 1.6 x 10-3 
seconds (upper value, ice-crystals originating 
from sublimation) and 3 x 10°? seconds (lower 
value, rime). At temperatures below ~40° C 
which are most frequently encountered at 
altitudes above 12 km the relaxation-time for 
ice-crystals is about 1.8 x 107? seconds. There- 
fore ice-crystals and ice-grains may be justifi- 
able regarded as being “conductive” within 
that range of temperature in which the elec- 
tricity of thunderstorms probably originates. 


Electrification of ice by contact in an 
electric field 


The following arrangement was investigated: 
A particle of rime is polarized by a vertical 
electric field through which it falls. Let it 
overtake small crystals of a shape shown, for 
instance, in the micrographs of H. KOHLER 
[13 a] or V. SCHAEFER [15]. For a first approxi- 
mation we consider the particle of rime (fig. 3) 
replaced by a conductive sphere of greater radius 
and the crystal by one of smaller radius. 
The greater sphere of radius R, be charged in 
proportion to the gradient thus 


a =pxER; (3) 


where B denotes the “specific” charge. An 
infinitesimal change of charge is 


de, = dB x ER? (3 a) 
The charge transmitted by contact from the 
greater sphere to the smaller one of radius Ry 
is proportional to the radial field strength at 
the point of contact: 


&=YyXER, (4) 


where y denotes a numerical factor which, 
2 


BAR Ta 7 
when the ratio —° ~ — is equal to — =1.65 
Rı 20/106 6 


Tellus VI (1954), 4 


371 


mae Se 
mm 


Fig. 3. Soft rime and soft hail. 


and may be determined for other ratios from 
fig. 4. 

At the surface of the drop polarized in the 
electric field E the radial field strength E, 
at an angle # is 


E,=E (3 cos ÿ +B) (5) 


The maximum field strength is obtained for 
d»=0,i.e.for the centre of the circle representing 
the projection of the sphere on a horizontal 
plane. The probability of a collision between a 
small sphere and the bigger one is regarded 
as constant all over the circle representing the 
projection of the bigger sphere. Then the mean 
angle at which a collision occurs is #=45°, 
and the mean field strength becomes 2.12 E. 
Thus at the point of contact with the sphere of 
the specific charge 6 the mean field strength is 


(5 a) 


Now, considering a large number of con- 
tacts during each of which a very small pro- 
portion of the charge on the big sphere is 
transmitted to the small sphere, the loss of 
charge per contact may be derived from 
(3 a) as à 

de, df (sb) 


rau ne 
dn dn L 


PACE) 


The loss of charge on one hand is balanced by a 
gain of charge by the small sphere identical 
with (4). Combining this result with (5a) 
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Fig. 4. The numerical value y for contact between two 
spheres of radius R, and Rj. 


Assuming the drop has an initial specific 
charge equal to zero (n =o) 


(6a) 


The specific charge tends towards the limit: 


ee 2.12, ER, (6 b) 
this limiting value being reached after a large 
number of contacts. The number of contacts 
(n) is equal to the number of crystals in the 
path of the falling particle of rime and within a 
certain interception area. This area is somewhat 
smaller than the area Ria of the projected 
sphere (factor of proportionality = y), because, 
as already mentioned above, some of the 
smaller crystals will be caused to deviate from 
their path by the rime. 

Taking into account the aerodynamical de- 
viation of the small sphere when flying towards 
the bigger one, the following procedure may 
be adopted in accordance with ALBRECHT [16] 
and SELL [17]. The deviation in a current 
of air is determined by a numerical factor 
whose physical meaning is determined by the 
ratio of two numbers, namely the path over 
which the small sphere is decelerated divided 


by the radius of the big sphere: 
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SXIR; ze) 
er: ( R, (7) 
where R, =radius of the “small sphere” 

s  =specific gravity of the small sphere 

vjo =relative velocity of the small 
sphere with respect to the big 
sphere with radius R, 


gy om 
8 sec? 
MES x sec 
nm =viscosity of air = 1.61 x 107? LS 
cm 
at — 20° C 


The relative interception-area @ of a sphere 
was investigated by Sell. @ denotes the ratio 
of a distant circular area from which the 
small particles will just reach the sphere, to 
the cross-section of the latter. Sell’s results 
may very approximately be interpreted by 


(8) 


e.g. Ô will have a value of 4 for a crystal (sphere) 
of 20 4 radius and S.G. = 0.9, travelling at a 
relative velocity of 1.6 m/s against soft hail 
of ımm radius. (Only very small particles for 
which 6 is less than 0.1—i.e. particles of less 
than 2.3 w radius in the above example— 
will deviate and not touch the sphere.) The 
relative interception-area is given by @=0.94. 
94 % of the particles travelling in a cylinder 
of the same cross-section as the sphere will 
collide with it. Thus the 20 u-crystals will fly 
towards the sphere on a nearly rectilinear path. 

Precipitate of radius exceeding ı mm, falls 
more and more slowly than in proportion to 


p=tanh 0.43 6 


TE; Ces 
its size, in other words, because En in (7) 


diminishes with growing R, the vol 
therefore the relative interception-area will 
become smaller. For hail with a radius of 
10 mm Ô is 2.3, and the interception-area for 
small particles of 20 u radius becomes 0.76. 

In connection with these aerodynamical 
discussions the mean angle of collision could 
now be investigated more closely for the 
particles hitting the big sphere. This would 
lead to a more accurate determination of the 
numerical factor given as 2.12 in (sa). The 
probability of a collision is not uniform over 
the projection of the sphere. On the one 
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hand the particles will be deviated in their 
flight from the central region of the projection 
and on the other there will no longer be any 
collision at the fringe areas. The angular limits 
are therefore not 0° and 90° but somewhat 
closer together, for instance ALBRECHT’S in- 
vestigations [16] allow their determination 
in the case of a cylinder, and Albrecht gives the 
diameter of that section of the cylinder with 
which the small spheres just collide. For the 
very small spheres of 2.3 u radius in the ex- 
ample the limiting angle is approximately 
40° for ö=o.ı and 76° for ö=2 under the 
above conditions. This would therefore have 
to be considered in investigating the collision 
of very small particles. 


The rate of electrification 


The rate of charging according to (6a) de- 
pends on the number of ice crystals (small 
spheres) impinging on the particle of rime 
(large sphere). Let the number of crystals 
per cm? in the interception-area pxR?; be N 
and let the relative distance of fall be v,)xt 
in time f. Then the number of ice crystals 
in impinging on the particle of rime equals 
the number of contacts n, where 


n= No Ri nviot (9) 
Combining (9) and (6 a), the specific charge 6 


of the rime-particle is 


= 2,12 (eo Te 1) (10) 


where the time constant T of the charging 
process is given by 


I 
T= 
(a-y-~) x NR} x 110 (1) 


The numerical constants y (eq. 4) and y (eq. 8), 
and relative velocity v,) depend on the size 
of the rime. 

The time constant may be illustrated by a 


‘numerical example. Let the radius of the ice 


crystal be Ry=50 wu (that is equal to 5 x 10° 
cm), the number N per cm? be 0.2, and the 
radius of the rime be 2 mm. Thus according 


to fig. 4 a becomes 0.025 and y becomes 


26; According to (8) @ is 0.95, and as the 
relative velocity v,, is about 230 cm/s, 

Tellus VI (1954), 4 

5— 410071 


if 
T= 
(n x 1.6 x 0.95) 0.2 X 52x 10-8 x 230 190 sec 


i.e. the time constant is therefore 3 min 10 sec. 
For crystals (spheres) which exceed 20 u in 
radius the product yg is about 1.5 and therefore 
the time constant is given approximately by 


(11 a) 


(ype oe Sen Be 

NR XU a6 
on the assumption that the number N of ice 
crystals per cm? is invariant with time. This is 
to be expected if the rime corn remains in 
its position, in other words, if its terminal 
velocity when falling, is equal and opposite 
to the velocity of the convection current of 
the air. If this were not the case, a conventional 
time-constant would not suflice to describe 
this process. 

The relative distance travelled by the ice 
crystals with respect to the rime-particle 
in a time equal to that given by time constant 
T shall be called the specific relative rise H. 
Let the relative distance travelled in time T 
by the first of the crystals impinging on the 
rime be h. Then (eq. 10) gives the specific 
charge 


(10 a) 
where 


I 
Be (x:y:p)xNR; (13) 


Consequently the specific relative rise depends 
primarily on the number of ice crystals per 
unit volume and the square of their size. 
Probable or possible values must now be 
assigned to these as no measurements are 
known. 

The number of ice crystals per cm? that 
can form under certain atmospheric conditions 
depends largely on the temperature, but is in 
no way fixed by it. At the same temperature 
there may be a variation of several orders of 
magnitude depending on the origin of the 
air. From numerous investigations on the 
number of ice crystals per unit volume the 
measurements of AUFM KAMPE and WEICK- 
MANN [2] and V. I. SCHAEFER’S researches 
[18] will be used hereafter. In artificially 
produced clouds ice crystals form only sporad- 
ically (some crystals per m?) at temperatures 
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Fig. 5. Relation between number of ice crystals and 
temperature. 


between -3° and -10°. Occasionally only 
supercooled water droplets and an absence of 
ice crystals may be observed during laboratory 
tests on cloud air down to —18°. In the 
temperature-region from —10° to —25° the 
number of crystals increases by two orders of 
magnitude for every decrease of -5°; i.e. 
if the temperature decreases from -10° to 
—25° the number of crystals increases by about 
5 orders of magnitude. 

Fig. s shows the data which will be used in 
the following. 

The size of ice crystals will now be investigat- 
ed in relation to their determined conditions, 
e.g. suppose an ice crystal originated at — 10° 
or -ı5° and travelled at a fixed vertical 
velocity, say 5 m/s, through cold regions of 
known air-velocity till at a definite time it 
crosses, ‘for ‘instance, ‘the —25° or =30° 
isothtrm. Further simplifying assumptions 
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must be made. Firstly let the forming ice 
crystal be spherical growing from zero-radius. 
On the other hand ice crystals can originate 
by freezing of supercooled droplets; the 
spheres of radius Ry which may therefore be 
present at the beginning of the process, will be 
neglected. Further assumptions are as follows: 


1) Diffusion will be considered on stationary 
spheres. If the radius does not exceed 100 y 
the calculated diffusion is not influenced 
by the motion of the sphere. 


The temperature of crystals is assumed to be 
equal to the temperature of the surrounding 
air. The influence of condensation heat on 
the temperature will be disregarded. The 
decrease of temperature of an ice crystal 
traversing a temperature field of constant 
gradient in a convection current of 5—15 
m/s is invariant with time. In the numerical 
example given below this temperature de- 


| baa 


crease is one degree for between 10—30 


seconds. 


The vapour-pressure of the supercooled 
droplets determines the vapour-pressure in 
the cloud. The known difference between 
vapour-pressures of water in solid and 
liquid states is the effective difference of 
vapour pressure between ice crystal and the 
surrounding air. 


The value of Ry necessary for eqs (11) and (13) 
will be too great as a result of this simplifica- 
tion. Because the number N of crystals per 
unit volume varies at a given temperature 
between extremely wide limits, as seen above, 
there will in any case be certain limits for the 
wanted value NR’. This makes the discussion 
of the results of the investigation possible. 
Finally it has to be considered that the radius 
R, of an ice crystal can change discontinuously 
namely if an ice crystal collides with a super- 
cooled droplet. Droplets with a radius from 
60 u to 100 u are to be found in cumolo- 
nimbus before reaching the 0° isotherm [2, 
p- 67]. If they are supercooled to a temperacure 
of -10° they solidify within ı to 2 seconds 
after the absorption of an ice crystal, which 
can easily be calculated. 
The increase of matter by diffusion on a 
sphere with radius Ry is 
dm AS MApR 
dt Rie 


(14) 
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where 


D=diffusion-coefficient. For water vapour 
in air at 16° C=289 K° D=0.285 cm?/s, 
D being approximately proportional to 
the absolute temperature 


M=molecular weight (18 g/mol) 
Ap=difference of vapour pressure (kg/cm?) 
R=gas constant (84.8 kg cm/mol degree) 


T,=absolute temperature (°K). If the above 
value of D is used, T,=289° K and 
the expression becomes independent of 
absolute temperature. 


The difference of vapour pressures A p can 
be calculated exactly. Table II gives values 
published by M. RosıtzscH [19]. Because 
these values are unwieldy for an integration 
they will be replaced by an empirical rela- 
tionship 
proto Je pts T_ (ko lem?) (rs) 
where T_ denotes degrees below zero (-° C). 
Table II shows a comparison between exact 
values and those calculated empirically accord- 


ing to (Is), and they show reasonable agree- 
ment with each other. 


Table II 


difference of vapour pressures PA 


a: 
42 exact calculated by (15) 
| mm merc. mm merc. | kg/cm? 
3 0.104 0.087 7.19% 107% 
6 0.163 0.153 2.08 
12 0.198 0.198 2.69 
18 0.175 0.176 2.40 
24 0.134 0.132 1.80 
30 0.092 0.094 1.275 
36 0.061 0.062 0.85 


‘In eq. 14 mass (m) is replaced by volume x 
specific weight s, and the element of time dt 
by a temperature increment dT_ because 
the ice crystal is carried through a temperature 
field of gradient h (deg/cm) with a velocity v. 


From (15) 
ROlRo= A: Tite T JT (ré) 
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where the constant 


_ DM:4x1075 
sRT hv 


(16 a) 


For h=6.67 x 1075 deg/cm and v=soo cm/sec 
A becomes equal to 0.278 x 10%. If the crystal 
originates at 1°, its size at temperature T_, 
is given by 
Tes 
R2=2A4f f(T_)aT- (17) 


Tax 
The integral may be expanded into an infinite 


series 


z 


DAs Reps = 


| (18) 


where z=0.1068 T_. Some numerical values 
of the integral are given in table III. 


ff 1 D2)d T= =o. 1068-2) 2223 


2. ze 
ese 315.28 = 


Table III 
Pec | 5 | Io | 15 | 20 | 25 | 30 
LE 
if 1210, |e On lee7a.o) 9.104. | er28) 7749 
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Figure (6) gives the radius Rg of a spherical 
ice crystal originating at the temperature indi- 


210" 


LC 
Fig. 6. Relationship between the size of ice crystals 
reaching the —15, —20, —25, and -30° C isotherms 
and their temperature of formation. 


4500 4000 3500 3000 m 


Fig. 7. Relation between temperature and (number of 
ice crystals) x (radius)?. 


cated by the abscissa; and which crosses the 
—15°, —20°, —25°, and -30° isotherms-at 
the three given velocities. Thus an ice crystal 
originating at —15° has a radius Ro of 40.9 u, 
28.9 u, and 23.7 u at respective convection 
velocities of 5, 10, and 15 m/sec on reaching 
the —20° isotherm; in consequence the time 
interval is 150, 75, and so seconds. The tem- 
perature gradient h =6.67 x 1075 deg/cm corre- 
sponds to the decrease of temperature with 
height in thunderclouds at temperatures below 
—10°. This value considers the entraining of 
the air surrounding the cloud and is conse- 
quently a bit greater than the decrease of 
temperature corresponding an adiabatic curve 
of damp air. 

From R° and the number N of ice crystals 
depending on the temperature (figure 5) 


AE. 
[Ran 


is obtained by graphical integration as a func- 
tion of temperature. These values of NR? are 
shown in figure 7. They hold for a vertical 
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velocity of 5 m/s and the greatest as well as 
the smallest number N of ice crystals. These 
values are inversely proportional to the veloci- 
ty of the crystals carried through the tem- 
perature field, and halved at 10 m/sec. 

The size distribution of the ice crystals 
according to figure 8 was calculated at -28° C 
and 16.6 crystals per cm’. The radii are tabu- 
lated in steps of 5 u. The number AN and the 
ice content per cm? is indicated respectively. 
The total ice content is 0.686 g/cm and this is 
small compared to the water content of 
supercooled droplets. The assumption that 
the vapour pressure of the supercooled droplets 
determines the rate of growth is reasonably 
substantiated. Figure 7 can now be interpreted 
as follows: If a particle of soft hail is floating, 
i.e. if its terminal velocity when falling is 


equal and opposite to the vertical velocity of | 


convection, it passes an air current carrying N 
crystals per cm? with a surface of |R;} |. In this 
case the appropriate values of the specific 
relative rise H and the time constant T are 
given at the right hand ordinate. The mean 
values of NR} at —22° C give a time constant 
of 5 minutes. These values give some idea 
of the speed with which electrification can 
take place. Moreover the heights above the 
o° isotherm are given on the abscissa for 
which a temperature gradient of 6.67 x 1075 
deg/cm is appropriate. 


Discussion 


The results so far obtained will now be 
checked against some measurements compiled 
by Ross Gunn [20]: 


1) The maximum charge in a drop is +(20 
= 30) ESU per gram (measurements at 
ground level) 


2) The mean charge of a drop is approxi- 
mately proportional to + 2.8 R® The 
highest charge is about + 10 R?, where 
R is the radius of the liquid drop (measured 
at ground level) 


3) The charge per drop is about + 0.27 ESU 
at 2,300 m altitude (+10.3°C), -0.135 
ESU at 3,800 m altitude (+2.4° C), and 


+ 0.052 or —0.062 ESU at 6,600 m altitude | 


(—9.9° C). 
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4) Only occasionally does the field gradient 


exceed 500 — 
cm 


The maximum field strength of 9 thunder- 


in large cloud regions. 


V 
storms had a mean of 1,300 —. 
cm 


According to (6b) the final charge of 
rime is: 
e=-2.12 ER, 


where R, denotes the radius of the rime. If 
we replace R, by the radius R, of the liquid 
drop of equal weight, then R’xs=R?. Ac- 
cordingly the charge per gram of water is 


0.001 
CeO 2072072507507 60'51 
; R 
Fig. 8. Size-distribution of ice crystals and their water 
content on reaching the — 28° C isotherm in a convection 
current. Number of crystals per cm® 
N =16.6 
NR? =68.9 10 ® em 
water content #  —0.686 g/cmÿ. 
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In table IV charge per gram as well as the 
charge in terms of (drop radii)? are compiled 
for various field gradients. The figures are 
based on a mean drop radius of 0.3 mm 
which may be taken as representative for 
thunderstorm rain near the earth’s surface. 


Table IV 
held Mean drop radius = 0.8 mm 
BAGUE soft hail (s=0.6) |soft rime(s—0o.125) 
V — charge in ESU’s | charge in ESU’s 
cm 
per gram| x FR? |per gram| x R? 
350 10.5 3.48 29.5 9.90 
500 14.8 4.98 42 14.15 
1,000 29.6 9.96 84 28.3 


The figures given by Gunn: 30 ESU/gram 
and 10 R? per drop are confirmed for soft 


; : V 
hail at 1,000 — and for soft rime at 350 —- 
cm cm 


The charge per drop for which Gunn gives a 
maximum figure of +0.27 ESU is used in 
table V to compute the drop radius. This 
value refers to very heavy rain. 


Table V 


drop-charge - 0.27 ESU 


field Rain drops originating from 
gradi- > 7 
ent soft hail (s=0.6) | soft rime (s=0.125) 
i DE enti Me oe ba 
aa radius (mm) weight radius (mm) weight 
solid fliquid| ”® | solid |liquid | > 
350 | 3.31] 2.801 -92 Bead 1.66 19.1 
Booalw2.77, 23a 5307 zeny 1.385 Ti tS 
1,000 1.06| 1.65| 19.0 1.96 0.98 3.96 


Droplet weights of 10—20 mg are rare but 
peculiar to regions of very heavy rain. Thus 
Gunn’s experimental results fit in with the 
calculations as far as negative charge is con- 
cerned. The following process applies to 
positively charged drops. Small positively 
charged ice crystals recoiling from rime after 
contact will finally cling to rime, e.g. to a 
growing particle of soft rime, thereby con- 
veying their positive charge to the soft rime. 
The ice crystals are then frozen in by super- 
cooled water droplets. A small percentage 
of charged ice crystals in the soft rime “skele- 
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Table VI 


Lightning | <—0'—15’—> Z 4— 25 —40" 


| Mature stage | 


15 —30 — | + 


| Cumulus stage | 


ro rs —> 


Convection 
A 


Precipitation begin 
Ü 

Initial radar echo 

.[. -4° or -21° .|. —24 


No radar 


A 
Ra 
Le] — 


2.5 


ton” is sufficient to give the soft rime the same 
charge held by soft hail. Thus the charges 
will finally be distributed in such a manner, 
that positively charged rime will exist beside 
negatively charged rime in a layer following 
on the first negatively charged layer, i.e.: 


1. layer: negatively charged soft hail or small 
smooth hailstones, 


2. layer: negatively and positively charged 
rime with a net charge approxi- 
mately nil when taken over a wider 
space, but possibly having different 


mean rates of fall. 


3. layer: positively charged ice crystals. 


It will now be attempted to incorporate this 
process in the chronological development of 
thunderstorms that can be derived from the 
results of investigations mentioned above. 


- Single turret thunderstorm 


To the extent table VI is based on the state- 
ments of Byers, it is valid for a temperature 
of 20°C and a pressure of 900 mb at the 
bottom of the cloud. The o° isotherm lies 
within the cloud at 4,900 m. Rain reaches the 
ground 10—15 minutes after the origination 
of the radar-cloud. No rain can reach the 
ground as long as the updraft continues, 
neither can soft rime and soft hail—which are 
characteristic of thunderstorms according to 
J. Kürrner [21]—appear at the altitude of 
the 0° isotherm. A downcurrent is a prere- 


+— 273 mir $—15 —> 


Maximum rainfall 


rth’s surface | | Downdraft over entire cell | 


® 
infall at ea 
® 


_—_—, = 


| Dissipating stage | 


30° — 


Radar top: —20° and lower or: -—28° ./. —29° 
Mean high: 11000 m 


quisite to precipitate reaching the earth, at 
least in that part of the cloud from which the 
precipitate originates. 

A crystal—whose growth by condensation 


Was shown in figure 6—grows from approxi-, 


mately 60 u radius taking up supercooled 
droplets, considerably faster than by conden- 
sation. It takes a relatively long time, though, 
before soft hail will descend to the ground 
in a convection current of 5 m/s. F. H. Lup- 
LAM [22] has investigated the growth of 
crystals partly by condensation and partly by 
coagulation. To simplify matters it will be 
assumed that precipitate of 60 radius will 
grow by coagulation with droplets of at 
least 15 u mean radius, and in this connection 
the process of condensation may be neglected. 
Let unit volume contain w grams of water 
as small droplets whose rate of fall may be 
neglected compared with that of the growing 
precipitate. Then 
4 sıdR,=w-@-v,dt (20) 
Under these assumptions the relative inter- 
ception-area according to (8) is practically ı 
for water droplets and less than ı for hail and 
soft hail. Thus 


a po 
ts wov, 


10 


dR, (21) 


This leads to the time interval to, in which 
the precipitate of spec. weight=s, grows from 
radius Rjo to Ru. The rate of fall v, of the 
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precipitate can be expressed by Schmidts 
formula: 
10° ; 
Nee B (22) 


— Es 
FER? 


The constants A and B, valid for an altitude 
of 6,000 m, are given in table VI. 


Table VI 


Hail: The values are valid for irregular forms. For 
regular forms reduce B with a factor 0.7. 


A | B | Sy 
dropse Fosse: 0.682 370 1.0 
OUP hal... 0.576 781 0.83 
Sena Nail. RES ar 0.317 1,350 0.60 


The time fo, is given for an interval, over 
which w, s, and m may be regarded as in- 
variant, by 

Ru 


45 A 
t, = —— — ——_ 1/2 —6 
nase | pe BR} i 10 (23) 


During the time #9 the growing precipitate 
rises through a height S: 


th 


45 
S= vias J vadt=vtq,— 72 (Rii— Rio) (24) 


0 


where v denotes the velocity of the convection 
current. As a numerical example, the rise is 
related in figure 9 to the time for droplets, 
hail, and soft hail growing by coagulation in 
an atmosphere containing 5 x 10° gram water 
per cm? as small droplets. Rise and time are 
reversely proportional to the water-content. 
The weights of precipitate in figure 9 are 
indicated alongside the curves. A water 
droplet rises about 500 meters from where 
its size corresponded to 60 y radius. Then it 
falls and, on exceeding a weight of 100 mg, 
it may disintegrate and initiate a new process 
with a multitude of droplets, dropping only 
slowly below its origin as long as the con- 
vection current persists. Within the first 15 
minutes taken from its beginning as 60 u 
precipitate neither hail nor soft rime drop 
below their place of origin, but this is con- 
tradicted by experience. Therefore the down- 
drift of the air must start soon after the first 
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1000 
Fig. 9. Height to which precipitate rises while growing 
by coagulation in a convection-current. With 

SS 1001. KG SSO > aay 
(A droplets, B hail, C soft hail, irregular form). 
W508, Sem? 
Weights are denoted by numbers, for C very approxi- 
mate. 
R=beginning of radar sight. 


strong precipitate is formed. According to 
scheme p. 12 this might take place after 5 
minutes. The downdrift of air is propagated 
to the upper regions, so that a soft hailstone 
originating in the convection current will 
hardly exceed an altitude of 1,500 m above 
its origin. The timing would then be as 
follows: A soft hailstone originating at - 5° C 
and growing in the vertical air current by 
condensation will have a radius of nearly 
60 u when it reaches the -ı5° C isotherm 
and its weight has increased to 10,000 times 
this initial value, when it reaches the — 25° C 
isotherm, 1,500 m above its ‘‘coagulation- 
origin”. The size of the precipitate and 
number per unit volume determine the 
possibility of radar soundings. It is assumed 
that the radar cloud appears 4—6 minutes 
after the beginning of the 60 u precipitate. 
This fixes the rise at 300 to 600 m for droplets, 
and at about 1,200 m for soft hail. According 
to the above scheme these would be the 
limits for the beginning of the “cumulus 
stage”. 

During the next 10°—15 minutes the down- 
ward movement of the air containing super- 
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cooled water droplets and ice crystals will 
increase more and more. Even.now contact 
between ice crystals and soft hail takes place 
in the same way as described above, but the 
effective area originally lying between the 
— 18° and — 26° isotherms will move down to 
higher isotherms. The generation of charge 
will cease when the ice crystals have dis- 
appeared. This limit is given by their subli- 
mation when the vapour pressure of the down- 
drifting air in the dissipating stage sinks 
below the vapour pressure of the ice by en- 
trainment with the surrounding dry air. 
Elsewhere [23] it was shown that in a 
region where charge is generated the electric 
field grows logarithmically, i.e. at first slowly 
and then very rapidly, according to a process 
intensifying with increasing field strength. 
If applied to the above case this calculation 
leads to plausible assumptions for the ap- 
pearance of a strong field within a short 
time. Assuming, for instance, that soft hail 
of 2.3 m/sec velocity drops within about 3.5 
minutes through a soo m thick layer of uni- 
form space-charge, thereby increasing this 
charge, then the field will be increased to 
1,500 times its initial value, if the rate of 
precipitation is 1.25 mm water per minute 
and the time constant is 5 minutes. At alti- 
tudes in which charge is generated the carth’s 
electric field is very weak: this is intensified 
immediately on the formation of a first heavy 
precipitate, because the small cloud-droplets 
have a slight negative charge leading to the 
formation of an initial field of about ı V/cm 
within approximately 3 minutes. [23] Conse- 
quently a strong field could be built up 
within 6.5 minutes. This is in good agree- 
ment with the statement of Byers [2, p. 62], 
who had observed no electric field when 
first flying through a thundercloud, but 8 
minutes later he found a strong positive field 
associated with heavy precipitation. 
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In this connection a possibility may be 
pointed out to account for the third and smaller 
region of positive space charge, mentioned 
above, which is between the + 10° and o° 
isotherms. This positive space charge does 
not always appear, and it is confined to a 
relatively small region. Its effectiveness may 
be roughly represented by a sphere of 500 m 
radius [1]. The space charge is relatively dense. 
According to G. C. Simpson’s soundings the 
mean density is 2.4x10-5 ESU/cm?. Ac- 
cording to Ross Gunn [20], it is 15.5 x 1075 
ESU/cm?. Ross Gunn who obtained this 
figure by using an aeroplane gives about 
0.5 ESU per drop on the + 10.3° level. Nega- 
tively charged hail having passed a region 
of negative charge will then enter a region 
with a negative field gradient. At temperatures 
above 0° hail develops a moist and slippery 
surface. This makes it possible for small 
water-droplets—e.g. having a radius of 15 u 
—to slip off the corn-shaped bottom of the 


hailstone thus leaving it again. Assuming that” 


the hailstone does possess a conical shape, 
these considerations may be transformed for a 
spherical shape. Water droplets will then 
leave the sphere at an angle # to the vertical 
axis, the final charge of the hailstone will be 


(25) 


ea = 3, cos PER? 


ie. it will equal 0.52 when 9 =80°. The time 
constant of this process is very small. As- 
suming small water-droplets amounting to 
some tenths of a gram per m? its order of 
magnitude will be 10 sec which is sufficient 
to discharge the hailstone and re-charge it 
with the opposite sign. This would increase 
the existing negative field. These considera- 
tions can, however, not be extended by cal- 
culations without further experimental in- 
vestigations. 
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Global Radiation Resulting from Multiple Scattering in a Rayleigh 


Atmosphere* 
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Abstract 


Using Chandrasekhar’s solution for the radiative transfer problem in a plane-parallel model 
of the earth’s atmosphere, under the assumption of perfect, conservative scattering according 
to the Rayleigh law, expressions are derived for the relative global and sky radiation received on 
a horizontal surface, as a function of normal optical thickness and inclination of incoming 
parallel radiation. The corrections representing the effect of ground reflection are also expressed 
m terms of certain functions already defined by Chandrasekhar. These expressions, which 
include the effect of all orders of scattering, are then used to compute the absolute global and sky 
radiation as a function of wavelength and solar zenith distance, based on the solar extraterrestrial 
energy curve as given recently by Nicolet. The results show that, under these assumptions, the 
global radiation received at the surface of the earth should remain essentially constant in spectral 
distribution for a wide range of solar altitudes, and that the pure sky radiation should have two 
maxima, centered at 3,300 À and 4,100 À respectively, for solar zenith distances ranging 
between o° and 45°. 

The theoretical work is compared with Bernhardt’s results obtained by means of certain 
simplifying assumptions about the successive orders of scattering. 

The absolute global and sky radiation is integrated over the range 0.29 u to 4.00 u, and it is 
shown that in general the computation is in good agreement with existing measurements 


By DIRAN DEIRMENDJIAN and ZDENEK SEKERA, Department of Meteorology, 


under clear sky conditions. 


I. Introduction 


A knowledge of the amount and of the 
spectral distribution of the radiation received 
by the earth’s surface is essential in all problems 
dealing with the radiational effects in the 
atmosphere. Several attempts have been made 


* U.C.L.A. Department of Meteorology, Papers in 
Meteorology No. 23. This research was supported in 
part by the United States Air Force under Contract No. 
AF 19 (122)—239, monitored by the Air Force Cam- 
bridge Research Center. 


in the past to evaluate this quantity from 
theoretical considerations. The main difficulty 
encountered in all of these arises from the 
complicated nature of the radiative transfer 
in the atmosphere, which, in the most general 
case, is rather inaccessible to mathematical 
solutions in explicit form. Thus the efforts 
have been mainly limited to the evaluation of 


this radiation in certain simplified models of 


the earth’s atmosphere, assuming Rayleigh 


scattering only. But even in this case the dith- 
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culties connected with an exact evaluation of 
the effects of scattering of orders higher 
than the first are such that only the secondary 
scattering could be included so far (cf. L. V. 
KING, 1913) and even that in the form of 
approximative solutions only. 

Only very recently however, an ingenious 
method developed by S. CHANDRASEKHAR, 
(1950) especially suited to problems of radia- 
tive transfer in stellar and planetary atmos- 
pheres, has made it possible not only to in- 
corporate, in an exact solution, all orders of 
scattering, but also to introduce the effects of 
ground reflection. In the present study, 
Chandrasekhar’s exact solution is used in 
order to derive an expression for the total 
solar energy received by a horizontal unit 
area of the earth’s surface, in such a wa 
that the relative importance of the sky radia- 
tion and the radiation due to earth reflection 
can be appreciated accurately, as a function 
of wave length and sun elevation. Numerical 
computations are then carried out using this 
expression and the scattering functions in- 
volved, which were available to the authors 
(SEKERA, BLANCH, 1952; SEKERA, ASHBURN, 
1953). Finally, these results are applied to a 
recent determination of the sun’s extraterres- 
trial energy curve in order to obtain—perhaps 
for the first time—the exact theoretical spectral 
distribution of solar energy reaching the earth’s 
surface after passage through a perfectly scat- 
tering atmosphere. As will be seen below, the 
theory, in so far as it can be compared with 
observation, gives remarkably consistent and 
accurate results, notwithstanding the highly 
idealized model used. 


2. Formulation of the problem and its 
theoretical solution 


Because of its small depth as compared to 
the earth’s radius, for the purposes of our 
problem the atmosphere may be represented 
by a plane parallel one of infinite lateral ex- 
tent but of finite depth. The density will be 
assumed to depend only on height above the 
ground, varying in accordance with mean 
conditions in the actual atmosphere. Further- 
more, because of the large distance of the sun 


from the earth, the solar radiation will be 


represented by parallel radiation, of a known 
spectral distribution, incident on the topmost 
layer of the atmosphere. The problem then 
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consists in finding the inward flux of radiant 
energy through a unit area of the earth’s 
horizontal surface, as a function of the zenith 
distance of the sun, and of the wavelength 
and intensity of the incident radiation. 

The problem can be simplified by separating 
this net flux into two distinct types of radia- 
tion: (a) a flux of parallel, direct solar radiation, 
modified only by the transfer through the 
atmosphere, called sun radiation; (b) a flux 
of non-parallel, diffuse radiation originating 
from the various scattering centers of the 
atmosphere, such as gas molecules, dust 
particles, etc., coming from the whole sky 
and called sky radiation. At present this sky 
radiation can be computed only under the 
assumptions (i) that the scatterers are of 
negligibly small dimensions when compared 
to the wavelength of the radiation, (ii) that 
all the radiation received by the scatterers is 
reradiated in all directions without modifica- 
tion of the original wavelength (ic., case of 
conservative perfect scattering). These are 
the assumptions under which the well known 
Rayleigh theory of scattering (STRUTT, 1871; 
RAYLEIGH, 1899), is applicable, which provides 
the laws governing the intensity of the scattered 
light as a function of direction and wavelength 
as well as its state of polarization. 

The sun radiation is proportional to the 
so-called reduced flux. The incident solar 
energy is simply assumed to be attenuated or 
reduced by the amount scattered in all di- 
rections by the scatterers, without, however, 
losing its character of a parallel radiation. If 
nF, be the net flux of monochromatic 
radiation of wavelength A, incident at the top 
of the atmosphere under the zenith angle ©,, 
then the reduced flux zF, through a unit area 
normal to the incoming radiation, at a height 
z, above the earth’s surface, can easily be 
shown to be 


nF, = nFye tH (1) 


where u, = cos @y; the non-dimensional pa- 
rameter tT, called the normal optical thickness, a 
function of the wave length À, stands for the 
integral 


t= f bide (2) 


where BB; is the monochromatic attenuation 
coefficient, a function of the density and of the 
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refractive index of the air, and z is the distance 
measured vertically above the earth. For a 
given distribution of air density with height, 
the normal optical thickness (2) corresponding 
to Rayleigh scattering can be determined 
(DEIRMENDJIAN, 1952). Then, at any fixed 
level z cm above sea level, the optical thickness 
becomes a function of the wavelength only, 
and can therefore serve as the dimensionless 
parameter defining the wavelength of the 
incident radiation. 

The net flux of parallel radiation through a 
horizontal unit area of the earth’s surface, 1.e. 
the sun radiation at the wavelength A, will 
be given by the expression 


Sa = nF, Mo — LF gage “/# (3) 


which reduces to (1) when the sun is exactly 
at the zenith (©, = 0). 

The sky radiation, as defined above, repre- 
sents the flux of the radiation scattered by 
all the scattering centers in the atmosphere. 
Not only is the direct solar radiation itself 
scattered, but also part of all the radiation 
reaching the earth’s surface which is reflected 
back into the atmosphere. As will be seen 
below, under certain assumptions this re- 
flection will result in an increase in sky ra- 
diation which can be computed separately. 

Considering first the standard problem, 
i.e., the one without reflecting boundaries, 
the monochromatic radiation received from 
an elementary cone of solid angle dw and 
from a direction making the angle © with 
the vertical will contribute to the net flux 
through a horizontal unit area by the amount 
I cos © dw, where in general, I = I (Oo, Mo, 
©, y) is the specific intensity of the diffuse sky 
radiation. The sky radiation is then given by 
the integral 


H; = f Ios Odo (4) 


where the integration is to be performed over 
all the solid angles of the outer hemisphere only, 
i.e., it is assumed that the surface element is 
receiving radiation from one side only. Since 
dw = sin OdOdp (p denotes the azimuth of the 
particular direction), the expression (4) can 
also be written as 


1 2% 


D | Sf Tudu de, (u =cos0) (5) 
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The diffuse sky radiation is in general polarized 
and thus the specific intensity I can be ex- 
pressed as the sum of the intensities of two 
components polarized in two directions per- 
pendicular to each other, for example, nor- 
mal (J,) and parallel (}) to the vertical plane 
through the direction given by the zenith 
angle © and the azimuth y. From the general 
equations of radiative transfer (with all or- 
ders of scattering included) according to the 
method described by CHANDRASEKHAR (1950, 
1951), the expression for the intensities I, and 
I, may be written in the form 


I; = 1} + A; (u, Mo) cos (p — Po) + 
+ B; (u, Mo) cos2 (P— Go), (i = br) (6) 


where u, and y, specify the direction of the 
incident radiation and the functions A; and 
B; are linear combinations of six special func- 
tions of u and u, only (CHANDRASEKHAR, 1950; 
SEKERA, BLANCH, 1952). If the expression (6) 
is substituted for the components of J in (5), 
ie. placing I=1, + I, it is evident that the 
azimuth dependent terms—i.e. the terms con- 
taining cos (p—@p) and cos 2 (y — 9,)— will 
vanish upon integration with respect to 9. 
Thus only the azimuth independent functions 
I; and I} of (6) will contribute to the inte- 
gration in (5), as might be expected from 
symmetry considerations. These functions J} 
depend on the flux of the incident radiation 
according to the equation (cf. CHANDRA- 
SEKHAR, 1950, pp 265, 276) 


Wen) m 


where T° represents the transmission matrix 
with elements Ti (i, k = 1, r). Hence substitu- 
tion of (7) into (5), with the assumption of a 
neutral incident radiation, i.e. FR= F, = Fy/2, 
yields the explicit integral for the sky radiation 


1 


3 UP, 0 0 0 0 
WM = aoa ( orte Tr + Ty + 1%) du (8) 


or with the use of the functions y; and ves 


conveniently introduced by CHANDRASEKHAR 
(1950, p. 276), 


I I 
Hy = HF to} 51 (u) + 57 (ita) = er" | (9) 
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A comparison with (3) shows immediately 
that the expression ‘ 


TE 
5 Foto [yi (to) + Yr (Mo)] = Ha + Sz (10) 


represents the monochromatic global (sun + 
sky) radiation where the functions y; (jg) /2 
and y,({9)/2 have a clear physical meaning: 
They represent the monochromatic global radiation 
resulting from the components polarized parallel 
and normal respectively to the vertical plane, for 
unit incident flux (2F = 1) of neutral radiation, 
when there is no reflection and ©, = 0. 

In the expressions (9) and (ro) the effect of 
the ground reflection has not been taken into 
consideration. In general, such a reflection is 
equivalent to an additional illumination of 
the atmosphere from below, hence a corre- 
sponding increase of the diffuse sky radiation 
should be expected. In the special case of 
reflection according to Lambert’s law*, this 
increase can be expressed in terms of the same 
y; functions mentioned above. The addi- 
tional intensity 1” +1f of diffuse light due to 
the ground reflection is also azimuth inde- 
pendent and has the form (CHANDRASEKHAR, 
ELBERT, 1951) 

AusF, 
T+ T= a Write) Fra) 


[2 — y) — 7, (u)] (11) 


where A is the albedo and s is a quantity de- 
pending upon + only, but independent of fu, 
connected with the y; functions by the relation 
(cf. CHANDRASEKHAR, 1950, p. 280) 


s=1— ff) ty) ndu (12) 


The sky radiation resulting from ground 
reflection only, can now be easily computed 
by using the expression (11) for the intensity I 
in (5). The integration, carried out by using 
(12) for s, yields the result 


(far 1) ududp = 


Eu. As 
= nu) + du} (3 


* i.e. the reflected radiation is assumed to be unpolar- 
ized and independent of direction, its total outward flux 
being always a constant fraction A, (0 S$ A < 1) of the 
incoming flux, called the albedo of the reflecting surface. 
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and if this is added to (9) one obtains a simple 
expression for the sky radiation which in- 
cludes the effect of the ground reflection, 
namely 

|yi(to) + Yr(Ho) 


Hy = SER lie teats =) 


_ - un! (14) 


Finally, when (14) is in turn added to the 
expression (3) for the sun radiation, a still 
simpler expression results, representing the 
monochromatic global radiation G; when 
the carth is reflecting with an albedo A, i.e., 


L (flo) + Yu 
Gum Hit Sim Fuge) * Ft (rs) 


Thus, the global radiation can be easily 
computed provided the functions y;(u) and. 
yı(u) are known. The evaluation of these 
functions is not so simple, however, requiring a 
solution of two systems consisting each of 
two non-linear integral equations for two 
pairs of functions (X;, Yı, X,, Y, in Chandra- 
sekhar’s notation). The solution of these 
equations can be obtained by successive itera- 
tions, and has been performed (SEKERA, 
BLANCH, 1952) for three normal optical thick- 
nesses, namely T=1.00, 0.25, and 0.15 corre- 
sponding at sea level to the wavelengthso.3 121, 
0.4364 u, and 0.495 u respectively. During 
the computation, it was found that the correct- 
ed second approximation, corresponding phys- 
ically to the inclusion of secondary scattering 
only, gives values of sufficient accuracy for 
T £ 0.10, so that the y; functions in this range 
were computed without iteration by direct 
analytical computation (FRASER, 1952, 1953; 
SEKERA, ASHBURN, 1953). For even smaller 
values of t (i.e. tS 0.01) when the second 
approximation leads to inaccurate numerical 
computations, the first approximation cor- 
responding to the primary scattering only 
provides a sufficient accuracy. In such a case 
the expressions for the intensities I; and I, 
can be easily derived directly from the equation 
of radiative transfer (cf. CHANRASEKHAR, 1950 
p- 44), which for the primary scattering 
reduces to the form 


I , 
1, = Le 9) = Le PU Bitte PIE (19 


where I and F denote one column matrices of 
Stokes parameters for the intensity of diffuse 
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sky radiation and for the flux of extraterrestrial 
incident radiation, respectively, and P denotes 
the phase matrix for Rayleigh scattering [cf. 
CHANDRASEKHAR, 1950, p. 42, equ. (220)]. For 
a neutral incident radiation (FR=F,=} Fo), if 
in (16) only the azimuth independent terms 
are considered, the equations for I} and I} 
asume the form 


di; 


li ae f r) 


=—Ji( 4) (i (17) 


where 
Ji(r, u) 


0 


= 3 —T/Lo ©) — u? > + 
À Fete [a (1— ne) (0-1) 
+u2(t+4,)] 
= + T oT + 2 
Jr (7, 4) 32 Foe [x ml 


The solutions of (17), satisfying the boundary 
condition 


G=lr) 


can be easily computed, and when added, the 
following expression results: 


I (0, U, ®) =} 


Te 


To 


IP +I =- 


ee LoFol3 - 


} } 
eilt — e” Tu 


pi EN 


+ (345-1) 42] (18) 
The sky radiation H;" corresponding to pri- 
mary scattering only is then given by (s) 
with the expression (18) substituted for 1. 
The integration yields the final result in the 
form 


Hj” = = Fe THe {(3 — w>)F4(T, Uo) + 
+ (3 03 - 1) F(t, 0) 


(19) 


where the functions F,, F, are of the form 
I 

à 3 — — 
u 


Ho u 
(20) 


DIE 
70 
(i=1, 2, 3, 4) 


and related closely to the exponential integral. 


5 


(Cf. CHANDRASEKHAR, 1950, App. I, p. 375.) 


ii 
ul [ 


I — EXP 
\ 


F(t, Mo) =H 
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There is no difficulty involved in the compu- 


tation of these functions since they are given — 


by very simple recurrence formulae. For very 
small + the functions F; can be easily computed 
from the series expansions in powers of r. 
The first term in these expansions gives the 
approximation 


((=2; 3,5 


so that (19) reduces to the very simple form 


Hi)== CF oat (21) 


valid for o <1 <<": 


Equations (19) and (21) give the sky radia- 
tion for the case of zero albedo. To get the. 
correction term in the case of reflection, one 
must add to this the corresponding sun radiation 
(3) and multiply the sum by the factor As/(1- 
— As) appearing in (13). If greater accuracy is 


desired, one can examine the behavior of the = 


ratio H;°/H3 as To, which should approach 
unity asymptotically, and then correct the 
H}"’ values accordingly. 

The expressions (3), (14), (15), and (19) for 
the sun, sky, and global radiation contain the 
monochromatic solar extraterrestrial flux 2F a 
as a constant factor. If this flux be taken as 
unity, then the above expressions represent 
what is sometimes known as the relative fluxes. 
On the other hand if the actual extraterrestrial 
flux is substituted for Fo, then the expressions 
(3), (14), (15), (19) will represent the absolute 
fluxes of sun, sky, and global radiation. 

Because of the instrumental difficulties in 
the measurement of the sky and sun radiation 
within narrow spectral ranges, almost all 
the measurements of these quantities, so far 
available, have been made by means of a 
bolometric instrument registering the total 
radiation in all wavelengths. The measured 
values can thus be compared with quantities 
derived from the absolute fluxes of the sky, 
sun or global radiation by integration over all 
wavelengths. These quantities will be denoted 
by S, H or G respectively, with 


s= f'su, H= he Hida, G= fi Gidà (22) 


0 


| 
A 
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and called the integrated sun, sky, and global 
radiation. . 


3. Numerical results 


Relative Fluxes. The relative fluxes S;, H,, 
and G; were computed first and the result of 
the computation is given in Tables I, II, and 
III. The tables cover four values of the variable 
Ho 1.€. 1.00, 0.60, 0.10, and 0.02 corresponding 
respectively to a solar zenith distance of 0°.0, 
53°.1, 84°.3 and 88°.8; and twelve values of 
t, the normal optical thickness. The particular 
value of #,=0.02 (O,=88°.8) was chosen for 
the following reason: In an infinite plane- 
parallel atmosphere, sunrise and sunset condi- 
tions would correspond to a zenith distance of 
90°, in which case the theoretical S, and Ha 
would be entirely unrealistic, on account of 


- the complete attenuation of the incident 


radiation (infinite optical path). However, 
from the dioptrical tables of Link and SEKERA 
(1940), the air mass and thus the attenuation 
factor along the light path in the actual 
spherical atmosphere, when the sun is at the 
local horizon, can be easily evaluated. These 
tables also take into account the bending of 
the parallel radiation caused by refraction. The 
attenuation arrived at in this way should then 
be considered as the maximum finite attenua- 
tion of the incident radiation in the plane- 
parallel model. The corresponding zenith 
distance turns out to be very close to 88°.8, 
which therefore is the value to be used in a 
plane-parallel atmosphere, whenever condi- 
tions equivalent to sunrise and sunset in the 
spherical atmosphere are to be approximated. 

The parameter 7 is used in preference to the 
wavelength for obvious reasons of general 
applicability. The following list shows the 
wavelengths corresponding at sea level to the 
values of t used in this paper, as obtained 
from a recent computation (DEIRMENDJIAN, 


1952): 


387 
I.00 0.25 0.15 
À at 
sea level 0.312 0.4364 0.495 
(microns) 
T 0.06 0.02 325) To 
Wat 
sea level 0.618 0.809 I.000 
(microns) 
T 10.5, OG Hier BAS TOR 2 T6 TO 
À at 
sea level 1.500 2.000 2.500 
(microns) 
T OS TOM OC TO EEK a=? 
Amat 
sea level 3.000 3.500 4-000 
(microns) 


The values for S; can be computed easily 
from a table of the exponential function. 
The computations for G, and H; were carried 
out as outlined in the previous section, ice. 
in the range 0.15 St S1.00 the y functions 
for multiple scattering were used as given by 
the Tables (SEKERA, BLANCH, 1952); for T= 
0.06 and 0.02 the analytically computed y; 
functions (FRASER, 1952, 1953; SEKERA, ASH- 
BURN, 1953) were used, implying the in- 
clusion of secondary scattering only; and for 
(ANGE "3.2510 = T= 8.5 S10" the quantity 


H;” was computed directly from (19) and then 
(1) 


: 5 VE 
corrected by using the ratios = shown below: 
i 


The figures in brackets indicate the extrapola- 

tion for lower t’s, justified both by physical 

considerations and by the fact that the quantity 
H; 


log (1- En) becomes a linear function of 


(1) 
2 


H; 

particular interest since, for the first time, they 
exhibit exactly the contribution of the pri- 
mary scattering to the sky radiation from all 


are of 


log + when t 0. These values of 


H‘)/H, (albedo = 0) 


ae 1.00 | 0.25 | 0.15 | 0.06 | 0.02 | OS TO SE 10 EE ECOLE 
Lo 
| 1.00 .403 734 .812 .904 .960 [.98] [.994] [1.00] 
0.60 .352 699 .785 .887 .95 [.98] [.99] [1.00] 
0.10 270 667 .765 OT .95 [.97] [.99] [1.00] 
0.025 .66 75 .88 .95 [.97] [.99] [1.00] 
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Table I. The relative sun radiation S, as a function of normal optical thickness 7 (see text for corre- 
sponding wavelengths) and cosine of the solar zenith distance 


! 
8.5X |16.5X 2x | 2.15 xX | 1.03% | 5-0 azar 
lo T. |. 1.002 || (0,25 0.15 0.06 | 0.02 ap ae ae Pore ose 2058 oe 
1.00 3679| .7788| .8607| .9418| .0802| .9915 | .09835| .9995 | .9998 | .9999 | .9999 | 1.0000 
0.60 1133] .3955| .4673| .5429| .5803| .5915 | .5984 | .5995 | .5998 | .5999 | .5999 | .6000 
0.10 0000] .0082| .0223| .0549| .o819] .09185| .0984 | .0995 | .0998 | .0999 | .0999 .1000 | 
0.02 0000| .0000| .0000| .0010| .0074| .OIZI .0184 .0195 .0198 .0199 .0199 .0200 


orders of scattering. It is seen that for in- 
frared radiation of À = 2.00 there would be 
less than 0.5 % error if the secondary and 
higher order scattering is neglected. 

As can be secn from Tables I and II, at any 
fixed sun elevation, while the relative sun 
radiation S; increases asymptotically to [Wo 
with decreasing 7 (increasing wavelength), the 
relative sky radiation decreases in the same 
direction, in such a way that for very small 
values of t, H; is proportional to $7 as per 
(21), becoming independent of solar elevation. 
The sum of these two quantities, i.e. the rela- 
tive global radiation G; is shown in Table 
III which, like Table II, is divided into three 
parts corresponding to no reflection and 
Lambert reflection with albedo 0.25 and 0.80 
respectively. The relative global radiation 
generally increases with decreasing t, but 


there is not as big a difference between short 
and long wavelengths as in the case of the 
sun radiation alone. For the case of an albedo 
of 0.80 and u,=1.00 (sun at the zenith) the 
variation with t is almost reversed, so that the 
maximum in G, occurs near T=0.25, at which 
point it exceeds the corresponding extra- 
terrestrial flux by almost 4%, and then 
decreases asymptotically to wu. However, 
when comparing this interesting result with 
conditions in the real atmosphere, it should 
be borne in mind that the model used involves 
a number of unrealistic assumptions: a plane 
parallel atmosphere of an infinite lateral ex- 
tent with each point of its infinite lower 
boundary reflecting according to Lambert's 
law with the high albedo of 0.80. 

These features of the relative sun and sky 
radiation given by the theory are illustrated 


Table II. The relative sky radiation H, for the same values of the parameters as for Table I, showing 


also the values corrected for Lambert reflection with an albedo of 0.25 and 0.80 


al ootlle ss AP .10.568| 20.02 8.5% To ue 2.15%. |11.03x 25.032508 
US 10 Ko 107 10-4 LORS TOM FOR 
Albedo = 0.00 
1.00 .2918| .1096| .0694| .0291| .0099| .0043 .0008 .0003 .0001 .0000° | .0000 .0000 
0.60 ‚2117| .1006| .0659| .0285| .0099| .0043 | .0008 | .0003 | .ooor .0000° | .0000 .0000 
0.10 .0296| .0415| .0371| .0223| .0090| .0041 | .0008 | .0003 | .ooor .0000° | .0000 .0000 
0.02 .0050| .0078| .0086| .0090| .0063| .0035 .0008 .0002°| .0001 .0000° | .0000 .0000 
Albedo = 0.25 
1.00 3748| .1514] .0980| .0422| .0146| .0064 | .00125| .0004 | .ooo1 .0001 .-0001 .0000 
0.60 ‚2525| .1239) .0822| .0363| .0127| .0056 | .ooro | .0003 | .ooo1 .0001 .000I .0000 
0.10 -0334| .0438| .0389] .0234| .0095| .00425| .0008 | .0003 | .ooor .OOOI .000I .0000 
ef Q =. 
0.02 .0056] .0082| .0089| .0091| .0064| .0035 | .0008 | .00025| .ooor .OOOI .000I .0000 
Albedo = 0.80 
1.00 .6590| .2589| .1675| .0724| .0225| .o110 | .00215 .0007 | .0003 | .oo017 .0001 .0000 
0.60 3925| .1839| .1221| .0540| .o190| .0084 .0016 .0005 .0002 .000I .0001 .0000 
Q Q 
0.10 -0461| .0498| .0434| .0258| .o104] .00475| .0009 | .0003 | .ooor .0001 .000I .0000 
0.02 .0078| .0091| .0095| .0095| .0065| .0036 | .0008 | .0002 | .ooo1 .OOOT .000I .0000 
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Fig. 1. Curves representing 
the variation of the log of  !© DETTE gia yee 
relative sun (S,) and sky 
(H,) radiation for given 
solar zenith distance Op, 
with the log of the normal 
optical thickness t. The 
corresponding wavelengths 
would increase from left 
to right on the diagram. 
Dash lines represent the 
quantity As/(r-As) for 
albedo 0.25 and 0.80 (see 
text). 
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Fig. 2. Spectral distribution 
of solar extraterrestrial 
fluxes (topmost full line) 
as given by Nicolet, and 
resulting sun radiation on 
a horizontal surface at sea 
level, for different solar 
zenith distances. The hatch- 
ed areas represent the addi- 
tional radiation received 
from the whole sky due 
to multiple Rayleigh scat- 
tering, and the blackened 
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in the diagram of Fig. 1, which was drawn 
by fitting smooth curves to plots of log H; 
(albedo=o) and log S;, versus log t, taken 
from Tables I and II. The asymptotic behav- 
ior of Ha as ro is clearly seen from the 
corresponding curves. Furthermore, the maxi- 
mum in relative sky radiation is seen to shift from 
shorter toward longer wavelengths as the zenith 
distance of the sun increases. 

Absolute Fluxes: It is well known that there 
is considerable uncertainity regarding the 
exact spectral distribution of the extraterrestrial 
flux of solar energy. A good analysis of up 
to date knowledge was presented recently 
by Marcel Nicolet (1948, 1951 a, 1951 b). 
There exist some more recent direct meas- 
urements of the solar spectrum, especially 
in the ultraviolet and near ultraviolet region, 
made with the use of spectrophotometers 
carried on rockets (JOHNSON ET AL., 1953), or 
by the Bouguer-Langley method applied to 
measurements from ground observatories (PET- 
TIT, 1932; STAIR, 1952). These show certain 
minor differences from Nicolet’s distribution, 
mainly in the near ultraviolet region. How- 
ever, in the absence of more extensive and 
reliable rocket measurements, we have pre- 
ferred Nicolet’s figures which are consistent 
with present knowledge on the solar con- 
stant, and are sufficiently detailed for our pur- 
poses. 
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Nicolet gives values of the net flux of the 
radiation in watt per square meter per 100 A, 
as received at the top of our atmosphere, 
covering the solar spectrum between 2,200 
and 70,000 A (see the uppermost curve of 
Fig. 2, marked Fo). If we use the sea level 
wavelengths corresponding to the twelve 
values of t used in the computation, the rela- 
tive fluxes in Tables I, II, and III would cover 
the spectrum from 3,120 to 40,000 A (this 
range was extended to 2,900 A by a small 
extrapolation toward shorter wavelengths). 

In order to obtain as detailed a spectrum of 
global and sky radiation as possible it was 
necessary to interpolate our numerical re- 
sults for the small wavelength intervals used 
by Nicolet, i.e for intervals of 100 A up to 
10,000 À. In the range 1 u SA<4u, we used 
values at intervals of o.5u as originally com- 
puted, since the solar energy curve is quite 
smooth in this region. The interpolation was 
carried out graphically by plotting H, on 
large scale logarithmic paper; by careful 
curve fitting it was possible to obtain inter- 
polated values to three significant figures, 
sufficient for the desired accuracy. These 
interpolated values of H; added to the Sj 
values computed from existing tables of e* 
yielded the G, values representing the rela- 
tive global radiation. A check was made by 
plotting G, curves and interpolating by the 


Table III. The relative global radiation G,, tabulated on the same basis as Table II 


er 85x | 16.5x152%x | 2.15% 7032| 50 ae 
Mo Tu 72.002 17.0.2527 15.027572 50.06 10,02 vn A Soe Toes ro 105 en 
Albedo = 0.00 
1.00 6597| .8884| .9301] .9709| .9901 -9958 .0991° | .0998 .9999 .9999 .9999 | 1.0000 
0.60 3250] .4961| .5332|) .5714| .5902| .5958 | .5992 | .5998 | .5999 | .5999 | .5999 | .6000 
0.10 0296| .0497| .0594| .0772| .0909| .09595| .0992 | .0998 | .0999 | .0999 | .0999 .1000 
0.02 0050| .0078| .0086| .o100| .0137| .0166 .O192 .0197°| .0199 .0199 .O199 .0200 
Albedo = 0.25 
1.00 7427 9302| .9587| .9840| .9948| .9979 .9996 .9999 .9999 | 1,0000 | I,0000 | 1.0000 
0.60 3658] .5194| .5495| .5792| .5930| .5971 .5994 .5998 .5999 .6000 .6000 .6000 
0.10 0334| .0520| .0612| .0783| .0914| .0961 .0992 .0998 .0999 .1000 .1000 .1000 
0.02 0056| .0082| .0089| .o10o1| .0138| .0166 | .0192 .0197°| .0199 .0200 .0200 .0200 
Albedo = 0.80 
1.00 | 1.0269] 1.0377| 1.0282| 1.0142] 1.0054| 1.0025 | 1.0005 | 1.0002 | 1.0001 | 1.0000 | 1.0000 | 1.0000 
0.60 5058| .5794| .5894| .5960| .5993| .5999 .6000 | .6000 .6000 | .6000 .6000 .6000 
0.10 0461 0580] .0657| .0807| .0923| .0966 .0993 .0998 .0999 .1000 .1000 .1000 
0.02 -0078| .0091] .0095| .0105| .0130| .0167 | .0192 | .0197 | .o199 | .0200 | .0200 .0200 
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Table IV. The absolute global radiation n(F,G) , in watt per m’ per 100 À as a function of wavelength 
in microns and solar zenith distance ©,, based on Nicolet’s values for nF y, 


7 Ov o° gen 84°.3 poe Où 53.1 84°.3 
2 2.14 1.00 0.10 .52 18.22 10.52 122 
30 2.59 1.25 Ee SS 19.07 11.03 1.31 
31 4.02 1.96 18 .54 19.18 11.12 1.34 
32 4.90 2.44 22 55 18.39 10.71 1.3 
33 6.20 SATA 28 .56 19.20 11.19 1.40 
34 6.35 3.26 29 57 18.99 11.09 1.41 
35 6.02 3.59 32 .58 19.08 De TO 1.43 
.36 7:33 3.85 34 .59 18.53 10.86 I.4I 
37 7.50 3.98 36 .60 18.34 10.76 1.42 
.38 7.56 4.06 37 .65 16.42 9.70 1.35 
.39 8.21 4-44 .4I .70 14.56 3.65 1.26 
40 13.06 7.13 .66 .75 12.87 7.05 1.15 
41 15.02 8.2 .78 .80 11.42 6.80 1.04 

2 15.33 8.47 .82 .85 10.38 6.19 .97 
43 14.47 8.05 80 .90 9.36 5.59 .88 
44 17.01 9.50 .96 .95 8.45 5.05 .8I 
45 18.82 10.56 1.09 1.00 7.65 4.58 .74 
46 19.59 11.04 1.16 1.50 2.98 1.79 .30 
47 19.65 TRS 1.19 2.00 1.19 DE nz) 
48 20.08 11.42 1.2 2.50 47 .28 .05 
49 18.75 10.72 FATO 3.00 .26 .16 .03 
50 19.36 DER AUT 1.25 3.50 .16 .10 .02 
51 19.43 11.20 1.28 4.00 .10 .06 .OI 


same method. As for the albedo correction, 
the interpolation was obtained by means 
of a curve of the quantity As/(1 -As) which 
closely parallels the H; curve for ©,=0°, 
as can be seen in Fig. 1. The physical meaning 
of this quantity is evident from the expression 
(13); when multiplied by (10) it gives the 
additive correction to the sky radiation 
contributed by the ground reflection. 

A sufficiently detailed spectrum of the abso- 
lute global z(F,G): and sky radiation x(F,H);, 
as obtained by using Nicolet’s values for 2Fo,, 
appears in Tables IV and V respectively for 
three solar positions*. It should be emphasized 
that these computations involve the assumption 
of perfect scattering in a Rayleigh atmosphere, 
hence the reader should bear this in mind 
when comparing with any actual observa- 
tions. The irregularities appearing in the region 


of wavelengths shorter than 6,000 A correspond 


to the Fraunhofer absorption regions indi- 
cated by Nicolet’s data. A subtraction of 
corresponding values of Table V from those 


* The units of watt/m? were adopted for the net 
absolute flux following Nicolet. To convert this to 
Langleys or calories/cm?, minute, multiply by the 
conversion factor 14.3335 x 1074. 
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of Table IV yields values for the absolute sun 
radiation z(FyS)». 

The corresponding values including the 
ground reflection with an albedo of 0.25 
(considered rather high for mean planetary re- 
flection) are not shown, but their variation 
with wavelength is similar to that of the sky 
radiation. The relative importance of each 
component of the global radiation can be 
readily seen from the diagram of Fig. 2, in 
which the wavelength scale is logarithmic. 
This shows the spectral distribution of sun, 
sky, and reflected radiation for the solar ze- 
nith distances of 0°, 53°.1 and 84°.3 respec- 
tively, as indicated. The lowermost curve 
in each case represents the absolute sun radia- 
tion, (FS), received through unit hori- 
zontal area at sea level of the earth’s surface. 
This radiation is poor in blue and ultraviolet 
light, and becomes even poorer in this respect 
the closer the sun is to the horizon. The next 
curve above, enclosing the hatched area, 
represents the contribution from the sky 
radiation, 2(FyH),. This shows clearly how 
the light from the entire sky, according to the 
Rayleigh theory, restores, so to speak, most 
of the short wave radiation lost by attenuation 
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Fig. 3. Variation of the spectral distribution of sky radiation with solar zenith distance, based on the extra- 

terrestrial flux curve of Fig. 2 and the Rayleigh-Chandrasekhar theory of multiple scattering. The full curves 

correspond to zero earth reflection and the dash curve shows the total sky radiation for a zenith sun in the case 
of reflection with albedo 0.25. 


in the parallel radiation, so that the global 
radiation reaching the earth’s surface is quite 
similar in spectral distribution (though not in 
magnitude) to the solar extraterrestrial curve. 
This similarity is conserved for a wide range 
of solar zenith distances, i.e. during the greater 
part of daylight hours. The additional sky 
radiation in the case of earth reflection with 
an albedo of 0.25, represented by the 
blackened area on top of the hatched area, 
is seen to be similar in distribution to the sky 
radiation. 

The theoretical spectrum of the absolute 
sky radiation and its variation with solar 
position appears in the diagram of Fig. 3. 
This is drawn on the same basis as that of 
Fig. 2, except that the flux scale has been 
expanded. In the family of curves representing 
the case without earth reflection, the one 
corresponding to the sun at the zenith shows 
two pronounced maxima at about 3,300 and 
4,100 A respectively and a secondary one at 
about 4,500 A, with minima in between which 


look like broad absorption regions. This 
feature of the sky radiation profile results 
from the character of the solar extraterrestrial 
flux curve (uppermost curve of Fig. 2) and 
the relative sky radiation values computed 
in the present work (Table II). It is seen from 
Fig. 3 that the maximum flux of sky radiation 
gradually shifts toward longer wavelengths 
with increasing solar zenith distance, so that 
for O,=84°.3 there is only one broad peak 
around 4,600 An interesting feature is 
that the spectral distribution of the sky radiation 
changes very gradually with increasing solar 
zenith distance, showing a persistent pre- 
ponderance in blue and near ultraviolet 
during most of the daylight hours. 

In order to demonstrate the effect of re- 
flection according to the Lambert law, the 
dash curve is included in Fig. 3, representing 
the spectrum of the sky radiation in the case 
of reflection with an albedo of 0.25, when the 
sun is at the zenith. It is seen that the contribu- 
tion from such an assumed reflection law is 
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quite large, and its effect is to accentuate the 
maxima and minima. 

In all cases under the assumptions of pure 
Rayleigh scattering, there is very little sky 
radiation in red and infrared wavelengths. 
It should also be noted that there is consider- 
able sky radiation in the near ultraviolet 
wavelengths not shown in Fig. 3. The compu- 
tations could not be extended to wavelengths 
shorter than about 3,000 Ä, due to the fact 
that the tables (SEKERA, BLANCH, 1952) of 
the functions y, and y, do not include values 
corresponding to T > 1.00. 


4. Integrated radiations 


It would be quite interesting to compare 
the theoretical results obtained by the method 
described above with the actual spectrum of 
sky radiation observed at the earth’s surface 
on a very clear day. At least in tbe visible 
range, except perhaps for small differences 
around 6,000 i (Chappuis region of ozone 
absorption), theory and observation should 
check quite well. Unfortunately, due perhaps 
to difficulties in instrument design, few or no 
observations exist providing a sufliciently 
detailed knowledge of monochromatic radia- 
tion from the whole sky. Observations so far 
available, including the more recent ones, 
generally have been concerned with sky and 
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global radiation measured by bolometric 
methods, which means a summation of radi- 
ation fluxes of all wavelengths. Such data 
will allow only a very rough check of the 
theory, and must be compared with the inte- 
grated values of the sun, sky, and global radi- 
ation represented by the expressions (22). In the 
present case, the integration had to be limited 
to the range 0.29 u SAS 4.00 u, valid at sea 
level. This limitation may not effect the re- 
sults appreciably since in the atmosphere 
below the main ozone layer, the solar energy 
at wavelengths shorter than 3,000 Ä is almost 
completely cut off. 

The results of the integration (by the tra- 
pezoidal rule) are given in Table VI, in watt 
per square meter, where the symbols Hp, 
H:5 and H.,, stand for the integrated absolute 
sky radiation for the case of no reflection, and 
reflection with albedo 0.25 and 0.80 respec- 
tively; Go, G..55 and G.,, for the corresponding 
global radiation; and S for the integrated 
sun radiation. Referring to Fig. 2, the quantity 
S represents the areas under the lowermost 
curves, H, that of the hatched regions, and 
H:5 that of the hatched plus black regions. 
For comparison one can use a restricted solar 
constant i.e. the area under the zFy, curve of 
Fig. 2 in the range .29 u SAS 4.00 u, com- 
puted by the same method, which amounts 


Table V. The absolute sky radiation ı(F,H), in watt per m? per 100 Ä on the same basis as Table IV 


u 0° BL 84°.3 | 7 9 0° DSL 84°.3 
.29 1.20 0.78 0.10 50 1.39 DS .076 
.30 32 .90 0.12 51 120 722 TB 
Say: mo MRS 18 52 1.12 1.06 .66 
132 1.99 I.49 22 53 1.08 1.03 66 
.33 2.25 1.75 28 54 I.OI .06 64 
-34 2.06 1.65 2 55 .89 .86 58 
.35 2.02 1.65 32 56 BO, .83 58 
.36 1.92 1.61 34 57 80 Tl 55 
. 1.77 1.52 35 58 75 .72 53 
33 ao 1.40 36 59 68 66 49 
.39 1.58 1.39 39 .60 63 60 46 
.40 2277 2.04 62 .65 41 40 2 
-41 2.38 2.14 71 .70 26 26 22, 
Ag 222 2.01 Wiz 75 18 18 16 
-43 1.9I 1.74 68 80 12 12 IT 
-44 2.04 1.87 79 85 08 08 08 
45 2.06 1.9I 86 90 06 06 06 
.46 1.98 1.83 88 95 .04 04 04 
.47 ara: 1.69 86 1.00 02 02 02 
.48 1.70 1.60 83 1.50 00 00 00 
.49 1.46 1.38 79 
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Fig. 4. Comparison of the ratio S/H of integrated sun and sky radiation (as a function of air mass or sec©)) . 
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with observed values. The full curves marked A=o and .25 represent the theoretical values for albedo zero 
and 0.25 computed in this paper (see text). 


to 1,390.55 watt/m? or 1.99 cal/cm? min. 
Comparing this with Gy at 9, =0°, i.e. 1,330.7 
watt/m?, one can say that under these condi- 
tions, the sun and sky radiation received at the 
carth’s surface from a pure Rayleigh atmos- 


Table VI. The zntegrated absolute sun, sky and 
global radiation, S$, H, and G respectively, in watt 
per m? over the range .29 4 < À < 4.00, for four 
positions of the Sun. The subscripts 0, .25 and .80 
refer to values corresponding to reflection with 
albedo 0.0, 0.25 and 0.80 respectively. The ratios 
H/G for different albedos are also shown for con- 
venience. The solar constant in the range used is 
1390.55 watt/m? 


O, | où | Sal 84°.3 | 88°.8 
VE LS ee 57.0 50.9 23.8 7.8 
TÉL 78.4 62.9 25.6 8.0 
lee Rebate 134.7 94.5 28.9 8.6 
Se ete 27807 728.6 86.7 10.1 
Go MEL 1,330.7 779.5 110.5 17.9 
Cros ere 1,352.1 791.5 11223 18.1 
G.so 1,408.4 823.1 115.6 18.7 
HolGo 0.0428 | 0.0653 | 0.2156 | 0.4377 
H 95/G 95 0.0579 | 0.0795 | 0.2275 | 0.4440 
ten Gase 0.0956 | 0.1148 | 0.2502 | 0.4595 


phere is about 95 % of the solar constant. In 
the extreme case of reflection with an albedo 
of 0.80, the integrated global radiation re- 
presents 101 % of the solar constant, an 
interesting result accounted for by the use of 
Lambert’s law, already discussed earlier. The 
integrated global radiation drops off with 
increasing zenith distance, slowly at first, and 
more rapidly as ©, increases. The dependance 
is such that log G, varies almost linearly with 
log uo- 

From the values of Table VI, one can also 
appreciate the relative magnitude of the inte- 
grated sky and sun radiation, and how this 
varies with zenith distance. When the sun is 
at the zenith, Hy is only 4 per cent while S is 
96 per cent of Gy; but when 0,=84°.3 
the percentages are 214 and 784 respectively. 
Table VI also includes the ratios H/G. 

A good set of early observations of the 
integrated sky radiation or ‘sky brightness’ 
were made by Moore and ABBOT (1920) of 
the Smithsonian Institution in 1917, at Hump 
Mountain, North Carolina, 1462 meters above 
sea level. Values of net integrated fluxes of 
sun and sky radiation in Langleys as well as 
their ratio S/H, are given as a function of 
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sun elevation. From an inspection of the 
values of Table VI it is seen that this ratic, 
should drop with increasing solar zenith 
distance. Also from physical considerations, 
S/H should increase for all solar positions, if 
the normal optical thickness is decreased, i.e. 
if the observer is higher in the atmosphere. 
The nature of this ratio is seen better from the 
diagram in Fig. 4, in which the ordinates 
represent the ratio S/H (logarithmic scale) 
and the abscissae are marked in sec ©, or the 
‘air mass’. The two top smooth curves marked 
A=o and A=.25 represent the theoretical 
variation of S/H) and S/H; respectively. 
The five dash curves marked with the dates 
August 28, September 12, October 14, 15, 
and 28 result by merely drawing straight lines 
between the points corresponding to some of 
Moore and Abbot’s observations (1920). In 
general the order of magnitude and slope of 
the observational data agree with the theory, 
however, most of the observed values are 
lower than the theoretical ones, contrary to 
expectations (i.e. considering the height of 
the observatory above sea level). One possible 
explanation would be that the atmosphere 
above Hump Mountain contained a large 
number of non Rayleigh scatterers during the 
days of observation. 

This discrepancy with the theory appears 
even more clearly from the curve marked 
1913, which represents two days’ observations 
made during that year by the Smithsonian 
Institution at Mt. Wilson (cf. ABBOT ET AL., 
1932, page 274), at an elevation of 1,737 
meters. Since Mt. Wilson is higher than Hump 
Mountain, the corresponding values of S/H 
should be larger than at Hump Mountain. 
Note, however, that the opposite was ob- 
served. The explanation might be connected 
with the fact that the volcanic dust thrown 
into the upper atmosphere during the June 
1912 eruption of Mt. Katmai in Alaska, was 
still visible in California in September 1913, 


‚as pointed out by the Smithsonian authors 


(cf. loc. cit., p. 268). 

One more observation is shown by the 
curve marked “NIC”, representing best condi- 
tions as observed by NICOLET and DOGNIAUX 
(1951)* at Uccle, Belgium, practically at sea 

* Nicolet furnishes values of the ratio H/G (cf. Nrco- 


LET, DOGNIAUX, 1951, p. 13) from which the S/H 
values can be deduced. 
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level. This falls well below the theoretical 
curve for sea level, and furthermore the slope 
shows some disagreement with that of the 
theoretical curve. This might be connected 
with a systematic dependence of the turbidity 
of the atmosphere above Uccle, on the solar 
zenith distance. Note, however, that the 
volcanic dust condition at Mt. Wilson is 
more effective than the turbidity near sea 
level in reducing the values of S/H. 


5. Discussion 


Can one make any definite statements 
about the relation between the observed ratio 
S/H and the turbidity of the actual atmos- 
phere? We do not believe one can, and 
disagree with Moore and ABBOT’s conclusions 
(1920), to the effect that turbidity or the 
existence of large particles in the atmosphere 
necessarily results in a lowering of the ratio 
S/H. This quantity in itself is not sufficient to 
tell us anything about the causes of deviation 
from the theory, since it represents the ratio of 
two integrals, both strongly dependent on 
the scattering mechanisms going on in the 
real atmosphere. Perhaps a more definite 
statement could be made with respect to 
volcanic dust in the high atmosphere: The 
opaqueness of this aerosol was demonstrated 
by the marked increase in normal optical 
thickness of the atmosphere in the years 
immediately following the Katmai eruption. 
Comparing the corresponding curve (marked 
1913) and the theoretical one for zero albedo 
in Fig. 4, one notes that they are remarkably 
parallel, and since the vertical scale is loga- 
rithmic this indicates a constant factor of 
proportionality, about 0.30 in this case. If one 
assumes that the sky radiation was almost 
unaffected, this suggests that the volcanic dust 
acted as a neutral filter in the high atmosphere, 
absorbing radiation uniformly at all wave- 
lengths. This of course is a highly tentative 
conclusion. 

In general, a study of the ratio S/H as 
observed at different places shows remarkable 
agreement with the theory, considering the 
very restrictive assumptions made in the 
latter. 

The quantitative results of this paper can be 
compared with the results obtained by F. 
BERNHARDT (1952, 1953). In his two papers, 
Bernhardt attempts to evaluate the effect of 
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multiple scattering by the procedure suggested 
by F. Linke (1942), and based on. the assump- 
tions (a) that the effect of higher order scatter- 
ing is additive to that of the primary scatter- 
ing, and (b) that the ratio of the increase in 
radiation due to the k-th order of scattering 
to the increase due to the (k-r) th order of 
scattering, is constant. If I® denotes the 
specific intensity of diffuse sky radiation due 
to primary scattering, and I® its increase 
due to the k-th order scattering, then the 
specific intensity including the effect of all 
orders of scattering can be written, according 
to the above assumption, in the form 


[=E1® =I (1-)-", 


k=1 
I 18 Te) 
where q= TOR (CT (ee constant. 


Now in the exact theory which was used in 
our discussion, the effects of the higher order 
scattering are included by means of successive 
iterations of four pairs of simultaneous, non- 
linear integral equations, defining the X- 
and Y-functions, mentioned under equation 
(15) above. The relation between successive 
orders of scattering is obviously not as simple 
as indicated in assumptions (a) and (b) above; 
this relation can be seen clearly from the 
values of successive approximations of the X- 
and Y-functions appearing under Table IX 
(SEKERA, BLANCH, 1952). Furthermore, in the 
computation of intensities due to the primary 
scattering, by Hess and LINKE (1942), and in 
the computation of the secondary scattering 
by BERNHARDT (1952, 1953), as well as in 
the earlier work of Kine (1913), the polari- 
zation of the scattered light has been com- 
pletely ignored and this leads to non-negli- 
gible differences as shown by CHANDRA- 
SEKHAR (1950, p. 264; cf. also VAN DE Hursr 
[1949] p- 76, 77). Finally, the computations 
made by Hess and Linke, and as extended by 
Bernhardt, were performed with the use of 
mean values of the transmission coefficient g:* 
and of the extraterrestrial flux xF, in the 
three spectral ranges: À < 0.4, 0.4 £ À < 0.6, 
0.6 E21 =300. 


* A quantity corresponding to exp (-T) in our nota- 
tion. 
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It is thus quite difficult to interpret properly 
the differences between Bernhardt’s and our 
results. In order to eliminate the differences 
resulting from the use of a different extra- 
terrestrial Aux, we had to deduce a mean 


relative sky radiation H; from our monochro- 
matic values covering the above ranges, 
which was then multiplied by the correspond- 
ing extraterrestrial fluxes used by both Hess 
& Linke and Bernhardt. The resulting abso- 
lute fluxes of sky radiation are compared with 
Bernhardt’s values in Fig. 5, which shows 
their variation with solar zenith distance Oy. 
The largest systematic difference appears in 
the curves marked 3, where Bernhardt’s values 
exceed ours by about 60 per cent for all solar 
zenith distances shown. Now in this range 
the third order scattering is quite negligible, 
hence the corresponding values of the ratio 
H}”/ Hy; appearing under section 3 above really 
represent the ratio Hj”/(Hj?+ Hi”) from 
which the ratio Hj”’/(H}” + H3”) can be easily 


deduced. This is about 0.11 for the beginning À 


of the range (A=0.809 u) and 0.04 for the end 


Exact theory : 
as ae Bernhardt 

I 290<A<.400u 
2- 400s \<.600pn 
3- 600 <X<.900u 
4- Total 


m col/cm?, min. 


70° 


80° 


90° 


Fig. 5. Variation of mean sky radiation for the three 

spectral ranges shown, with solar zenith distance. Full 

lines represent the exact theory, dash lines the approxi- 
mation carried out by Bernhardt (see text). 
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of the range (A=0.618 u), with a mean value 
of 0.075, comparable with Bernhardt’s value 
of 0.077. Similarly the transmission varies in 
the same range from 0.94 to 0.98, with a mean 
of 0.96 as compared with Bernhardt’s value of 
0.95. It is clear therefore that the difference in 
the curves 3 of Fig. 5 can only be accounted 
for by a difference in the respective sky radi- 
ation due to primary scattering. 

This can be easily verified by evaluating 
mean values of H;'’ from our computations, 
for the same three spectral ranges, and com- 
paring the absolute fluxes obtained with those 
given by Bernhardt as corresponding to pri- 
mary scattering. Our values, obtained by 
graphical methods of finding the mean Hi, 
are compared with those of Bernhardt below: 


Spectral range | A<o.4 |o.4SA<0.6] 0.6</<o.9 
Absolute fluxes | 
Big tate) 15.33 32.91 8.48 mcal/ 
cm? min 
Bernhardt’s 
values...... 16.48 34.30 13.27 mcal/ 
cm? min 
Differencein%| 7.5 % 420060150510 


The values shown in the last column above 
_ corroborate our statement that the large differ- 
ence occuring in the curves 3 of Fig. 5 is due 
to a different estimate of the sky radiation 
from primary scattering. Since our values 
for Hj" in the third range do not show any 
systematic deviation and are continuous when 
plotted next to the values of H, for shorter 
wavelengths (cf. Fig. 1); and since these Hj"’ 
values approach asymptotically the theoretical 
value when To, we are inclined to believe 
that the values used by Bernhardt for the 
primary scattering are in error, at least in 
the third range. 

_ Another discrepancy is evident if one con- 
_ siders that Bernhardt’s values in the first two 
columns above also exceed ours somewhat, 
however the corresponding sky radiation 
due to all orders of scattering computed by 
- Bernhardt’s method is lower than ours for 
small solar zenith distance and exceeds our 
values for O, >50°, as seen from the curves 
marked ı and 2 in Fig. 5. Thus it would 
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seem that the evaluation of multiple scattering 
eftects from those of lower order scattering 
by a simple assumption such as mentioned 
above, and the omission of the polarization 
of the scattered radiation lead to systeme“. 
deviations from the correct theoretical value: 
the method results in underestimating the 
sky radiation for high sun’s elevations and in 
overestimating this radiation for low sun 
elevations. 

Finally, with respect to actual measure- 
ments of monochromatic sky radiation, one 
of the rare works of this kind available is 
that carried out by Perrir (1932) in the late 
twenties. He was mainly interested in finding 
the ratio of ultraviolet to green direct solar 
radiation received at the ground, and in the 
determination of the solar extraterrestrial flux 
in the short wavelengths. From this work 
the only information which could be used 
here for comparison is his scanty measure- 
ments of sky radiation at .320 u. For Pasa- 
dena, California (elevation 258 m.) the ratio 
H;/S1 at this wavelength is given as 1.00 
but unfortunately the exact solar altitude is 
not indicated, except that the sun was near 
the zenith. From our Tables IV and V, the 
ratio at sea level for this wavelength at ©,=o° 
turns out to be 1.99/2.91 or about 0.66. How- 
ever, if one considers that this ratio increases 
rapidly with ©,, and should increase also with 
the increase of atmospheric turbidity, Pettit’s 
value seems very reasonable. 

Pettit also shows the variation of sky radia- 
tion at .320 u with solar altitude (cf. Fig. 6 in 
PETTIT, 1932) and finds a simple linear de- 
pendence. This does not quite check with the 
theory and the linear extrapolation to a zenith 
sun is not justified. The slope of the corre- 
sponding theoretical curve, if plotted on the 
same basis, would decrease gradually becoming 
horizontal for a zenith sun, as should be 
expected from physical considerations. The 
seemingly straight character of the plotted 
points might be due to the gradual increase 
in haze mentioned by the author, toward 
noon and early afternoon, which makes the 
afternoon plots fall above the forenoon ones. 

A more complete and carefully made set 
of such observations giving the ratio H;/S, 
in narrow spectral bands would be highly 
desirable. It would be useful not only in 
checking the Rayleigh model for a pure 
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atmosphere, but also in learning more about 
large particles existing in the real atmosphere, 
since the ratio H:/S; seems to be highly 
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sensitive both to variations in wavelength 
and the size and character of the large particles 
present. 
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A Theory of the Ventilated Psychrometer Based on the 
Thermodynamics of Irreversible Processes 


By D. J. BOUMAN, Royal Netherlands Meteorological Institute 


(Manuscript received April 23, 1954) 


Abstract 


It is shown that the thermodynamics of irreversible processes as developed by PRIGOGINE 
(1947) and DE GROOT (1951) offer the possibility of deriving the psychrometric equation 


without using pseudothermostatic reasoning. 


I. The traditional theory 


As a specimen of the theory of the ventilated 
psychrometer one may quote MONTGOMERY 
(1948). In this paper the argument runs liter- 
ally (except for some slight modifications in 
notation): 

“It follows from the definition of wet- 
bulb temperature that a mass M of dry air 
and a mass Mx of water vapor are imagined 
to cool from T to T,, and a mass M(x,,— x) of 
water to evaporate at the temperature T,,. 
The enthalpy of the dry air decreases by Mc,, 

T-T,) and that of water vapor decreases 
by Mxc,,(T- Ty). The enthalpy of the added 
water increases by ML,(x,-x). Since the 
process is adiabatic and isobaric, the total 
change in enthalpy is zero. Hence 


(Co + X€p,) (T- Tw) =Lul&w- x) I (x) 


the basic and complete form of the psychro- 
metric formula (for a ventilated psychrometer).”’ 

Others have given derivations of (1) which 
differ only in minor points from that of 
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MONTGOMERY, their main difference being 
the substitution of the word heat for enthalpy 
which illustrates the loose use often found in 
meteorological literature of these words which, 
however, denote quite different thermody- 
namical concepts. . 

At first sight there is nothing wrong with 
arguments of this kind. But a little reflection 
shows that the whole way of reasoning is at 
variance with the basic postulates of ther- 
modynamics. For it is known that classical 
thermodynamics are only valid under very 
stringent restrictions. Indeed the application 
of thermodynamics is limited to: 


a) closed systems 

b) equilibrium conditions 
c) quasi-static processes 
d) reversible processes 


Ad a. The system under consideration is clear- 
ly not closed. Therefore MonTGOMERY’s 
remark “the total change in enthalpy 
is zero” is not valid in this form. For 
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one never knows how much enthalpy 
is added to the system by the ventila- 
tion. 


Ad b. The system under consideration lacks 
the most essential condition of equili- 
brium viz. the uniformity of tempera- 
ture through the system. Even the final 
state is not an equilibrium but a sta- 
tionary state. But classical thermody- 
namics offer no tools to treat stationary 


States. 


Ad c. The process which occurs in a ventilated 
psychrometer cannot be considered as a 
quasi-static process, since it is known 
that it is impossible to ventilate in- 
finitesimally slowly. In fact the ven- 
tilation should preferably exceed 3 or 


ea SCC. 


Ad. d. As it is impossible to apply a reversed 
ventilation the process cannot be con- 


sidered as reversible. 


Of course there is nothing to say against 
MONTGOMERY’S reasoning if it is interpreted 
not as a theory of the ventilated psychrometer 
but as a theory of the concept of wet-bulb 
temperature as a property of the air. The real 
problem is to prove the identity of a theoreti- 
cally defined wet-bulb temperature with the 
temperature of the wet-bulb thermometer. 

Some authors have given treatments in 
which they try to avoid one of the above 
mentioned four defects. The most important 
paper on this subject is one by van MIEGHEM 
(1941), who used the thermodynamics of 
open systems. But a theory in which the four 
defects are avoided simultaneously has as far 
as the present author knows never been given. 
The main purpose of this paper is to offer 
such a treatment with the aid of the recently 
developed thermodynamics of irreversible 
processes. A short summary of this new theory 
will be given in the next section. For details 
the reader should consult the books of Pri- 
GOGINE (1947) and DE Groot (1951), where the 
subject is dealt with in full detail. 


2. Introductory remarks on the thermody- 
namics of irreversible processes 


In the thermodynamics of open systems the 


differential d2 of any quantity Q is split up 
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in two parts, the exterior part d.Q and the 
interior part dQ. 


dQ=d,Q+d4,Q (2) 


The exterior part is defined as the change in 
Q as far as it is a consequence of the exchange 
of matter with the surroundings, the interior 
part is defined by equation (2). By definition 
a pseudo-adiabatic process is a process for 


which 
d;Q=0 (3) 


The last equation might be interpreted as 
follows: except for what happens as a conse- 
quence of exchange of matter with the 
surroundings, the system under consideration 
exchanges no heat with the surroundings. 
Strictly speaking there is no first law for 
open systems. The often quoted formula* 


dQ=dE+PdV-hdM=dH-VdP-hdM (4) 
is not a generalisation of the first law for 
closed systems to open systems, but is to be 
considered as a definition of dQ, a quantity 
which for open systems is often not measur- 
able. In fact other definitions may be used e.g. 
the definition proposed by TOLHOEK and DE 
GROOT (1952) ‘ 


dQ* =dE+PdV-Sh,d.M, (5) 

k 
Comparing (4) with (5) one infers that the 
ambiguity in the definition of dQ has no 


consequences for the equation of the pseudo- 
adiabatic process, which reads in both cases 


d;Q = d;Q* =d;E+Pd;V=o (6) 
The second law for open systems offers no 
difficulties and has already been given by 
GiBBs (1875— 78) 
TdS= dE + PdV—ZuydMy =dH- VdP - SupdMy 
k k | 
(7) 
As far the entropyproduction o, defined by 


o= ea (8) 


* List of symbols in section 5. 
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classical thermodynamics can only assert 
520 (9) 


but according to the new theory o can always 
be written as 
(10) 


The quantities J4, which must be time-deri- 
vatives of some function of state, are called 
flows or fluxes. It appears that they are flows of 
matter, electrical currents, flows of heat, etc. 
The quantities X4 are called conjugated forces 
or conjugated affinities. Apart from some factor 
they are temperature gradients, electrical 
potential gradients, concentration gradients, etc. 
In the case of the presence of chemical reac- 
tions one of the X, is the affinity as it is ex- 
tensively used by the Belgian thermodynam- 
icists of the school of DE Donper. 

For a large class of problems one may 
assume linear relations between fluxes and 
forces, called phenomenological relations: 


o=L] 4X4 
A 


(11) 


The coefficients L4p are called phenomenological 
coefficients; in most cases the only way to 
compute them is by kinetic reasoning. They 
appear to be of the same kind as e.g. diffusion 
coefficients. 

In 1931 Onsacer proved with the tools of 
statistical mechanics his famous reciprocity 


law 
(12) 


which later on (1945) was discussed by Casr- 
MIR. 
At first sight this formalism does not seem 
to be of much use. But PRIGOGINE (1947) and 
DE Groot (1951) have shown in their books 
that the reciprocity law (12) is a mighty 
weapon which in many cases enables one 
to arrive at far-reaching conclusions without 
the knowledge of the exact values of the 


Ja en 


Lap=Lpa 


' coefficients L 4p. 


3. The ventilated psychrometer 


The psychrometer is a polythermal system, 
consisting of two phases, a gaseous phase and a 
liquid phase, having different temperatures. 
The following convention will be used: 
quantities referring to the gaseous phase are 
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denoted by symbols with one accent, quantitie 
referring to the liquid phase by symbols witl 
two accents, whereas quantities referring to 
the whole system are denoted by symbol: 
without accents. With the subscripts 1, 2 
and 3 the components dry air, water vapor and 
liquid water are indicated, thus e.g. M,, stands 
for the mass of water vapor in the gaseous 
phase. 


Furthermore the abbreviation 


aM, dM; 
= in ae 


(13) 


is used. Equation (13) expresses the conserva- 
tion of matter. 

As for the distinction between exterior and 
interior changes, it should be kept in mind 
that this distinction refers always to the whole 
system even in an equation which is set up 
for a part of the system. 

Now apply the first “law” (4) first to the 
whole system and afterwards to the two 
subsystems formed by the gaseous and the 
liquid phase. Since the process is pseudo- 
adiabatic and isobaric, one gets 


DOME 
ve ee (14) 
LOU Hee RL, 
dt dt ate dt dt zu (5, 
d; Q’ d;H é 7 d;M x d;H f 7 
ae 


Equation (14) may be considered as the exact 
formulation of MONTGOMERY’s statement “the 
total enthalpy does not change”. One observes 
that this is only valid as far as exterior changes 
are left out of consideration. Adding (15) 
and (16) together one gets 


d;Q’ 4,0% d,H’ | d;H" 
bi dt dt 


(hh) Jy 


from which together with (14) and 


H=H'+H" 
it follows that 
d} j d; hi 7 2 
EL (th), (m) 
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an equation which is often called Prıcocıne’s 
equation since PRIGOGINE discussed it exten- 
sively in his book (PRIGOGINE, 1951). As a 
consequence of (17) one gets 
d;Q' +4,Q” #diQ (18) 

an inequality which is one of the often un- 
expected freaks of the thermodynamics of 
open systems. 

Now apply the second law (7) to both sub- 
systems 


Sd; EL : 
a: 
NS Pe 
1 dt > dt EU Je 
(19) 
d;S d;S’ d;S" li d;H' U d;H" 
on ea ne 
Ju, us 
u: 


Using (14) and the abbreviations 


d;H' I 
= 8); Sn (21)s 
el, 
= TF eT (22) 


the entropyproduction o may be written as 


o=hXı+JeXg (23) 
which has the form (10). 

It is interesting to remark that van MIEGHEM 
and Durour (1948) in their joint paper 
derived an equation (VII, 42) which is sub- 
stantially equivalent to (19). But since they 
did not use the Onsacer theorem (12) but 
only the inequality (9) they failed to produce 
the psychrometric equation and could only 
prove 


Dern (24) 


in which Tz denotes the dewpoint. The 
fertility of the reciprocity law (12) is thus 
clearly demonstrated by the fact that on 
account of (12) we shall now proceed to deduce 
the said equation. 

The phenomenological equations (IT) are 


Dre ]E2 BEOIUIMENIN 


Jı=Lufı +L12X3 | (25) 


J2= Lan Xi + Los Xe 


and the reciprocity law (12) yields 


Lis = Lan 


(26) 


Therefore 


o=L,X%:+21,X,%,+1,%,; (27) 
It is a well known theorem in classical thermo- 
dynamics, that the state of equilibrium cor- 
responds with maximum entropy. In the ther- 
modynamics of irreversible processes there is 
an analogeous theorem which states that the 
stationary state corresponds with minimum 
entropy production. Since T” is given by the 
exterior circumstances and T” is sought, the 
stationary state is given by 


do 
aT" =0 (28) 


Application to (27) gives, making use of (25) 


and (26) 


do dX 
OX. IX, 
+ (2 Ie 121%) Jr" 2h apt 


OX. 2 af REN 
pe Jh, Js=0 


+2 Ja DT" 
or 


Ji=h Ja (29) 
(Here use is made of the well known formula 


oe ee 7) 


N 


Going back to the definitions (13) and (20) 
one gets from (29) 
d;H’ 


he 
= 


dt 


d;M, 
a (30) 


Now H’ is a function of T’ and the masses 
M, and M, of the gaseous phase. Therefore 


dH! OH' dT’ OH’ dM’ OH’ dM: 
de oT & ‘aM. dt OM. ode 
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By definition 


oH ae » 
am (32) and Zm=G (33) 
whereas 

d;M! 

En, (34) 


Substituting (31, 32, 33, 34) in (30) the last 
equation becomes 

U ER di, 

en 


Now the temperatures T’ and T” are suffi- 
ciently close to each other to substitute in (35) 


u UE 


L being the latent heat of evaporation at some 
suitable intermediate temperature. Further- 
more 


= Me (37) 
and 
d;M, = ’ dx’ 
dt SM dt (38) 


Substitution of (36), (37) and (38) in (35) 
yields: 

rl LE 

ner (39) 
Now one may integrate the last equation 
over a time t in which a quantum of air is in 
thermal contact with the wet surface of the 
psychrometer. During this time 1 the air 
would in absence of exchange of matter with 
the surroundings cool from T’ to T” and 
increase its water content from x’ to xy. 
Thus one arrives at the psychrometric equation 


Ebel L(g =X) (40) 


which is equivalent to (1). 


. 4. Concluding remarks 


A. Montcomery (1948) thinks it is very 
important that in the psychrometric equa- 
tion (1) or (40) one should take the latent 
heat of evaporation at the wet-bulb tem- 
perature. The present theory shows that 
this is not the case. This was already known 
to VAN MIEGHEM (1941) who in his paper 
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clearly demonstrated the approximations 
involved in equation (40) even if it is 
interpreted as an equation for the theo- 
retically defined “wet-bulb temperature”. 
This paper was apparently known to but 
not fully appreciated by MONTGOMERY. 


B. The general theory of irreversible processes 
shows that the theories criticized in par. 
I in many cases may give the correct 
results. But no general rule can be given 
which enables one to decide apriori in 
which circumstances the pseudo-thermo- 
static theories, as they are called nowadays, 
give correct results. Thus one is never 
apriori certain that one may use pseudo- 
thermostatic reasoning unpenalized. The 
present paper shows that the formalism 
developed by PRIGOGINE and DE Groot 
makes the use of pseudo-thermostatic 
methods unnecessary and provides very 
powerful tools, which are easy to handle 
as well as a sound basis for theoretical 
reasoning. 


5. List of symbols 
NB.: Specific means always: per unit of mass. 


Cp = heat capacity of total mass in gaseous 
phase at constant pressure. 


E = internal energy 
= tlowsor ax 
LAB = phenomenological coefficients 
L = latent heat of evaporation 
M’ = total mass gaseous phase 
M: = mass of dry air 
M! = mass of water vapor 
M = mass of liquid water 
P = pressure (overall pressure) 
dQ = absorbed heat 
S = entropy 
T = absolute temperature 
Tz = dewpoint 
V = volume 
X = force or affinity 


c, = specific heat of gaseous phase at con- 
stant pressure 


d, = exterior change 
d; = interior change 


h = mean specific enthalpy 


h;, = specific partial enthalpy of water vapor 
h = specific enthalpy of liquid water 
t = time 


x’ = specific humidity 
x, = maximum specific humidity at wet- 
bulb temperature 


,, dy = chemical potential (specific partial Subscript 1 = dry air 
thermodynamic potential, specific 2 = water vapor 
partial GiBBs function, specific par- 3 = liquid water 
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The Finnish Radiosonde With Rolled Hair as hygrometer 


By LARS RAAB, Swedish Meteorological and Hydrological Institute; Stockholm 


(Manuscript received II June 1954) 


In the last years the rolled hair (also called 
the Frankenberger-hair or Pernix-hair) with 
very good response at low temperatures has 
been used at the humidity sensing element 
in the German morse-radiosonde. In two 
papers SCHULZE (1952) shows from measure- 
ments made at Flensburg how extremely use- 
ful the new method is in measuring humidity. 
The time-constant of the rolled hair hygro- 
meter at normal temperatures is five times 
smaller than that of the unrolled hair and the 
hair is also functioning at as low temperatures 
as —70° C, where the time constant is about 
500 sec at 100 per cent humidity. 

Using the values of the time-constant ob- 
tained by Schulze the time-lag corrections are 
directly applicable at temperatures above 
—20° C, and with some computational work 
applicable also below - 20° C. One disadvan- 
tage, according to Schulze, is that the “drying” 
effect of the hair causes a rise in the time- 
constant which depends on the time of the 
radiosonde passing a very dry layer in the 
atmosphere. As yet no simple method exists 
of applying the corrections in these cases. 

Not all radiosondes are suitable for applying 
the rolled hair in the hygrometer. Especially 
the audiofrequency modulated types where the 
representation of the meteorological elements 
into electrical units are made by means of a 
resistance are as yet not suitable without using 
complicated mechanical constructions. The 
simplest mechanical mounting of the hair is 
obtained when mounted to move one side 
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of a small plate condensor elastically suspend- 
ed. If the condensor is in the tank circuit of an 
radio-frequency oscillator, radio-frequency 
modulation is obtained by means of the hair 
and in relation to the humidity changes. 
This is made in the Finnish sonde. The re- 
duced strength of the rolled hair is of no im- 
portance to the reliability of the measure- 
ments. Breaks in the hair are avoided when 
the sharp bends of the mounted hair are fixed 
with some cellulose paint. 

The Finnish sonde is easily adopted for 
double humidity measurements. The figures 
show the result of a sounding made at Stock- 
holm - Bromma, with two hygrometers, one 
with normal hair and one with rolled hair. 
The humidity are to the left in the figures 
(heavy line for rolled hair, thin line for nor- 
mal hair). At the ground a strong inversion 
keeps the humidity high and equal to 80 per 
cent as indicated by both the hygrometers 
(Point 1). Half the way up in the inversion 
the humidity is very rapidly sinking to low 
values indicating a subsidence layer reaching 
up to 700 mb. The sounding is made in a high 
pressure center with cold air mowing south 
and penetrating the lowest layers. With the 


normal hair this subsidence process at low 


temperature is not clearly shown. 

Above 700 mb, humid air from south is 
moving northwards over western Scandinavia 
and according to the rolled hair, humidities 
above 60 per cent are shown in all layers up 
to the tropopause. Here the normal hair gives 
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Fig. ta and ıb: Sounding with double hygrometers. 
Temperature curves to the right, humidity curves to the 
left (heavy line rolled hair, thin line normal hair). The 
descent curves are indicated with dashed lines. 


no information. Interesting is the marked 
inversion at 550 mb. The rolled hair rapidly 
reaches high humidities (95—100 per cent) 
when passing an observed altocumulus layer. 
From point 12 the humidity decreases but the 
normal hair is still indicating rising humidity 
(point 12—14). The normal hair in fact gives 
a false sign of the humidity lapse-rate. 

In the “stratosphere both ascent and descent 
(dashed lines) curves are shown. At point 
18 the sonde enters the dry stratosphere with a 
humidity according to the descent curve being 
equal to 20 per cent. Between points 18 and 
19 the humidity decreases at a slow rate but 
more rapidly between 19 and 20. This be- 
cause of the increasing temperature with in- 
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crease of the rate of response. Between points 
22 and 23 the humidity stops falling indicating 
a somewhat higher humidity value (30 per 
cent) higher up in the stratosphere (!) than 
in the bottom layer just above the tropopause. 
The very low humidity values in the lower 
stratosphere may be the result of subsidence 
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in the meridional northward flow of air 
from the tropical tropopause levels. 

Even if the correction methods are as yet 
not always useful for the stratospherical parts 
of the sounding, useful conclusions can be 
drawn from both the ascent and descent cur- 
ves in this region. 
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The Density Ratio at the Boundary of the Earth’s Core 


By MARKUS BÄTH, Meteorological Institute, University of Uppsala 


(Manuscript received 8 April 1954) 


Abstract 


A theoretical study is made of the behaviour of a P wave incident from the mantle on the 
boundary of the earth’s core, with special reference to phase reversals. The results are given 
graphically in Fig. 2, which gives the conditions for given values of the angle of emergence e 


r 


and the density ratio o’/o or the ratio of the elastic parameters at the core boundary. 


L 


The dependence of the phase changes on the density ratio is clear; for instance for 9’/9=1 (no 
density discontinuity) there is a change of phase both for e<10° and for e>55°, whereas for 
o’/o > 1.7 there is a change of phase only for e<6.2°. Observations of PcP waves around epi- 
central distance = 30° (Fig. 3) prove conclusively that there is a first-order density discontinuity 


2 


at the core boundary. They also show that 


>0.34. À method for determining the density 


+2 


ratio by more complete observations of PcP is indicated. 


I. Introduction 


It is generally assumed that the density has a 
first-order discontinuity at the surface of the 
earth’s core at a depth of 2,900 km. See e.g. 
BULLEN (1947), pp. 212—219. There are no 
direct data concerning the density jumps at 
the discontinuities within the earth, of which 
the outer core boundary is the most important. 
BULLEN (1947) computes the densities within 
the earth on two different hypotheses, one 
giving a density jump from 5.66 g/cm? to 
9.7 g/cm, i.e. a density ratio of 1.71, at the 
core boundary; the other giving densities of 
5.72 g/cm? and 9.1 g/cm’, 1.e. a ratio of 1.59. 
Later these figures have only been slightly 
changed as consequences of other hypotheses. 


In 1941 KUHN and Ritrmann published 
their hypothesis about the earth’s interior, 


which entailed a density curve with no dis- 
continuity whatsoever at the boundary of 
the core. A sharp change of density at the 
core boundary seems, however, to be most 
probable. GuTENBERG (1951) says that “no 
hypothesis can be considered a good approxi- 
mation which disagrees with the result that 
the boundary of the outer core is very sharp.” 
“Waves having lengths of ten km or even 
less are reflected at both sides of the boundary.” 
During the last years new ideas have been 
proposed about the earth’s interior, especially 
the nature of the boundaries of the two cores. 
I am not going into detail about these things 
here, as I have done that in another paper 
(see References). It may only be added here 
that Professor A. RITTMANN in a lecture in 
Uppsala on September 12, 1953, about the 
KUHN-RITTMANN hypothesis admitted the 
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possibility of a first-order discontinuity of the 
density at the outer core boundary as a conse- 
quence of a phase transition in silicate under 
pressure. 

In disentangling the problem of the density 
behaviour at the core boundary there are at 
least two ways, i.e. investigations of (a) the 
amplitudes of reflected waves, (b) the phase 
changes upon reflexion at the core boundary. 

Amplitudes and travel times of seismic 
waves reflected and refracted at the core 
boundary have already been investigated to a 
large extent. On the other hand, phase changes 
do not seem to have attracted very much 
attention. The reason is natural, because phases 
(compression or dilatation) are generally 
difficult to determine from seismograms, ex- 
cept for the first pulse (P or PKP). But this 
does not seem to be enough reason not to 
investigate the matter, as it may be of some 
theoretical interest anyway. On the whole, 
quite a lot of information could be obtained 
from studies of the direction of initial motion 
of the different seismic waves. 


II. Theory 


The wave functions for the P waves on the 
mantle side of the boundary, i.e. incident and 
reflected P waves, for the refracted P wave, 
and for the reflected S wave (see Fig. 1) are 
respectively 


eg Aeit —4z-wt) = A,ei* +22 — wt) 
D = A'eix = a'z — wt) 
y= Bee + Bz - ot) 


There is still no proof that an S wave exists 
in the outer core. & is connected with the 
angle of emergence e by the relation 


CAUSE 


The notations are the same as used by JEFFREYS 
(1926). Quantities referring to the transmitted 
P wave are denoted by accents, reflected 
‘waves by the suffix 1. The elastic equations of 
motion which must be satisfied by the wave 
functions are 


PE 
0 = (1 + 2u) v?p 

dv 

nn reas! 
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Fig. 1. P wave incident on the core boundary from the 
mantle. 


They give the following relations 


2 
1402 = (2) (1 +02) 


cp 


enabling the calculation of x’ and ß for given 
values of «. cp and cs are the velocities of P 
and S waves just outside the core, and cp 
the velocity of P waves just inside the core. 
o and 0° are the densities. The boundary condi- 
tions (here between a solid and a fluid) are 
that the vertical displacement w and the 
elastic stress components pz, and p>, are con- 
tinuous across the surface z = o. The general 
expressions for these quantities are 


22 
pzx= M pa Sees 
DDR EC" 


This gives us three equations with the ratios 
of the amplitudes A,/A, B,/A, A’/A as the 


unknowns. These equations are respectively 
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nr 


2 
Solving for A,/A we obtain A,/A=1- = 


where 

4x8 | o(B + 1)°x 
(ST ott). 
This means that we have A,/A= o for X=2 
4aBoa' + o'(B%+ 1)? —_ 


a ee 


resp. For given values of the wave velocities 
cp, cp, and cs this relation may be written as 


JCA ae 


This means that the angle of emergence e 
for which f = 1 and A,/A=o depends upon 
the density ratio at the core boundary. 


The curve represented by the equation 


X=1+ 


respectively or for 


f (o'le, e)-1=0 


is given in Fig. 2 (curve aa), assuming cp = 
13.7 km/sec, cp = 8.0 km/sec, and ¢;=7.25 
km/sec. For e = 0° this equation gives an infinite 
value of 0'/o. The limit for e= 90° of the den- 
sity ratio 0.10 115 „equal to} cpp, 1.c ait 
The curve has a minimum with 9’/o0 © 0.4 at 
e= 25°. As e is generally not measured, a 
second curve (bb) is also drawn in Fig. 2, 
giving the relation between e and the angular 
distance A between epicentre and station for a 
surface focus. The relation e — A is determined 
by the equation 
Cp dT 


COS =a 


R dn’ 


where R=the radius of the core (3,470 km) 
dT 
and 7 has been computed from the travel 


times of PcP. For a deep-focus earthquake 
this curve e (A) is not directly applicable but 
needs a correction as a given e then always 
corresponds to a somewhat smaller a. 


MARKUS BATH 


0.5 


Fig. 2. Curve aa represents the equation f (0’/o, e) =1, i.e. 
A,/A =o. 


, 
Curve a’a’ represents the equation fi es: ) = il. 
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The ordinate scale is the same for o’/o as for 


A+2u 


Curve bb gives e (A) for a surface focus. 


2 


Substituting 
°1+2u 


the relation 


for p'/o by means of 


A+2 À 
(cp/cp)? = oe -0'/o 


v 


in the equation 


tele, =a 
this becomes 
CE A el ia eee 
(B21)? A+2u I+a® (pF 1)2«’ 


which may be written as 


Hy 
fi (; +2u’ ) a 


It is given as curve a’a’ in Fig. 2. 

The curves aa and a’a’ in Fig. 2 correspond 
to A,/A=o, the areas above the curves to 
A,/A>o, and the areas below the curves to 
A,/A<o. Naturally only the parts of the 
diagram where o’/o > 1 are of practical 
importance. 

The behaviour of the P wave on reflexion 
at the outside of the core with regard to 
phase changes is immediately clear from 
Fig. 2 for any value of 0’/o and e (or a). If for 
instance there is no first-order density discon- 


Tellus VI (1954), 4 


! 


= At 


Lee DaNSLeYsRATIO AL THE BOUNDARY OF THE EARTHS CORE 


tinuity, Le. e’/o=1, we get a reversal of phase 
OBEN 82) and tones ss" (4230: 
an accurate calculation gives 54.65° as the 
limiting value of e) and there is no reversal 
of phase between these limits. On the other 
hand, assuming a density ratio of 0’/o >1.7, 
there is a phase change only for small e-values 
(e<6.2°), but no phase change for other e- 
values. This is the marked difference between 
this case and cases of smaller density ratio, 
and particularly the case of continuity. 

Two special cases of the general equation 


A,/A=1—2/X 


are of interest. 

(1) e=o°, ie. grazing incidence. In this case 
&=0, X=1, and A,/A= - 1, independent of 
o'/o. This means a total reflexion with reversal 
of phase. 


(2) e= 90°, i.e. normal incidence. In this case 
= ©; 

0" « d 

lim X=1+ > - 


Ce 


Xx=00 


and 


D 


0 
T+— 
0 
cp/cp = 13.7/8.0=1.7 and o'/o=1.7 according 
to one of BuLLEN’s (1947) hypotheses. For these 
values of the ratios, A,/A=o for normal 
incidence, i.e. no PcP wave is obtained. For 
another value of o’Jo we get A,/A 40, but 
most probably rather small. 
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II. Observations and Discussion 


It is often difficult to determine from seismic 
records whether a wave begins as compression 
or as dilatation, especially for the later phases. 
Furthermore the direction of motion of P 
and PcP cannot be directly compared (except 
for grazing incidence), as they start in different 
directions from the hypocentre, and the initial 
motion may be different along different di- 
rections. This difficulty may be overcome 
at least for deep-focus earthquakes by com- 
paring PcP and PcS which start in nearly 
the same direction. 

It may prove helpful to investigate the 
direction of motion of PcP (compression or 
dilatation) for a given shock at a number of 
stations in about the same azimuth but at 
different epicentral distances. If at some 
distance a change compression to dilatation 
or vice versa is found for PcP, this may either 
be due to a change of the initial motion due 
to different direction from the hypocentre, 
or it may be due to a change on the reflexion 
at the core. If for a larger number of shocks, 
investigated in this way in different azimuths, 
there is a change compression to dilatation or 
vice versa for PcP always at the same epi- 
central distance the reason is certainly the 
reflexion at the core, but if the changes occur 
at different distances they must be ascribed to 
the earthquake mechanism. 

In Fig. 3 are reproduced a number of re- 
cords, also given in Table 1, with epicentral 
distances around 30° and with especially 
strong PcP. 

They are partly Central American earth- 
quakes at intermediate depth recorded at 


Table 1. Earthquakes with strong PcP. Reproductions in Fig. 3 


O = origin time, y, À=latitude and longitude of epicentre, H=focal depth (n normal), M = earthquake 


magnitude, A 


geographical epicentral distance. 


à O 4 H : Seismograph 
a | ee | GMT 7 % km zz station A 
I 19 May 1934 | 104737 14 SN | or %W| 120 61% Tinemaha 2 8° 
2 19 May 1934 Pasadena ST TN 
3 20 Jan. 1939 | 204027 13 %Ni or ZW 70 61%, Pasadena 31:8 
4 23 Sept. 1943 150044 15 Nı ot %W| r1o 634 Pasadena 30.75 
5 16 July 1948 | 071230 | 14 %N]| 92 W| 100 64 Pasadena 30,7% 
6 25 Mar: 1952 || 033511 | 34.3 N | 2374 —E n 6 Kiruna 53:04 
7 23 June 1953 | 015312 | 36 N| 25 E| x00 51% Kiruna 31.9° 
8 23 June 1953 Uppsala 24.4 
9 zu, Oe loser | rasyose 38.35 N1E20.8 À n 6% Kiruna Moy 
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Fig. 3a. Records of PcP around A =30°. The numbers on the records refer to Table ı. 
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Fig. 3b. Records of PcP around A =30°. The numbers on the records refer to Table 1. 


Pasadena and Tinemaha (Benioff electro- 30.0° to 29.8°) for the same angle of emer- 
magnetic vertical seismometers), and partly gence (e) at the core as for a surface focus. 
Mediterranean earthquakes recorded at Ki- The depth is therefore of no consequence 
runa and Uppsala (Grenet-Coulomb elec- here. The variations of A,/A with À are con- 
tromagnetic short-period vertical seismome- tinuous (see the general theory), and the slight 
ters). A focal depth of 100 km means only a deviations from À — 30° in our cases are there- 
decrease of 0.2° in epicentral distance (from fore also of no importance. A comparison of 
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direction of first motion of P and PcP would 
require more data as already indicated. But if 
o’/o were=1, then A,/A=o for a = 30°, and 
this would be valid for any P wave incident 
on the core, independent of the incident 
amplitude A and independent of its phase. 
The fact, evident from our observations, 
that PcP may be very strong also for A = 30° is 
conclusive proof that o’/o > rat the core bound- 
ary. We also note that both P and PcP arrive 
as compressions both at Uppsala (a =24.4°) 
and at Kiruna (A=31.9°) for the earthquake 
on 23 June 1953 (No. 7 and 8). There is no 
change of phase of PcP when the distance 
A = 30° is passed. 

It may possibly be objected that PcP is not 
always so strong at A=30° as in the cases 
given here. This is true, as the cases presented 
here have exceptionally strong PcP. But a 
weak or possibly absent PcP at A=30° in 
another case does not invalidate our result, as 
such a circumstance cannot be due to the 
core but must be due to the earthquake 
mechanism. 

From the relation 


tp\* Arzu 0 
cp 40 


it is at least theoretically possible to explain 
the velocity discontinuity cp- cp in three differ- 
ent Ways. 


/ 4/ 


ian À 
(1) ; 1 and er 


been excluded. 


+1. This case has already 


A knowledge of the numerical value of 0’/o 
is required in order to decide between (2) 
and (3). From the study of records for À = 
30° (Fig. 3) and the graphs in Fig. 2 the 
results may be summarized as follows. 

o’/o <1 is excluded, leaving 0’/e>ı as the 
only possibility. This means that A,/A > for 
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A=30°. And this in turn entails that 


A+2 
0.34. According to BULLEN’s hypotheses 0’/o = 


A+ 24 + 

70, he. 5 à 0.59. 

It is true that certain simplifying assumptions 
underlie the theory but it is believed that 
corrections for them would only be of second 
order. 

The study of PcP could be used also for a 
determination of the ratio o’/o. We write 
the general equation for A,/A in the following 


way 
Al 2 ae 0 
iz ze) 


e is a given function of A for a given depth 
of focus. The velocities cp, cs, and cp are 
supposed known. The relation obviously 
means that a determination of A,/A should 
make it possible to compute 0’/o. But A,/A is 
strictly not the same as the amplitude ratio 
PcP/P at a given station, as this would pre- 
suppose equal radiation of energy in all 
directions from the focus. Instead it would 
be necessary first to determine the earthquake 
mechanism, i.e. the orientation of its fault 
plane and the direction of slip along this 
plane, by means of compression and dila- 
tation of P waves or by other means at a 
large number of stations. A, is the observed 
PcP amplitude and A has to be computed 
from observations of P. Another difficulty 
has been observed by MARTNER (1950) and 
ERGIN (1953), namely that PcP may be 
complicated by minor phases; see also ERGIN 


(1952). 


1.59—1.7ı and 


Acknowledgements 


Copies of the records at Pasadena and 
Tinemaha (Fig. 3) were obtained from the 
Seismological Laboratory in Pasadena, Cali- 
fornia. 

Miss E. Wırson, Uppsala, has assisted me in 
drawing the figures. 


414 


MARKUS BÄTH 


REFERENCES 


Batu, M., 1954: Jordkärnan — ett aktuellt geofysiskt 
problem (The Earth’s Core - A Modern Geophysical 
Problem). Föreningen för Filosofi och Specialvetenskap, 
Uppsala, Vol. 3 (in press). 

BULLEN, K. E., 1947: An Introduction to the Theory of 
Seismology. Cambridge Univ. Press, 276 pp. 

ERGIN, K., 1952: Observations on Recorded Ground 
Motion Due to P, PcP, S, and ScS. Bull. Seism. Soc. 
Am., Vol. 42, No. 3, pp. 263 — 270. 


ERGIN, K., 1953: Amplitudes of PcP, PcS, ScS, and 
ScP in Deep-Focus Earthquakes. Bull. Seism. Soc. Am., 
Vol. 43, No. 1, pp. 63 — 83. 


GUTENBERG, B., 1951: PKKP, P’P’, and the Earth’s Core. 
Trans. Am. Geophys. Un., Vol. 32, No. 3, pp. 373 — 
390. 

Jerrreys, H., 1926: The Reflexion and Refraction of 
Elastic Waves. Monthl. Not. Roy. Astr. Soc., Geophys. 
Suppl., Vol. I, No. 7, pp. 321 — 334. 

Kuun, W., and RITTMANN, A., 1941: Uber den Zustand 
des Erdinnern und seine Entstehung aus einem ho- 
mogenen Urzustand. Geol. Rundschau, Bd 32, pp. 
215 — 256. 

MARTNER, S. T., 1950: Observations on Seismic Waves 
Reflected at the Core Boundary of the Earth. Bull. 
Seism. Soc. Am., Vol. 40, No. 2, pp. 95 — 109. 


Tellus VI (1954), 4 


4 


Page Column 


133 

158 2 
159 I 
159 2 
160 i 
162 > 


Tellus VI (1954), 4 


AIS 


ERRATA 


Line Should read 
The legends to fig.6 and 7 have been interchanged 


s from bottom g=— = V2 
3 from top Vig Sev NO EN CNE ENV ic vorge 
4 from top uv=uvt etc. 
22 from top line symmetrical 
7 IY 
8 from bottom ‘ | = — 


LING TERS 


ALFVÉN, H.: On the Origin of Cosmic Radiation, 
232. 

aluminosilicates, accumulation, G. Arrhenius, 215. 

aminoacids, in rainwater, S. Fonselius, 90. 

analysis, objective, Cressmann, Gilchrist, 309. 

Arons, A. B.: Photographic Investigation of the 
Projection of Droplets by Bubbles Bursting at 
a Water Surface, 1. 

ARRHENIUS, G.: Origin and Accumulation of 
Aluminosilicates in the Ocean, 215. 

ARRHENIUS, O.: Chemical Denudation in Sweden, 
326. 

Atlantic Ocean, carbon dioxide exchange, Nelson- 
Dingle, 342. 

— —, drift currents, H. Stommel, 203. 

— -, salinity, F. C. Fuglister, 46. 

— —, temperature, F. C. Fuglister, 46. 

atmospheric chemistry, E. Eriksson, 302. 


Baltic Ice-Lake, B. Kullenberg, 221. 

baroclinic waves, N. A. Phillips, 273. 

barotropic model, Staff Members, Institute of Mete- 
orology, University of Stockholm, 139. 

Batu, M.: A Study of T Phases Recorded at the 
Kiruna Seismograph Station, 63. 

—:The Density Ratio at the Boundary of the Earth’s 
Core, 408. 

BERKOFSKY, L.: The Exact Determination of the 
Effective Domain of Dependence of a One- 
Dimensional Numerical Prediction Formula, 165. 

BLANCHARD, D. C.: Photographic Investigation of 
the Projection of Droplets by Bubbles Bursting 
at a Water Surface, 1. 

Bouman, D. J.: A Theory of the ventilated Psychro- 
meter based on the Thermodynamics of irre- 
versible Processes, 399. 

BRUNBERG, E. Ä.: On the Interpretation of the 
Diurnal Variation of Cosmic Rays, 73. 


BRUNBERG E. A.: Variations of the Cosmic Ray 
Intensity During Magnetic Storms, 254. 

bubbles, at a water surface, C. F. Kientzler, A. B. 
Arons, D. C. Blanchard, A. H. Woodcock, 1. 

Burt, W. V.: Specific Scattering by Uniform 
Minerogenic Suspensions, 220. 


Carbon dioxide, A. Nelson-Dingle, 342. 

charge generation in thunderstorms, Müller-Hille- 
brand, 367. 

chemical denudation, Arrhenius, 326. 

chemistry, atmospheric, Eriksson, 302. 

circulation, general, Starr, 268. 

—, meridional, Phillips, 273. 

climatic fluctuation, Mather, 287. 

clouds, cumulonimbus, Starr-Malkus, Ronne, 351. 

—, lee-wave, Larsson, 124. 

computor, electronic, Staff Members, Institute of 
Meteorology, University of Stockholm, 139. 

cosmic radiation, diurnal variation, Brunberg, 
Dates: 

— —, during magnetic storm, Brunbere, Dattner, 

254. 

electric field, Lundquist, 260. 

, origin, Alfvén, 232. 

ray intensity, Brunberg, Dattner, 254. 

variations, Brunberg, Dattner, 73; Brun- 
berg, Dattner, 254; Lundquist, 260. 

CRESSMAN, G. P.: An Experiment in Objective 
Analysis, 309. 

cross-section, Gulf Stream, Worthington, 116. 

cumulonimbus clouds, Starr-Malkus, Ronne, 351. 

currents, drift, Stommel, 203. 

—, surface, Post, 59. 

cyclogenesis, Namias, 8. 


DATTNER, A.: On the Interpretation of the Diurnal 
Variation of Cosmic Rays, 73. 

—: Variations of the Cosmic Ray Intensity During 
Magnetic Storms, 254. 


Tellus VI (1954), 4 


INDEX 


DEIRMENDJIAN, D.: Global Radiation Resulting from 
Multiple Scattering in a Rayleigh Atmosphere, 
382. | 

density ratio, Bäth, 408. 

denudation, chemical, Arrhenius, 326. 

dependence, domain of, Berkofsky, 165. 

domain of dependence, Berkofsky, 165. 

Doporto, M.: The Cellular Solution of the Dynam- 
ical Equations and the Actual Perturbation of 
the Westerlies up to 16 km at Valentia, 32. 

drift current, Stommel, 203. 

droplets, Kientzler, Arons, Blanchard, Woodcock, r. 

Duvepai, T.: A Quantitative Prognostic Method 
for Thickness Charts Using Advective Tendencies, 
192. 


Earth’s core, Bath, 408. 

eclipse 1945, Jerlov, Olsson, Schuepp, 44. 

EGNER, H.: Composition of Atmospheric Precipita- 
tion in Sweden, 261. 

Ekman, V. W., in memoriam, 201. 

electric field, cosmic rays, Lundquist, 260. 

— —, vertical, Miiller-Hillebrand, 367. 

EzassEen, E.: Numerical Solutions of the Pertur- 
bation Equation for a Linear Flow, 183. 

EMANUELSSON, A.: Composition of Atmospheric 
Precipitation in Sweden, 261. 

energy transfer, Phillips, 273. 

Eriksson, E.: Composition of Atmospheric Precipi- 
ation in Sweden, 261. 

—: Report on an informal conference in atmospheric 
chemistry held at the Meteorological Institute, 
University of Stockholm, May 24—26, 1954, 302. 

Eriksson, T. O.: A Measurement of the Winds 
Aloft During a Cold Front Passage, 27. 


Flow, linear, Eliassen, 183. 

fluid, dynamics, Long, 91. 

—, stratified, Long, 97. 

flux, momentum, Starr, White, 180. 

—, sensible heat, White, 177. 

Fonsetius, S.: Amino Acids in Rainwater, 90. 

forecasting, numerical, Scorer, 23; Staff Members, 
Institute of Meteorology, University of Stock- 
holm, 139. 

front, cold, Eriksson, 27. 

FUGLISTER, F. C.: Average Temperature and Salinity 
at a Depth of 200 Meters in the North Atlantic, 


46. 


General circulation, Starr, 268. 
geoid, form of, Hirvonen, 84. 


Tellus VI (1954), 4 


417 


GILCHRIST, B.: An Experiment in Objective Analysis, 
309. 

graupel, Müller-Hillebrand, 367 

gravimetric method, Hirvonen, 84. 

Gulf Stream, Worthington, 116. 


Heat, sensible, White, 177. 
Hirvonen, R. A.: The Gravimetric Method for 
Determination of the Form of the Geoid, 84. 


Ice crystals, Müller-Hillebrand, 367. 
irreversible processes, Bouman, 399. 


Jämtland mountains, Larsson, 124. 

Jertov, N.: Measurements of Solar Radiation at 
Lövänger in Sweden during the Total Eclipse 
1945, 44. 

Jorpan, C. L.: On Syonos Initial Stages of Typho- 
genesis, 198. 


KIENTZLER, C. F.: Photographic Investigation of the 
Projection of Droplets by Bubbles Bursting at 
a Water Surface, 1. 

Kiruna seismographic station, Bäth, 63. 

KULLENBERG, B.: On the Presence of Sea Water in 
the Baltic Ice-Lake, 221. 


Larsson, L.: Observations of Lee Wave Clouds in 
the Jämtland Mountains, Sweden, 124. 

lee wave clouds, Larsson, 124. 

linear flow, Eliassen, 183. 

Long, R. R.: Some Aspects of the Flow of Stratified 
Fluids II, 97. 

Lunpauist, S.: A Note on the Variation of Cosmic 
Ray Intensity Produced by an Electric Field, 260. 


Magnetic storm, Brunberg, Dattner, 254. 

MATHER, J. R.: The Present Climatic Fluctuation and 
its Bearing on a Reconstruction of Pleistocene 
Climatic Conditions, 287. 

mean motions, Thompson, 150 

minerogenic suspensions, Burt, 229. 

momentum, flux of, Starr, White, 180. 

MÜLLER-HILLEBRAND, D.: Charge Generation in 
Thunderstorms by Collision of Ice Crystals with 
Graupel, Falling Through a Vertical Electric 
Field, 367. 


Namias, J.: Quasi-Periodic Cyclogenesis in Relation 
to the General Circulation, 8. 

NELSON-DINGLE, A.: The Carbon Dioxide Exchange 
Between the North Atlantic Ocean and the 
atmosphere, 342. 

normals, in numerical forecasting, Scorer, 23. 


418 


Objective analysis, Cressman, Gilchrist, 309. 
O1sson, H.: Measurements of Solar Radiation at 
Lövänger in Sweden during the Total Eclipse 


194$, 44. 
one-dimensional prediction formula, Berkofsky, 165. 


Pernix-Hair, Raab, 405. 

perturbations, Doporto, 32; Eliassen, 183. 

Pants, N. A.: Energy Transformations and 
Meridional Circulations Associated with Simple 
Baroclinic Waves in a Two-level, Quasi- 
geostrophic Model, 273. 

pleistocene, climatic conditions, Mather, 287. 

Post, L. A.: A Practical Method for the Prediction 
of the Surface Currents of the Ocean, 59. 

precipitation, composition, Emanuelsson, Eriksson, 
Egnér, 261. | 

prediction, currents, Post, 59. 

prognostic method, Duvedal, 192, 

production, of vorticity, Van Mieghem, 170. 

psychrometer, Bouman, 399. 


Quasi-Geostrophic flow, Phillips,273; Thompson, 3 19. 


Raas, L.: The Finnish Radiosonde With Pernix- 
Hair, 405. 

radiation, global, Deirmendjian, Sckera, 382. 

— solar, Jerlov, Olsson, Schuepp, 44. 

radiosonde, Finnish, Raab, 405. 

RAETHJEN, P.: Dynamik der Zyklonen, book- 
review, 93. 

rainwater, content of, Fonselius, 90; Emanuelsson, 
Eriksson, Egnér, 261. 

Rayleigh atmosphere, Deirmendjian, Sekera, 382. 

Ronne, C: On the Structure of some Cumulo- 
nimbus Clouds which Penetrated the High 
Tropical Troposphere, 351. 


Salinity, North Atlantic, Fuglister, 46. 

scattering, Burt, 229. 

— multiple, Deirmendjian, Sekera, 382. 

SCORER, R. S.: The Use of Normals in Numerical 
Forecasting, 23. 

SEKERA, Z.: Global Radiation Resulting from 
Multiple Scattering in a Rayleigh Atmosphere, 
382. 

STAFF MEMBERS, INSTITUTE OF METEOROLOGY, 
UNIVERSITY OF STOCKHOLM: Results of Fore- 
casting with the Barotropic Model on an Elec- 
tronic Computor (BESK), 139. 


INDEX 


Srarr, V. P.: Commentaries Concerning Research 
on the General Circulation, 268. 

Starr, V. P., Warte, R. M.: Two Years of Mo- 
mentum Flux Data for 13°N, 180. 

SrarR—Matxus, J.: On the Structure of Some 
Cumulonimbus Clouds which Penetrated the 
High Tropical Troposphere, 351. 

STOMMEL, H.: Serial Observations of Drift Currents 
in the Central North Atlantic Ocean, 203. 

stratosphere, flux of heat, White, 177. 

suspensions, minerogenic, Burt, 229. 

symposium on the use of models in geophysical 
fluid dynamics, Long, 91. 

Syono, S.: On Syonos Initial Stages of Ttyphoge- 
nesis; reply, 199. 


Temperature, North Atlantic, Fuglister, 46. 
tendencies, advective, Duvedal, 192. 
thermodynamics, Bouman, 399. 

thickness charts, Duvedal, 192. 

Tuompson, P. D.: Prognostic Equations for the 
Mean Motions of Simple Fluid Systems and 
Their Relation to the Theory of Large Scale 
Atmospheric Turbulence, 150. 

— Note on the Integration of the Vorticity Equa- 
tion for Quasi-Geostrophic Flow, 319. 

thunderstorm, Müller-Hillebrand, 367. 

T-phases, Bäth, 63. 

tropical troposphere, Starr—Malkus, Ronne, 351. 

turbulence, Thompson, 150. 

typhogenesis, Jordan, Syono, 199. 


Waves, baroclinic, Phillips, 273. 

westerlies, Doporto, 32. 

Ware, R. M.: The Counter-Gradient Flux of 
Sensible Heat in the Lower Stratosphere, 177. 

— Two Years of Momentum Flux Data for 13°N, 
180. 

winds aloft, Eriksson, 27. 

Woopcock, A. H.: Photographic Investigation of 
the Projection of Droplets by Bubbles Bursting 
at a Water Surface, 1. 

WORTHINGTON, L. V.: Three Detailed Cross- 
Sections of the Gulf Stream, 116. 

Valentia, Doporto, 32. 

VAN MIEGHEM, J.: Transport and Production of 
Vorticity in the Atmosphere, 170. 

vorticity equation, Thompson, 319. 

— production, Van Mieghem, 170. 


Tellus VI (1954), 4 


